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1.  SUMARY 


Growth  rate  data  for  small  fatigue  cracks  in  fine  grained  (12  urn)  and 
coarse  grained  (30  um)  Astroloy  and  in  coarse  grained  (42  ym)  Waspaloy  were 
obtained  at  temperatures  of  25°C,  200°C  and  600°C.  Expressing  the  crack 
growth  kinetics  in  terms  of  the  linear  elastic  parameter  aK  resulted  in  small 
cracks  growing  faster  than  large  cracks  in  all  cases  except  that  of  the  fine 
grained  Astroloy  at  elevated  temperature.  This  difference  occurred  even 
though  the  crystallographic  nature  of  cracking  in  small  cracks  was  shown  to  be 
the  same  as  that  in  large  cracks,  based  on  metallographic,  selected  area 
electron  channeling,  and  fractographic  information.  Furthermore,  the  fracture 
morphology  of  small  and  large  cracks  also  exhibited  the  same  dependence  on 
temperature  and  applied  stress  level  indicating  that  the  crack  propagation 
mechanism  was  the  same  in  both  cases.  Consequently,  it  was  concluded  that  the 
observed  differences  in  kinetics  were  due  to  differences  in  mechanical  driving 
force  for  small  and  large  fatigue  cracks. 

Contributions  to  these  differences  in  mechanical  driving  force  were 
examined  using  high  resolution  scanning  electron  microscopy,  coupled  with  the 
stereoimaging  technique.  This  technique  enabled  the  quantification  of  local 
crack-tip  displacements,  strains,  and  crack  opening  loads  during  fatigue  crack 
growth.  Measurements  showed  that  the  crack-tip  strain  fields,  as  well  as  the 
nature  of  crack  closure,  were  significantly  different  for  small  and  large 
cracks.  These  observations  were  consistent  with  the  fact  that  the  linear 
elastic  fracture  mechanics  concept  of  small  scale  yielding  was  shown  to  be 
invalid  for  small  cracks  based  on  the  fact  that  measured  crack-tip  plastic 
zone  sizes  were  on  the  order  of  the  crack  size. 

In  view  of  this  situation,  the  local  crack-tip  measurements  were  used  to 
define  an  elastic-plastic  crack  driving  force  based  on  the  relationship  be¬ 
tween  crack-tip  opening  displacement  and  aJ,  the  cyclic  J-integral.  To 
facilitate  comparison  of  small  and  large  crack  data,  aJ  was  expressed  in  terms 
of  an  equivalent  aK,  termed  aKgq.  The  value  of  aKg^  for  small  cracks  was 
found  to  depend  on  significant  contributions  from  plasticity,  mixed  modeness 
and  crack  closure,  while  that  for  large  cracks  was  found  to  depend  only  on 


crack  closure.  Expressing  the  crack  growth  rates  of  small  and  large  cracks  as 
a  function  of  AKeq  resulted  In  kinetics  which  were  independent  of  crack  size, 
within  experimental  scatter. 

The  above  results  formed  the  basis  of  a  phenomenological  model  for  the 
growth  of  both  small  and  large  fatigue  cracks  in  which  aKeq  is  related  to  the 
applied  aK.  Since  AKgq  cannot  currently  be  computed  analytically,  the 
phenomenological  model  was  used  to  predict  the  fatigue  life  of  unflawed 
turbine  engine  discs.  The  initial  flaw  sizes  for  these  calculations  were  set 
equal  to  the  material  grain  size  —  a  size  typical  of  the  intrinsic  micro- 
structural  defects  expected  to  be  found  in  Ni-based  superalloys.  The  fatigue 
lives  predicted  using  the  phenomenological  model  were  found  to  be  in  better 
overall  agreement  with  measured  fatigue  lives  than  those  predicted  based  on 
the  classical  linear  elastic  fracture  mechanics  approach  and  large  crack 
data.  At  longer  lives,  the  latter  approach  overpredicted  the  measured  lives 
by  several  orders  of  magnitude.  This  overprediction  was  shown  to  be 
associated  with  the  the  occurrence  of  the  fatigue  crack  growth  threshold  for 
large  cracks,  AKth.  Based  on  these  results,  it  is  concluded  that  AKth  and 
associated  near-threshold  growth  rate  data  are  not  applicable  to  material 
selection,  design  and  damage  tolerant  analyses  of  turbine  engine  discs. 


2.  INTRODUCTION 


The  approach  to  insuring  the  safety,  reliability  and  cost  effectiveness 
of  high-performance  gas  turbine  engines  has  undergone  a  change  from  the  "safe 
life"  philosophy,  based  on  traditional  fatigue  crack  initiation  concepts,  to 
the  "damage  tolerant"  philosophy,  based  on  fracture  mechanics  concepts.  The 
latter  has  been  formally  implemented  by  the  United  States  Air  Force  in  two 
separate  programs.  The  Retirement-for-Cause  of  Turbine  Engine  Components 
Program  uses  fracture  mechanics  concepts  in  the  life  management  of  critical 
components  of  the  F100  engine  (2.1,  2.2].  The  Engine  Structural  Integrity 
Program  (ENSIP)  is  broader  in  scope  and  establishes  a  general  damage  tolerance 
design  specification  for  critical  components  of  future  engines  [2.3-2.51.  In 
both  programs,  crack-like  defects  must  be  assumed  to  exist  in  the  components 
at  the  time  of  production  and  analysis  and  testing  is  required  to  demonstrate 
that  these  defects  will  not  grow  to  a  critical  size,  corresponding  to 
catastrophic  failure,  within  the  lifetime  of  the  engine. 

The  above  approach  is  based  on  the  assumption  that  there  exists  a 
unique,  geometry- Independent  relationship  between  the  crack-tip  stress  inten¬ 
sity  factor,  defined  by  linear  elasticity,  and  the  race  of  fatigue  crack  pro¬ 
pagation.  Such  a  relationship  has  been  demonstrated  for  large  through¬ 
thickness  cracks  which  have  traditionally  been  tested  in  the  laboratory  (2.6- 
2.8(.  However,  more  recent  results  on  cracks  which  are  small  in  relation  to 
either  the  crack-tip  plastic  zone  size  or  the  material  microstructure  appear 
to  the  violate  such  a  relationship  [2.9-2.12].  These  "small"  cracks  appear  to 
grow  at  faster  rates  than  those  of  large  cracks  subject  to  the  same  elastic 
stress  intensity  factor  range,  aK.  Furthermore,  small  cracks  can  propagate 
at  aK  values  which  are  below  the  threshold  aK  value  exhibited  by  large  cracks. 

Unfortunately,  the  demarcation  between  "small"  and  "large"  crack 
behavior  is  not  well  defined.  Thus,  the  precise  impact  of  this  anomalous 
behavior  on  life  prediction  has  been  uncertain.  Nevertheless,  it  is  clear 
that  with  the  advent  of  improved  nondestructive  examination  techniques  this 
issue  will  increase  in  importantance.  For  example,  recent  analyses  have 
estimated  that  an  Improvement  in  NOE  resolution  of  about  a  factor  of  two  will 


result  in  an  order  of  magnitude  difference  in  the  life  predicted  using  crack 
growth  rates  from  "small1'  versus  "large"  cracks  [2.10]. 

The  above  situation  provided  the  impetus  for  the  current  program.  The 
primary  objective  of  this  program  was  to  define  the  limitations  of  linear 
elastic  fracture  mechanics  to  treat  the  growth  of  small  cracks  and  to  identify 
the  underlying  physical  mechanisms  responsible  for  the  anomalously  high  growth 
rates  exhibited  by  small  cracks.  The  anticipated  benefits  from  such  an  effort 
include  an  improved  basis  for  life  prediction,  as  well  as  more  realistic 
material  selection  criteria. 
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3.  MATERIAL  MICROSTRUCTURAL  AND  MECHANICAL  PROPERTIES 


Two  nickel-base  superalloys  typical  of  those  used  in  gas  turbine  engine 
discs  were  chosen  for  this  study:  1)  Astroloy,  examined  in  fine  grained  (FG) 
and  coarse  grained  (CG)  variants,  2)  and  Waspaloy.  The  thermo-mechanical 
processing  of  these  alloys  and  resulting  monotonic  and  cyclic  stress-strain 
properties,  as  well  as  creep  properties,  are  provided  in  this  section. 

3.1  Microstructural  Characterization 

The  Waspaloy  used  in  this  program  was  in  the  form  of  a  large  diameter 
turbine  disc  and  was  produced  by  Cameron  Iron  Works  using  a  cast-wrought 
route.  The  measured  chemical  composition  is  given  in  Table  3-1.  The  alloy 
was  melted  by  a  Vacuum  Induction  Melt/Vacuum  Arc  Refine  (VIM-VAR)  process  and 
heat  treated  as  follows: 


Solution 

2  hours 

1016C 

Oil 

Quench 

1st  Age 

4  hours 

850C 

Air 

Cool 

2nd  Age 

16  hours 

760C 

Air 

Cool 

The  Waspaloy  grain  size,  determined  in  accordance  with  ASTM  E112,  was 
42  urn  (ASTM  5.9).  The  microstructure,  which  consisted  of  equiaxed  grains, 
with  occasional  large  grains  up  to  200  ym  by  100  urn,  is  shown  in  Figure  3-1. 
The  y'  structure  was  also  examined,  and  was  a  bimodal  distribution  of 
spheroidal  particles,  typically  0.3  and  0.03  ym  in  diameter  (Figure  3-2). 

Both  the  FG  and  CG  Astroloy  were  produced  by  a  powder  metallurgy  route 
involving  argon  atomization,  powder  classification  and  blending  followed  by 
Hot  Isostatic  Pressing  (HIP)  into  billet  prior  to  forging.  The  powder  was 
supplied  and  HIP'ed  by  Wiggin  Alloys.  The  FG  discs  were  forged  by  Daniel 
Doncasters,  while  the  CG  discs  were  forged  by  Landish.  The  measured  chemical 
composition  is  given  in  Table  3-2.  All  discs  were  heat  treated  as  follows: 
I218°C  HiP  (30,000  psi) 

Solution  4  hours  1095°C  Air  Cool 

1st  Age  24  hours  650°C  Air  Cool 

2nd  Age  8  hours  760°C  Air  Cool 


TABLE  3-1 


CHEMICAL  COMPOSITION  OF  WASPALOY 


Specified  Composition 
Height  %  (ppm) 


Element 

Min 

Max 

Analysis 

Carbon 

0.02 

0.10 

0.023 

Si  1  icon 

- 

0.15 

0.037 

Manganese 

- 

0.10 

0.02 

Phosphorus 

- 

0.015 

0.003 

Sulphur 

- 

0.008 

<0.002 

Silver 

- 

(5) 

- 

Aluminum 

1.2 

1.6 

1.5 

Boron 

0.003 

0.01 

0.007 

Bismuth 

- 

(1) 

- 

Cobalt 

12.0 

15.0 

12.26 

Chromium 

18.0 

21.0 

19.52 

Copper 

- 

0.10 

- 

Iron 

- 

2.0 

0.97 

Molybdenum 

3.5 

5.0 

3.82 

Lead 

- 

(10) 

- 

Titanium 

2.8 

3.3 

3.29 

Zirconium 

0.02 

0.08 

0.062 

Nickel 

Remainder 

TABLE  3-2 


CHEMICAL  COMPOSITION  OF  ASTROLOY 


Specified  Composition 

Analysis 

Element 

Weight  %  (ppm) 

A3770 

A3771 

AAH2 

MIN 

MAX 

Carbon 

0.01 

0.035 

0.023 

0.025 

0.026 

Silicon 

- 

0.2 

0.016 

0.014 

0.021 

Manganese 

- 

0.15 

- 

0.010 

0.020 

Phosphorus 

- 

0.15 

0.003 

0.002 

0.003 

Sulphur 

- 

0.005 

<0.002 

<0.002 

<0.002 

Argon 

- 

(3) 

- 

- 

- 

Si  1 ver 

- 

(5) 

- 

- 

- 

Aluminum 

3.85 

4.15 

4.06 

4.09 

3.97 

Boron 

0.015 

0.03 

0.02 

0.02 

0.018 

Bismuth 

- 

(0.5) 

- 

- 

- 

Cobalt 

16 

18 

17.23 

17.12 

17.13 

Chromium 

14 

16 

14.92 

14.91 

14.84 

Copper 

- 

0.1 

- 

- 

0.02 

Iron 

- 

0.5 

0.07 

0.07 

0.05 

Molybdenum 

4.5 

5.5 

5.05 

5.05 

5.07 

Nitrogen 

- 

(50) 

- 

- 

- 

Oxygen 

- 

(150) 

- 

- 

- 

Lead 

- 

(10) 

- 

- 

- 

Titanium 

3.35 

3.65 

3.49 

3.53 

3.53 

Tungsten 

- 

0.05 

0.03 

0.02 

0.C2 

Zirconium 

Nickel 

0.03 

0.06 

Balance 

0.042 

0.040 

0.047 

The  CG  Astroloy  received  an  additional  grain  growth  treatment  of  4  hours 
at  1180°C,  before  being  re-solutioned  and  re-aged. 

The  grain  sizes  of  the  two  Astroloy  variants  were  12  ym  (ASTM  10)  for 
the  FG  material  and  30  ym  (ASTM  7)  for  the  CG  material.  The  general  micro¬ 
structure  and  y'  morphology  for  the  FG  Astroloy  are  shown  in  Figures  3-3  and 
3-4,  respectively.  Analogous  results  for  the  CG  Astroloy  are  shown  in 
Figures  3-5  and  3-6.  In  the  FG  material,  the  y'  existed  as  large,  monolithic 
particles  carried  over  from  forging,  (that  is,  not  broken  up  during  solution) 
with  a  fine  back-ground  dispersion  of  cubic  y',  approximately  0.2  ym  on 
edge.  In  the  CG  variant  the  monolithic  y'  was  broken  down  by  the  additional 
grain  growth  treatment.  The  structure  was  blmodal  and  cubic,  the  particles 
being  0.5  ym  and  0.1  ym  on  edge. 


3.2  Mechanical  Properties 

Tensile  tests  were  performed  under  both  monotonic  and  cyclic 
conditions.  The  specimens  utilized  are  shown  in  Figures  3-7  and  3-8.  The 
monotonic  tests  were  conducted  in  accord  with  BS4A4  Part  1  Section  3,  while 
the  cyclic  testing  was  conducted  in  accord  with  the  method  of  Pickard  and 
Knott  [3.1]. 

Tabulated  summaries  of  the  tensile  results  at  25°C,  200"C  and  600°C  are 
given  in  Tables  3-3  thru  3-5,  with  graphic  representation  given  in  Figures  3-9 
thru  3-11.  A  tabulation  of  the  strain  hardening  data  is  given  in  Tables  3-6 
and  3-7.  The  values  of  K  and  n  in  Table  3-4  were  determined  by  least  squares 
fitting  of  the  monotonic  stress-strain  data  to  the  following  equation: 


0 


n 


where  a  =  the  true  stress  and  ep  =  the  true  (logrithmic)  plastic  strain.  The 
cyclic  values  given  in  Table  3-7  were  similarly  obtained  using  the  following 
equation: 


Ae_ 

=  k 


Aq 

2 


n 


FIGURE  3-3.  Microstructure  of  FG  Astroloy 


FIGURE  3-4.  Y  Structure  of  FG  Astroloy 
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FIGURE  3-6.  Y  '  stnirti 


STRESS-STRAIN  PROPERTIES  OF  HASPALOY  AT  VARIOUS  TEMPERATURES 
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STRESS-STRAIN  PROPERTIES  OF  F6  ASTROLOY  AT  VARIOUS  TEMPERATURES 
(M  =  Monotonic  and  C  =  Cyclic) 
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STRESS-STRAIN  PROPERTIES  OF  CG  ASTROLOY  AT  VARIOUS  TEMPERATURES 
(M  =  Monotonic  and  C  =  Cyclic) 
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FIGURE  3-10.  Monotonic  (solid  lines)  and  cyclic  (dashed 
lines)  stress-strain  curves  for  FG  Astrolov 
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FIGURE  3-11.  Monotonic  (solid  lines)  and  cyclic  (dashed 
lines)  stress-strain  curves  for  CG  Astroloy 
at  several  temperatures. 


TABLE  3-6 


NONOTONIC  STRAIN  HARDENING  DATA 


Material 

Temp 

E 

k 

n 

Correlation 

Coefficient 

Waspaloy 

20 

211 

1207.8 

0.043 

0.98 

200 

215 

1051.4 

0.033 

0.99 

600 

181 

1099.2 

0.053 

0.99 

FG  Astroloy 

20 

214 

1162.3 

0.018 

0.92 

200 

214 

1218.1 

0.038 

0.94 

600 

185 

1231.6 

0.036 

0.98 

CG  Astroloy 

20 

222 

1191.1 

0.042 

0.99 

200 

205 

1068.5 

0.033 

0.99 

600 

168 

1245.8 

0.059 

0.99 

TABLE  3-7 


CYCLIC  STRAIN  HARDENING  DATA 


Material 

Waspaloy 


FG  Astroloy 
( A3771 ) 


FG  Astroloy 
( A3770 ) 


Temp 

E 

k' 

n* 

Correlation 

Coefficient 

20 

229 

1375.3 

0.064 

0.99 

200 

215 

1420.6 

0.075 

0.98 

600 

197 

963.7 

0.029 

0.69 

20 

216 

1372.8 

0.052 

0.96 

200 

207 

1591.5 

0.065 

0.98 

600 

192 

2604.1 

0.150 

0.95 

20 

225 

1464.8 

0.064 

0.97 

200 

201 

1539.2 

0.061 

0.98 

600 

198 

2423.4 

0.140 

0.96 

20 

220 

1274.3 

0.059 

0.99 

200 

221 

1642.0 

0.089 

0.99 

600 

183 

1480.0 

0.094 

0.94 

CG  Astroloy 


Constant  load  creep  strain  to  rupture  tests  were  conducted  at  600°C  for 
all  three  materials  to  define  the  stress  rupture  curve  between  50  and  500 
hours.  The  test  specimen  utilized  is  shown  in  Figure  3-12.  The  plotted 
rupture  data  are  shown  in  Figure  3-13. 

The  tensile  and  creep  properties  of  the  Waspaloy  and  the  FG  Astroloy 
were  considered  to  be  typical  for  materials  form  forgings  of  these  sizes.  The 
CG  Astroloy  exhibited  tensile  and  creep  properties  which  were  somewhat  below 
typical  levels,  but  nevertheless  above  the  specification  minimum  levels. 
Thus,  the  three  forgings  supplied  form  a  solid  basis  for  the  current  study. 

3.3.  References 

3.1  A.  C.  Pickard  and  J.  F.  Knott,  "Effect  of  Testing  Method  on  Cyclic 
Hardening  Behavior  of  Face  Centered  Cubic  Alloys,"  ASTM  Symposium  on  Low 
Cycle  Fatigue,  October  1985,  in  press. 
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FIGURE  3-12.  Creep  Rupture  Tes 
(dimensions  in  mm 


stress  (MPa  ) 


FIGURE  3-13.  Stress-rupture  data  for  Waspaloy,  and  FG  and  CG  Astroloy 
at  600°C . 


4.  FATIGUE  THRESHOLD  AND  LARGE  CRACK  GROWTH  BEHAVIOR 


4.1  Large  Crack  Testing  Procedure 

The  experiments  were  performed  to  a  Rolls-Royce  procedural  specification 
which  is  broadly  similar  to  the  ASTM  draft  specification  for  the  determination 
of  AKtn.  This  load  shedding  technique  utilizes  load  reductions  of  5%,  growing 
the  crack  at  least  three  times  the  calculated  maximum  plastic  zone  size  prior 
to  further  reductions  in  load.  As  the  threshold  is  approached,  the  load  re¬ 
ductions  are  decreased  in  magnitude  until  they  reach  the  limit  of  the  machines 
sensitivity.  The  threshold  condition  was  defined  as  no  detectable  growth  in 
107  cycles.  After  the  threshold  is  established,  the  alternating  load  is 
increased  at  constant  load  ratio  to  that  expected  to  yield  a  growth  rate  of 
10-7  mm/cycle  and  the  crack  allowed  to  grow  freely  under  constant  load 
ampl itude. 

Standard  compact-type  specimens  having  an  effective  width  of  26  mm  and  a 
thickness  of  13  mm  were  used  to  obtain  the  large  crack  data.  Threshold  and 
crack  growth  testing  was  performed  on  a  100  kN  capacity  Amsler  vibrophore 
electro-resonant  fatigue  machine  at  a  cyclic  frequency  of  130  Hz.  The  crack 
length  was  determined  by  a  O.C.  potential  drop  technique  as  described  by  Hicks 
and  Pickard  {4.1].  Crack  length  was  also  determined  optically  during  the  test 
as  a  back  up  to  that  determined  by  the  potential  drop  method. 

After  termination  of  the  tests,  the  specimens  were  heat  tinted  (where 
required)  and  the  final  crack  length  measured  at  five  points  across  the  crack 
front.  The  average  final  crack  length,  and  the  final  crack  voltage  were  used 
to  calibrate  the  potential  drop  to  a  crack  length  polynommal  function,  en¬ 
abling  an  accurate  post-test  determination  of  the  threshold  and  subsequent 
crack  growth  rates.  Crack  growth  data  were  not  determined  for  the  decreasing 
aK  phase  of  the  tests. 

Crack  growth  rates  were  determined  using  a  three  point  secant  method  as 
fol lows: 


da 

dN. 
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a. 


i  +  1 


N  . 


i  +  1 


a 

¥ 


i-1 

i-1 


[he  corresponding  value  of  aK  was  computed  using  a^. 


4.2  Near-Threshold  Fatigue  Crack  Growth  Rates  In  Large  Cracks 


The  effects  of  temperature  on  threshold  and  near  threshold  growth  rates 
are  summarized  in  Figures  4-1  through  4-3.  In  the  FG  Astroloy  material,  the 
threshold  conditions  is  independent  of  temperature,  but  with  fatigue  crack 
growth  rates  at  higher  stress  intensities  becoming  faster  with  increased 
temperature.  This  latter  feature  is  considered  to  be  a  manifestation  of  oxide 
induced  closure  [4.2] .  In  CG  Astroloy  (Figure  4-2)  a  trend  of  increasing 
threshold  with  decreasing  temperature  was  noted.  In  a  comparison  of  the 
coarse  and  fine  grain  variants  of  Astroloy  it  is  apparent  that  the  threshold 
values  at  600°C  are  quite  similar,  indicative  perhaps  of  an  equivalence  of 
surface  roughness.  The  effects  of  temperature  in  the  necklaced  variant  are 
shown  in  Figure  4-3.  In  the  absence  of  any  RT  data,  the  data  for  CG  Astroloy 
has  been  included.  A  comparison  of  the  necklace  data  with  that  of  the  FG  and 
CG  Astroloy  again  reveals  good  correlation  at  600C  with  increasing 
discrepancies  at  lower  temperatures.  The  variations  in  threshold  and  near 
threshold  crack  growth  with  microstructure  are  shown  in  Figures  4-4  and  4-5 
for  200°C  and  600°C,  respectively.  The  effects  of  load  ratio  at  600°C  are 
shown  in  Figures  4-6  and  4-7  for  FG  and  CG  Astroloy,  respectively.  At 
R  =  0.1,  thresholds  for  both  microstructural  variants  offer  comparatively  good 
agreement,  as  might  be  expected,  since  at  600°C  the  dominant  mode  of  material 
deformation  at  the  crack  tip  is  wavy  slip  producing  a  comparatively  smooth 
transgranular  fracture  surface  with  little  variation  in  fracture  surface 
roughness  between  microstructural  variants.  In  keeping  with  this,  the 
threshold  condition  will  be  dominated  by  oxide  induced  closure,  which 
manifests  itself  as  a  grain  size  independent  threshold.  The  high  load  ratio 
data  (R  =  0.5)  is  perhaps  a  more  representative  indication  of  the  true 
resistance  of  the  material  to  the  growth  of  a  fatigue  crack,  with  the  fine 
grain  material  being  marginally  superior. 

The  near  threshold  growth  rates  of  CG  Astroloy  in  air  and  vacuum,  are 
contrasted  in  Figure  4-8.  The  threshold  condition  in  vacuum  is  generally  re¬ 
cognized  as  being  higher  than  its  air  counterpart,  and  in  this  case  as  having 
a  significantly  higher  knee.  Work  by  Beevers  [4.3]  on  both  Waspaloy  and 
Astroloy  has  indicated  threshold  behavior  at  200°C  in  vacuum  to  be  variable, 
due  to  load  history  effects,  and  having  an  enhanced  tendency  towards  a  planar 
slip  deformation  mechanism.  The  behavior  at  600°C,  while  no  longer  dominated 
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FIGURE  4-1.  Temperature  Effects  in  Fine  Grain  Astroloy 
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.  Temperature  Effects  in  Coarse  Grain  Astrotoy 
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FIGURE  4-3. 


Effect  of  temperature  on  near  threshold  crack  growth 
behavior  of  necklace  Astroloy. 
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FIGURE  4-4.  Effect  of  microstructure  on  near  threshold  crack 
growth  behavior  of  Astroloy  (200°C). 
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FIGURE  4-5.  Effect  of  microstructure  on  near  threshold  crack 
growth  behavior  of  Astroloy  (600°C). 
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FIGURE  4-8.  Effect  of  environment  and  temperature  on  near  threshold 
crack  growth  rates  of  Coarse  Grain  Astroloy. 


by  such  a  tendency  towards  facetted  fracture  does  remain  similar  in  nature 
with  an  indication  of  a  higher  threshold  and  lower  growth  rates  in  the  Paris 
regime.  Although  the  range  of  growth  rates  explored  was  not  as  extensive  as 
the  200° C  vacuum  test,  corroboration  of  the  lower  growth  rates  in  vacuum  comes 
from  the  test  data  for  Waspaloy  at  55Q°C  !n  air  and  vacuum,  (Figure  4-9)  where 
growth  rates  in  air  are  a  factor  of  two  faster  than  the  vacuum  growth  rates. 

Crack  growth  data  for  Waspaloy  is  shown  in  Figure  4-10  along  with  the 
growth  rates  measured  for  coarse  grain  Astroloy  of  30  gm  grain  size.  The 
behavioral  trends  are  similar  to  those  seen  in  CG  Astroloy  (Figure  4-2), 
although  it  should  be  emphasized  that  the  threshold  conditions  measured  were 
lower  for  the  R1  test  than  the  600°C  test.  The  600°C  behavior  is  almost  iden¬ 
tical  to  the  6G0°C  behavior  of  the  30  m  grain  size  Astroloy,  with  threshold 
again  being  dominated  by  oxide  Induced  closure,  while  the  RT  qrowth  rates  are 
slightly  faster  than  the  30  urn  Astroloy  due  to  the  differences  in  yielding  be¬ 
havior,  and  the  fact  that  the  Waspaloy,  while  coarser  in  grain  size,  did  ..ot 
produce  facets  in  the  near  threshold  regime  to  the  same  extent  as  the  30  um 
grain  size  Astroloy. 

/ 

4.3  The  Influence  of  Grain  Size  on  Fatigue  Threshold 

The  room  temperature  CG  Astroloy  data  obtained  in  the  current  study  ex¬ 
hibits  significantly  lower  aKth  values  than  that  obtained  previously  for  a 
different  heat  of  CG  Astroloy  (4.4,  4.5j.  Several  factors  are  believed  to 
contribute  to  this  difference  as  follows:  First,  the  material  used  in  the 
previous  study  was  "As-HIP'ed",  and  particularly  prone  to  facetting  on  {111} 
planes  at  RT.  Secondly,  the  average  grain  size  in  the  previous  study  was 
50  um  as  opposed  to  30  um  in  the  current  study.  When  the  above  factors  are 
considered,  the  differences  in  the  measured  thresholds  are  not  unexpected. 
Figure  4-11  plots  the  measured  thresholds  from  this  and  other  studies  [4.4, 
4.6 1  of  Astroloy  against  grain  size.  The  correlation  between  points  is  very 
gcod,  and  certainly  puts  the  measured  values  into  perspective.  Although  the 
mechanism  by  which  grain  size  Influences  the  fatigue  crack  growth  rate  is  not 
completely  understood,  the  fractographlc  information  is  consistent  with  the 
operation  of  roughness-induced  and  oxide-induced  mechanisms  of  crack  closure. 
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FIGURE  4-9.  Effect  of  environment,  on  fat  * 
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4.4  Crack  Closure  in  Large  Cracks 


Crack  closure  results  are  of  interest  to  understanding  certain  trends 
exhibited  by  the  previous  crack  growth  rates  in  large  cracks.  In  addition, 
the  basic  differences  in  the  extent  of  crack  closure  developed  by  large  and 
small  cracks  need  to  be  taken  into  account  in  attempting  to  formulate  a 
driving  force  for  crack  growth  which  is  independent  of  crack  size. 
Fundamental  crack  closure  concepts  are  briefly  reviewed  below,  in  preparation 
for  a  summary  and  reinterpretation  of  literature  data  on  crack  closure  in 
Astroloy  large  cracks. 

From  the  basic  definition  of  crack  closure,  as  defined  originally  by 
El ber  ( 4. 7 | ,  we  have: 

AKeff  =  Kmax  ”  Kop  =  hK  ~  aKc1  (4_1) 

where:  A<cl  =  KQp  -  Kmjn 

where  the  definitions  of  the  above  parameters  is  given  in  Figure  4-12.  A1 
though  the  load  reguired  to  open  the  crack  (Pop)  is  the  basic  measurement  in 
any  crack  closure  study,  this  measurement  is  often  expressed  in  terms  of  any 
one  of  several  dimensionless  closure  parameters,  including:  Reff  =  P0p/Pmax  = 
K0p/Kmax,  or  AKgff/dK.  However,  based  on  the  definition  in  Eq.  (4-1)  these 
parameters  are  related  as  follows: 


AKeff/AK  =  (1  -  Reff)/(1  -  R)  (4-2) 

Recently,  Hudak  and  Davidson  [4.81  have  found  that  crack  closure  could 
be  described  by  the  following  relatively  simple  relationship  for  several  dif¬ 
ferent  materials  and  a  wide  range  of  fatigue  loading  variables: 

"•■effM  ‘  ‘  -  VW  «-3) 

where  KQ  is  a  material  constant.  Eq.  (4-3)  is  based  on  crack  closure  data 
obtained  using  the  stereoimaging  technique  which  enables  high  resolution  mea¬ 
surements  to  be  made  local  to  the  crack-tip  (see  Section  7.1).  This  technique 
allows  the  determination  of  the  load  required  to  fully  open  the  crack,  to 
within  a  fraction  of  a  micron  from  its  tip.  Nevertheless,  when  remote  tech¬ 
niques  are  used  to  measure  the  opening  load,  such  as  changes  in  specimen 
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FIGURE  4-12.  Definition  of  terms  used  to  characterize  fatigue 
crack  closure. 


compliance,  a  similar  relationship  is  also  found.  However,  the  latter 
techniques  often  detect  less  closure,  primarily  due  to  their  lack  of  precision 
with  regard  to  events  close  to  the  crack-tip. 

Combining  Eqs.  (4-1)  and  (4-3),  along  with  the  basic  definition  aK  = 
(1  -  R)Kmax.  gives: 


aKc1  =  K0(l  -  R)  (4-4) 

It  is  significant  to  note  that  the  above  observation  appears  to  hold  for  large 
cracks  subjected  to  a  wide  range  of  applied  aK  (or  Kmax)  and  R  values.  It 
should  also  be  noted  that  for  the  special  case  of  R  =  0  (and  Kmin  =0),  Eq  (4- 
4)  reduces  to: 


The  above  relationship  is  consistent  with  observations  on  a  variety  of 
materials  that  «op  is  relatively  constant  for  R  near  zero;  for  example,  see 
Ref.  [4.9]. 

Although  crack  closure  measurements  were  not  obtained  for  large  crack  in 
the  current  study,  available  literature  data  of  Soniak  and  Remy  on  necklace 
Astroloy  at  20"C  and  650°C  [4.10]  have  been  analyzed  using  Eqs.  (4-3)  and  (4- 
4).  These  results,  summarized  in  Figure  4-13,  indicate  relatively  good  agree¬ 
ment  with  the  above  equations,  although  considerable  scatter  is  evident  in  the 
elevated  temperature  data. 

The  data  in  Figure  4-13  show  measurably  less  crack  closure  (higher 
AKeff/AK  values)  at  650°C  than  at  25°C.  This  result  is  consistent  with  the 
crack  growth  rates  measured  in  the  present  study,  as  well  as  with  the  fracto- 
yraphy  in  Ni-base  alloys  which  shows  noticeably  flatter  surfaces  at  the 
elevated  temperature  (See  Section  6.3).  The  latter  observation  supports  the 
previous  hypothesis  of  Hicks  et  al  [4.11]  that  crack  closure  in  these  alloys 
is  controlled  by  fracture  surface  roughness.  Thus,  the  observed  temperature 
dependence  of  the  long  crack  growth  rates  appears  to  be  due  to  differences  in 
the  effective  driving  force  for  crack  growth,  as  well  as  to  expected  differ¬ 
ences  in  the  intrinsic  resistance  to  fatigue  crack  growth. 
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FIGURE  4-13.  Crack  closure  behavior  of  large  cracks  in  necklace 
Astroloy  at  various  temperatures,  from  Soniak  and 
Remy  [4.10]. 


The  intercepts  along  the  abscissa  In  Figure  4-13  correspond  to  aKc-j 
values  of  4-6  MPa/m  at  25°C  and  3-5  MPa/m  at  650°C.  These  results  are  dis¬ 
cussed  further  in  Section  7,  where  they  are  compared  with  crack  closure  mea¬ 
surements  on  small  cracks. 
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5.0  GROVfTH  OF  SHALL  CRACKS  AS  A  FUNCTION  OF  A K 


This  section  compares  small  crack  data  obtained  using  naturally  ini¬ 
tiated  surface  cracks  in  three-point  bend  specimens  with  the  large  crack  data 
of  Section  4.  These  results  are  also  compared  with  the  additional  small  crack 
data  from  Section  7  which  were  obtained  in  the  scanning  electron  microscope. 
This  initial  comparison  is  made  within  the  framework  of  linear  elastic 
fracture  mechanics,  although  it  is  recognized  that  this  approach  may  be 
inadequate  for  small  cracks.  Nevertheless,  it  does  provide  a  convenient  basis 
for  initial  comoarison  of  the  influence  of  material,  environmental  and  loading 
variables  in  a  fo.*mat  which  is  consistent  with  most  other  small  crack 
studies.  The  development  and  evaluation  of  a  more  appropriate  driving  force 
for  the  growth  of  small  cracks  is  presented  in  Sections  7  and  8.  Since  the 
success  of  any  small  crack  study  is  dependent  on  the  experimental  techniques 
employed,  particularly  at  elevated  temperature,  these  techniques  are  discussed 
below  prior  to  presenting  the  experimental  results. 

5.1  Small  Crack  Experimental  Procedure 

Data  on  the  growth  of  small  cracks  at  both  room  temperature  and  elevated 
temperature  were  obtained  using  10x10x70  mm  specimens  loaded  in  three  point 
bending  as  shown  in  Figure  5-1.  These  specimens  were  prepared  using  low- 
stress  machining.  In  addition,  to  eliminate  residual  stresses  which  may  have 
remained  after  machining,  the  gage  section  was  electropol ished  in  a  solution 
of  2 5  percent  nitric  acid  in  methanol  at  -30°C.  Typically  a  130  ym-thick 
layer  of  material  was  removed  during  electropolishing.  Induction  heating  was 
used  to  conduct  the  elevated  temperature  testing.  To  accommodate  the 
induction  heating  both  the  upper  and  lower  bend  fixtures  were  water  cooled. 
Contact  between  the  upper  fixture  and  the  specimen  was  achieved  using  a 
ceramic  tup  in  order  to  eliminate  localized  induced  heating  in  this  region  of 
the  specimen.  The  configuration  of  the  induction  heating  coils  was  also 
optimized  to  minimize  temperature  gradients  in  the  central  region  of  the 
specimen.  Thermal  transients  arising  from  room  air  currents  were  minimized  by 
placing  a  plexiglass  box  around  the  center  of  the  test  machine  as  shown  in 
Figure  5-1.  With  this  system,  the  temperature  in  the  region  of  the  specimen 
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(b)  Close-up  of  specimen  and  loading 
fixtures . 


(a)  Overview  of  test  system  and  induction  heater 


containing  the  small  cracks  was  constant  to  within  ±3°C  as  indicated  by 
results  in  Figure  5-2  from  thermocouple  measurements  obtained  for  a  prescribed 
temperature  of  650°C. 

All  cracks  were  naturally  initiated  from  the  electropol ished  surface  of 
the  specimen.  Generally  two  or  three  independent  cracks  occurred  and  were 
monitored  using  surface  replication  at  intervals  of  10,000  to  20,000  cycles. 
Replicas  were  subsequently  coated  with  vapor-deposited  silver  and  crack 
lengths  were  measured  using  a  light  microscope  at  magnifications  of  200X  to 
1000X.  Although  this  technique  involved  periodic  interruption  of  the  experi¬ 
ment,  and  cooling  in  the  case  of  the  elevated  temperature  experiments,  this 
inconvenience  was  outweighed  by  the  fact  that  reliable  data  were  obtained  on 
multiple  cracks  as  small  as  10  urn  in  length.  The  detection  of  the  smallest 
cracks  was  possible  due  to  the  fact  that  larger,  more  readily  detectable 
cracks  could  be  easily  found  and  back-tracked  on  previous  replicas.  The  fact 
that  cracks  naturally  initiated  at  either  voids  or  persistent  slip  bands  in 
favorably  oriented  grains  enabled  interactions  between  these  microstructural 
features  and  the  growth  of  small  cracks  to  be  studied. 

The  above  specimen  geometry  and  loading  mode  proved  to  have  the  follow¬ 
ing  advantages  for  small  crack  testing:  First,  due  to  the  stress  gradient  in 
the  bend  specimen,  cracks  initiated  at  the  surface,  rather  than  at  subsurface 
pores  and  inclusions;  whereas,  subsurface  crack  initiation  is  often  a  problem 
in  tensile  loading.  Secondly,  the  localization  of  the  maximum  bending  moment 
at  the  center  of  the  specimen  in  3-point  bending  facilitated  the  detection  and 
monitoring  of  microscopic  cracks  by  limiting  the  area  over  which  cracks  ini¬ 
tiated.  Finally,  the  simplicity  of  specimen  design  and  fixturing  is  both 
economical  and  promoted  efficient  heating  and  cooling  during  elevated  tempera¬ 
ture  experiments.  The  major  disadvantage  of  the  technique  was  the  fact  that 
the  average  crack  initiation  lives  were  longer  than  those  which  would  be  ex¬ 
pected  under  tensile  loading.  This  difference  is  believed  to  be  partly 
statistical  --  that  is,  due  to  the  fact  that  less  volume  of  highly  stressed 
material  occurs  in  the  bend  specimen,  thereby  reuucing  the  probability  of 
encountering  favorable  initiation  sites  in  the  region  of  highest  stress.  In 
addition,  at  high  applied  stresses  contained  plastic  deformation  near  the 
surface  of  the  bend  specimen  caused  compressive  stresses  to  be  generated  which 
reduced  the  local  mean  stress  and  tensile  stress  range  as  discussed  below. 
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5-2.  Temperature  profile  in  induction-heated  bend 
specimen  at  650°C. 


5.2  Stress  Intensity  Factor  for  Small  Surface  Cracks  In  Bending 


The  following  simplified  stress  intensity  factor  expression,  derived 
from  Pickard's  [5.1]  three-dimensional  finite  element  analyses  of  surface 
cracks  in  rectangular  bars,  was  used  to  analyze  the  growth  of  small  cracks  in 
the  bend  bar  experiments: 


aK  =  1.155  {1  -  0 . 5(a/W) }  {Smax(l-  R)/(l-  Ro)}  ✓wa  (5.1) 

where:  aK  =  the  stress  intensity  factor  range,  at  the  position  where  the 

crack  intersects  the  specimen  surface 

a  =  the  half-crack  length  (or  crack  depth  for  the  case  of 
a  semicircular  crack) 


W  =  the  half-height  of  the  specimen 

Smax  =  the  maximum  "outer-f 4ber"  stress  computed  for  three-point 
bending 

R  =  the  stress  ratio,  given  by  the  ratio  of  the  applied  minimum 
to  maximum  load 


R  =  the  local  stress  ratio  during  fatigue  loading,  given  by  the 

ratio  of  local  minimum  to  maximum  stress.  (As  discussed  below, 
this  ratio  can  differ  from  R) 


The  relative  simplicity  of  Eq.  (5.1)  compared  to  the  general  solution  given  in 
Ref.  [5.1]  arises  from  the  fact  in  the  current  experiments  the  relative  crack 
length  was  always  less  than  1  mm,  thus  a/W  <  0.2,  and  the  crack  aspect  ratio 
(a/c)  was  approximately  equal  to  one.  Under  these  conditions  Eq.  (5.1)  re¬ 
presents  the  original  finite  element  results  to  within  0.5  percent. 

The  assumption  that  the  crack  remained  semicircular  (a/c  =  1)  as  it  grew 
was  assessed  by  post-mortem  examination  of  the  fracture  surfaces.  The  deter¬ 
mination  of  the  crack  shape  was  relatively  easy  for  experiments  conducted  at 
elevated  temperature  since  the  fracture  surfaces  exhibited  a  color  gradient 
due  to  differences  in  oxide  thickness.  Results  from  these  microscopic 
examinations  are  summarized  in  Figure  5-3.  Although  considerable  scatter  in 
the  aspect  ratio  is  evident,  no  consistent  trend  emerged  with  respect  to  alloy 
composition  or  test  temperature,  indicating  that  the  nominal  crack  shape  is 
controlled  by  specimen  geometry  and  loading  mode.  Most  importantly,  the 
average  aspect  ratio  is  consistent  with  the  assumed  crack  shape  used  in  the 
analysis. 
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All  experiments  were  conducted  such  that  the  bend  specimen  was  nominally 
elastic.  However,  due  to  the  stress  gradient  in  the  specimen,  localized 
plastic  deformation  sometimes  occurred  in  the  surface  region  of  the  specimen 
where  cracks  initiated  and  grew.  (This  behavior  Is  analogous  to  the  localized 
deformation  which  can  occur  in  notched  geometries.)  The  primary  result  of 
this  contained  plasticity  at  the  surface  of  the  specimen  was  an  alteration  in 
the  local  mean  stress. 

The  above  process  was  quantified  using  strain  gage  measurements  on  an 
uncracked  CG  Astroloy  specimen  cycled  at  different  load  (strain)  ranges.  As 
shown  in  Figure  5-4,  the  local  stress  ratio  (R  )  decreased  and  became  negative 
as  the  applied  strain  range  increased.  The  values  of  Rq  in  the  CG  Astroloy 
bend  specimen  were  computed  from  the  measured  strains  using  the  monotonic  and 
stable  cyclic  stress-strain  relations  provided  in  Section  3.  Since  informa¬ 
tion  on  cyclic  stress  relaxation  could  not  be  included  in  these  computations, 
they  provide  an  upper-bound  estimate  of  the  local  stress  ratio  in  the  bend 
specimen.  The  additional  influence  of  cyclic  stress  relaxation  on  the  local 
mean  stress  can  be  estimated  from  the  measurements  of  Lindholm  et  al  (5.21, 
obtained  on  the  superalloy  B1900+Hf  tested  under  cyclic  strain-control.  Since 
uniaxial  tensile  specimens  were  used  in  this  study,  the  mean  stress  could  be 
measured  directly.  The  difference  between  these  results  and  those  obtained 
for  the  bend  bar  is  primarily  due  to  cyclic  stress  relaxation.  It  is  in¬ 
teresting  to  note  that  the  measured  results  are  In  good  agreement  with  values 
computed  from  the  Bodner-Partom  constitutive  relations  and  show  a  limiting 
case  of  fully-reversed  cycling  (R  =  -1). 

The  results  in  Figure  5-4  including  relaxation  were  assumed  to  apply  to 
all  of  the  Ni-base  superalloys  examined  in  the  current  program.  This 
assumption  is  consistent  with  the  fact  that  the  extent  of  mean  stress  relaxa¬ 
tion  is  known  to  primarily  depend  on  the  magnitude  of  the  cyclic  plastic 
strain  range  for  a  given  class  of  materials  (5.3).  The  value  of  strain  ex¬ 
perienced  in  the  bend  bars  for  various  materials  and  temperatures  was  assumed 
to  scale  with  material  strength  level  as  shown  in  Figure  5-5.  Thus,  for  a 
given  applied  normalized  stress  S  =  ^max^ys’  *-he  s^ra^n  was  obtained  from 
Figure  5-5,  and  the  corresponding  local  mean  stress  from  Figure  5-4.  Since 
the  crack  is  not  expected  to  grow  while  the  local  stress  is  compressive,  the 
local  stress  ratio  was  used  to  eliminate  the  local  compressive  stress  from  the 
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Figure  5-4.  Dependence  of  local  stress  ratio  (R0)  on  local  cyclic 
strain  range  (Ae)  for  three-point  bend  specimens. 
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computation  of  aK  by  reducing  the  applied  stress  by  the  factor  1/(1-  R  ) ,  as 
indicated  in  Eq.  5.1. 

5.3  Effect  of  Elastic  Anisotropy  and  Mixed  Hode  Loading  on  the  Stress 

Intensity  Factor  for  Small  Cracks 

As  shown  in  Section  6,  the  crystallographic  nature  of  the  cracking  en¬ 
countered  in  small  cracks  causes  cracks  to  be  inclined  at  an  angle  a  with 
respect  to  the  remote  loading  direction  --  for  example,  see  Table  6-1.  This 
inclined  crystal lographic  crack  growth  introduces  mixed-mode,  anisotropic 
effects  which  need  to  be  considered  in  analyzing  small  crack  data.  The  effect 
of  these  factors  on  the  stress  intensity  factor  for  small  cracks  was  assessed 
based  on  anisotropic  finite  element  (FE)  analyses  conducted  by  RR  under  the 
present  program  (see  Appendix),  as  well  as  the  anisotropic  boundary-integral- 
equation  (BIE)  results  of  Chan  and  Cruse  [5.4]. 

The  two  dimensional  BIE  results  of  Chan  and  Cruse  [  5 . 4 1  on  Mar-M200 
single  crystal  compact  specimens  are  shown  in  Figure  5-6.  Here  the  normalized 
stress  intensity  factor  is  given  as  a  function  of  the  projected  relative  crack 
size,  where:  a1  is  the  projected  crack  length,  b  is  the  specimen  thickness,  W 
is  the  effective  specimen  width,  and  P  is  the  applied  load.  The  similarity  of 
results  for  normal  cracks  (s  =  0°)*  with  two  crystal lographic  orientations 
(i  =  0’  and  45°)  indicates  that  the  effect  of  elastic  anisotropy  on  the  stress 
intensity  factor  is  insignif icant.  This  conclusion  is  also  supported  by  the 
fact  that  these  results  were  found  to  be  in  excellent  agreement  with  the 
elastic  isotropic  results  recommended  by  ASTM  for  the  compact  specimen 
geometry  [5.5|.  In  addition,  as  shown  in  Figure  5-6,  results  for  cracks  at 
various  inclined  crack  angles  3  were  also  found  to  be  similar  for  small  values 
of  s  and  small  relative  projected  crack  sizes.  These  results  suggest  that 
elastic  anisotropic  effects  can  be  ignored  for  conditions  encountered  in  the 
present  study  and  that  mixed  mode  effects  can  be  approximated  by  employing  the 
projected  crack  length  to  compute  the  dominant  Mode  I  stress  intensity  factor 
using  Eq.  5.1. 

Based  on  the  results  of  Chan  and  Cruse  [5.4|,  and  the  currently  observed 
crack  growth  inclinations,  it  is  estimated  that  the  stress  intensity  factor 

*  Note  that  0  in  Figure  5-6  and  0  in  Table  6-1  are  re'ated  as  follows: 

0  =  e  -  9o  . 
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incl ined 
i  (Ref.  5.4) 


for  Mode  II  will  always  be  less  than  the  35  percent  of  cne  Mode  I  value. 
Since  the  overall  elastic  driving  force  under  mixed  mode  loading  in  single 
crystals  Ni-based  alloys  has  been  shown  to  be  governed  by  the  strain  energy 
release  rate  15.6-5.7],  which  is  proportional  to  aKj2+aKjj2,  the  Mode  II 
contribution  will  be  less  than  6  percent  of  the  total  elastic  driving  force. 
This  contribution  has  been  neglected  in  analyzing  the  results  from  the  current 
study. 

An  obvious  limitation  of  the  above  results  is  the  fact  that  they  are  for 
long-thru  cracks  in  a  specimen  geometry  which  differs  from  those  used  in  the 
present  study.  Thus,  the  ultimate  objective  of  the  FE  analyses  conducted  by 
RR  under  the  present  program  was  to  examine  the  importance  of  material 
anisotropy  for  the  three-dimensional  case  of  surface  cracks.  The  approach 
taken  was  a  systematic  one  which  included  both  two-dimensional  and  three- 
dimensional  FE  modeling  of  cracks  in  anisotropic  bodies  and  of  cracks  in  local 
regions  of  anisotropy  w'*"hin  an  otherwise  isotropic  body.  Computational 
details  of  this  approach  are  provided  in  the  Appendix,  while  the  results  which 
are  significant  to  the  analysis  of  small  cracks  are  summarized  below. 

Results  of  the  finite  element  analyses,  for  several  different 
anisotropic  orientations,  revealed  an  increase  in  the  stress  intensity  factor 
compared  to  the  isotropic  case.  However,  as  shown  in  Figure  5-7,  on  average 
these  results  differ  py  only  5  percent.  This  same  figure  shows  that  this 
difference  is  small  compared  to  the  variation  in  stress  intensity  factor  which 
occurs  for  various  positions  along  the  semicircular  crack  front,  characterized 
by  the  angle  delta.  Thus,  the  issue  of  elastic  anisotropy  is  less  important 
than  that  of  selecting  the  proper  isotropic  stress  intensity  factor  along  the 
crack  front  to  characterize  the  crack  driving  force.  In  the  current  study, 
the  value  at  the  position  where  the  crack  intersects  the  specimen  surface  was 
used.  This  value  was  found  to  provide  the  best  overall  agreement  between 
fatigue  crack  growth  results  from  the  largest  surface  cracks  and  from  the  long 
thru  cracks. 

5.4  Effect  of  Temperature  on  the  Initiation  and  Early  Growth  of  Microcracks 

The  initiation  and  early  growth  of  microcracks  at  room  temperature  was 
relatively  easy  to  accomplish  in  the  three-point  bend  specimens.  For  example, 
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Figure  5-7.  Comparison  of  isotropic  and  anisotropic  Model  I  stress  intensity  factors  from 
finite  element  model  of  a  semicircular  surface  crack. 


when  the  maximum  (outer  fiber)  bending  stress  was  about  equal  to  the  materials 
yield  strength,  microcracks  of  15  to  20  um  in  total  surface  length  were 
obtained  in  about  105  cycles.  The  subsequent  growth  of  these  cracks  was  also 
relatively  easy  to  monitored  since  it  occurred  over  70  percent  to  85  percent 
of  the  total  cyclic  life  of  the  specimen.  However,  the  initiation  and  early 
growth  of  microcracks  at  elevated  temperature  was  considerably  more  difficult 
to  produce.  As  indicated  by  the  S-N  data  in  Figure  5-8,  applied  stresses  had 
to  be  increased  by  about  40  percent  In  order  to  produce  cyclic  lives 
comparable  to  those  experienced  at  room  temperature. 

In  order  to  facilitate  the  acquisition  of  data  at  elevated  temperature, 
microcracks  were  in  many  cases  initiated  at  room  temperature  before  increasing 
the  temperature  for  the  remainder  of  the  experiment.  However,  even  with  this 
procedure  the  growth  of  microcracks  was  erratic  at  elevated  temperature  in  all 
three  materials  examined.  An  example  of  this  erratic  behavior  is  given  by  the 
CG  Astroloy  data  summarized  in  Table  5-1.  Fatigue  precracking  at  room  temper¬ 
ature  was  conducted  at  S  =  1.21  and  R  =  0.1  for  160,000  cycles  to  insure  that 
growing  cracks  were  present.  Although  several  cracks  initiated  under  these 
conditions,  Table  5-1  only  provides  information  on  the  dominant  crack  which 
eventually  caused  the  specimen  to  fail.  This  particular  crack  initiated  from 
a  pore  on  the  electropol ished  surface  of  the  specimen.  Shortly  after  the 
temperature  was  increased  from  room  temperature  to  600°C,  the  crack  arrested 
and  remained  dormant  for  most  of  the  remaining  specimen  life.  During  this 
period  the  applied  stress  was  gradually  increased  in  steps.  Following  the 
final  increase  to  S  =  1.45,  the  crack  grew  rapidly  and  resulted  in  specimen 
failure  in  only  25,000  cycles.  Similar  behavior  was  also  observed  at  higher 
applied  stresses,  as  shown  in  Fiqure  5-9.  In  this  case,  the  crack  initiated 
in  only  30,000  cycles  due  to  the  high  stress  (S  =  1.8),  but  then  remained 
dormant  for  about  120,000  cycles  before  the  growth  rate  increased  several 
orders  of  magnitude,  resulting  in  failure  of  the  specimen. 

The  above  behavior  made  it  difficult  to  monitor  crack  growth  at  600°C 
using  the  replication  technique.  This  problem  was  eventually  overcome  after 
numerous  specimens  were  tested,  thereby  enabling  the  periods  of  initiation  and 
dormancy  to  be  estimated.  Using  this  experience,  the  replication  interval  was 
reduced  at  the  appropriate  time  so  that  adequate  data  could  be  obtained. 
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Figure  5-8.  Influence  of  relative  applied  stress  level 
(S  =  Smax/Sys)  and  temperature  on  the 
cyclic  life  of  10x10x70  mm  bend  specimens 
of  CG  Astroloy.  (NPC  =  non-propagating 
crack) 
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Figure  5-9.  Example  of  erratic  crack  growth  in  Ni-base  alloys 
at  600°C. 


The  cyclic  lives  obtained  using  the  three-point  bend  specimens 
(Figure  5-8)  are  longer  than  those  which  one  would  expect  in  gas  turbine 
discs.  This  difference  is  believed  to  be  due  to  differences  in  the  crack 
initiation  lives  between  small  laboratory  specimens  and  actual  components 
which  result  from:  1)  differences  in  surface  finish,  and  2)  a  statistical 
size  effect.  As  discussed  in  Refs.  1 5 . 8 J ,  actual  components  contain  a  variety 
of  manufacturing  defects  such  as  dents  and  scratches  which  often  serve  as 
mechanical  stresses  concentrations  for  crack  initiation.  Such  defects  are 
generally  not  present  in  laboratory  specimens  --  for  example,  the  specimens 
used  in  the  current  study  were  polished  both  mechanically  and 
electrochemical ly.  Furthermore,  the  relatively  small  size  and  mode  of  loading 
of  the  bend  specimens  resulted  in  a  limited  amount  of  highly  stress  material 
being  sampled.  Thus,  the  probability  of  this  highly  stressed  region 
containing  a  void  or  other  microstructural  defect  which  would  assist  in  the 
nucleation  of  a  crack  was  also  relatively  small.  Consequently,  many  of  the 
cracks  in  the  small  laboratory  specimens  initiated  at  persistent  slip  bands  in 
favorably  oriented  surface  grains,  rather  than  at  material  defects.  The  above 
arguments  are  consistent  with  the  data  in  Figure  5-10  which  compares  cyclic 
lives  of  the  three  point  bend  specimens  with  those  from  much  larger  tensile 
specimens  containing  simulated  bolt  holes  [5.9).  Note  that  the  lives  in  the 
latter  case  are  significantly  shorter.  Interestingly,  the  apparent  fatigue 
limit  in  the  bolt-hole  specimen  corresponds  to  the  stress  in  the  smooth 
specimen  which  produced  non-propagating  cracks  (NPC)  which  were  the  size  of 
one  or  two  grains. 

Although  the  precise  mechanism  responsible  for  the  long  periods  of  crack 
dormancy,  following  crack  initiation,  at  elevated  temperature  is  not  clearly 
defined,  several  mechanisms  can  be  postulated.  Initially,  the  possibility  of 
a  decrease  in  the  effective  crack  driving  force  due  to  oxide-induced  crack 
closure  [5.10,  5.11)  was  considered.  However,  the  importance  of  this 
mechanism  was  subseauentiv  discounted  when  the  same  erratic  crack  growth  be- 
havior  occurred  during  experiments  conducted  in  10  torr  vacuum  at  600T, 
where  oxidation  was  absent.  The  fractographic  observations  of  Section  6  show 
that  the  sharpness  of  the  crystallographic  facets  decrease  as  temperature  was 
increased  from  25°C  to  600°C,  thereby  suggesting  that  multiple  slip  becomes 
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Figure  5-10.  Comparison  of  cyclic  lives  of  notched  (kt  =  2.13) 
tension  specimens  and  smooth  (kt  =  1)  bend  speci¬ 
mens  in  Astroloy  at  600°C.  (NPC  =  non-propagating 
crack) 


elevated  temperature  since  the  fracture  surfaces  exhibited  a  color  gradient 
due  to  differences  in  oxide  thickness.  Results  from  these  microscopic 
examinations  are  summarized  in  Figure  5-3.  Although  considerable  scatter  in 
the  aspect  ratio  Is  evident,  no  consistent  trend  emerged  with  respect  to  alloy 
composition  or  test  temperature,  indicating  that  the  nominal  crack  shape  is 
controlled  by  specimen  geometry  and  loading  mode.  Most  importantly,  the 
average  aspect  ratio  is  consistent  with  the  assumed  crack  shape  used  in  the 
analysis. 
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activated  at  the  higher  temperature.  ft  was  also  observed,  through  periodic 
replication,  that  the  dormant  cracks  at  600°C  exhibited  a  time-dependent  in¬ 
crease  in  the  Mode  I  crack  opening  displacements.  The  resulting  crack-tip 
blunting  and  stress  relaxation  which  one  would  associate  with  this  process  is 
consistent  with  the  operation  of  multiple  slip  at  elevated  temperature. 
Although  details  of  this  process  remain  sketchy,  it  would  seem  to  be  related 
to  the  difficulty  encountered  in  growing  microcracks  at  elevated  temperature. 

5.5  Sources  of  Variability  In  Snail  Crack  Growth  Rates 

One  of  the  cnaracterlstics  of  small  cracks  is  their  large  variability  in 
crack  growth  rates,  relative  to  those  which  occur  in  long  cracks.  For  ex¬ 
ample,  variability  of  a  factor  of  ten,  or  more,  is  common  in  small  cracks, 
whereas  for  long  cracks  in  the  same  growth  rate  regime  a  factor  of  two  is 
typical  [5.12].  For  both  practical  and  fundamental  reason,  it  is  of  interest 
to  assess  whether  this  difference  arises  from  basic  differences  in  the  growth 
mechanisms  of  large  versus  small  cracks  or  from  differences  in  the  precision 
of  measurement  used  in  each  case. 

In  order  to  examine  the  above  issue,  as  well  as  validate  the  small  crack 
measurement  procedure  used  in  the  current  study,  a  series  of  replicate  experi¬ 
ments  was  performed  on  CG  Astroloy  at  room  temperature.  Information  was 
obtained  on  the  specimen-to-specimen,  lab-to-lab,  and  heat-to-heat  variabil¬ 
ity,  as  well  as  the  influence  of  crack  measurement  interval  and  data  analysis 
procedure  on  the  variability  in  small  crack  growth  rates. 

To  examine  the  reproducibility  of  data  from  lab-to-lab  several  small 
crack  experiments  were  conducted  at  both  SwRl  and  RR  using  the  same  test 
procedure  described  in  Section  5.1.  As  shown  in  Figure  5-11,  there  was  no 
discernable  difference  between  results  from  the  two  laboratories,  although  the 
variability  in  rates  was  in  both  cases  about  a  factor  of  ten.  Similar  agree¬ 
ment  is  found  between  the  current  RR  data  and  those  generated  previously  at  RR 
on  a  different  heat  of  CG  Astroloy  (5.13),  as  indicated  in  Figure  5-12.  For 
comparison,  Figure  5-12  also  gives  the  corresponding  long  crack  data  for  these 
same  two  heats  of  Astroloy.  As  discussed  previously  the  long  crack  data 
appear  to  be  dependent  on  grain  size,  whereas  the  small  crack  data  are  re¬ 
latively  independent  of  grain  size. 
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Figure  5-11.  Comparison  of  small  crack  growth  obtained 
at  SwRI  (□)  and  RP  (A). 
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As  shown  In  Figure  5-13,  the  variability  In  growth  rates  from  a  single 
small  crack  is  equal  to  that  exhibited  by  the  entire  data  set  consisting  of  7 
different  cracks  from  three  replicate  specimens.  Thus,  the  observed  variabil¬ 
ity  in  the  small  crack  data  is  not  due  to  a  lack  of  reproducibility  from 
specimen-to-specimen.  This  characteristic  variability  is  believed  to  be 
primarily  due  to  the  incremental  manner  in  which  cracks  grow  and  not  to  errors 
in  crack  length  measurement.  For  example,  cracks  often  appear  to  stop  growing 
for  tens  of  thousands  of  cycles,  even  when  high  resolution  measurement 
techniques  are  used.  Estimating  the  precision  of  crack  length  measurement 
with  the  current  replication  technique  to  be  1  um  gives  a  growth  rate  during 
these  apparent  dormant  periods  of  less  than  1X10"7  navVycIe. 

A  detailed  analysis  of  this  irregular  growth  was  performed  on  three 
specimens  of  CG  Astroloy  tested  at  20°C  in  order  to  examine  the  potential  role 
of  interactions  between  the  growing  crack  and  microstructural  features. 
Temporary  crack  arrests  occurred,  on  average,  each  95.00G  cycles  and  lasted 
for  20,000  to  30,000  cycles.  Results  from  one  of  the  three  specimens  which 
best  illustrates  the  variety  of  interactions  which  can  occur  are  shown  in 
Figure  5-14.  As  indicated,  although  arrests  were  sometimes  observed  at  either 
grain  boundaries  or  triple  points,  this  was  not  always  the  case.  In  fact, 
arrests  only  appeared  to  be  associated  with  such  interactions  about  25  percent 
of  the  time.  The  present  observations  also  showed  that  whether  or  not  the 
crack  changed  crystallographic  plane  as  it  crossed  the  grain  boundary  had  no 
bearing  on  the  ability  of  the  grain  boundary  to  arrest  the  crack.  Thus, 
contrary  to  what  one  might  expect,  the  relative  misorientation  of  adjacent 
grains  does  not  appear  to  influence  the  ability  of  a  grain  boundary  to  arrest 
a  small  crack.  The  frequency  ./ i t h  which  grain  boundary  arrests  were  observed 
in  CG  Astroloy  is  in  good  agreement  with  recent  observations  at  Rockwell 
[5.14)  on  Waspaloy,  Figure  5-15.  Interestingly,  the  Rockwell  data  show  that 
grain  boundary  arrests  occur  more  frequently  as  the  grain  size  decreases. 
This  factor  may  have  contributed  to  the  difficulties  encountered  in  obtaining 
small  crack  data  in  FG  Astroloy. 

An  obvious  limitation  of  the  above  analysis  is  the  fact  that  observa¬ 
tions  were  restricted  to  crick  growth  within  the  surface  grains.  If  a 
stronger  correlation  between  temporary  crack  arrests  and  grain  boundaries 
exists  it  would  have  to  involve  interactions  among  adjacent  grains.  In  this 
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Figure  5-13.  Comparison  of  small  crack  data  variability  from  3 
single  data  set  versus  the  entire  data  set  for 
CG  Astroloy  at  25°C. 
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Figure  5-14. 


Crack  length  (top-curve-right  scale)  and  crack  growth  rate 
(lower  curve-left  scale)  for  one  end  of  a  small  crack  in 
Coarse  Grained  Astroloy  at  20°C  (BF-17,  Crack  18).  The 
number  of  grain  boundaries  crossed  are  indicated  across 
the  top.  The  dashed  line  is  the  minimum  value  of  crack 
growth  rate  detectable  by  the  measurement  technique:  how¬ 
ever,  it  does  represent  temporary  crack  arrest.  Circled 
numbers  indicate  crack  arrest  at  mi crostructural  features 
as  follows: 

T  arrest  at  triple  point 

2  arrest  at  GB,  with  subsequent  growth  along  another 
di recti  on 

3\  arrest  at  GB,  then  growth  along  grain  blundaries 
4 '  arrest  at  point  slightly  beyond  a  GB 
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Figure  5-15.  Percentage  of  surface  crack  tips  stopped  by  grain 
boundaries  in  Waspaloy  and  coarse  grained  Astrcloy 
at  room  temperature. 


way  arrest  at  a  grain  boundary  below  the  surface  might  influence  what  was 
observed  at  the  surface,  thereby  accounting  for  the  majority  of  ca^es  when 
arrests  at  the  surface  appeared  not  to  be  associated  with  a  grain  boundary 

As  one  might  expect,  the  data  variability  can  be  reduced  by  averaging 
crack  growth  rates  over  larger  intervals  of  cycles  (or  crack  length)  as  shown 
in  Figure  5-16.  Interestingly,  using  larger  interval  to  analyze  the  data 
produces  variability  which  is  comparable  to  that  obtained  on  large  cracks  in 
this  growth  rate  regime.  This  observation,  combined  with  the  fact  that  large 
cracks  exhibit  this  same  irregular  growth  when  observed  at  high  resolution 
[5.15),  suggests  that  the  observed  difference  in  the  variability  of  crack 
growth  rates  between  large  and  small  cracks  is  due  primarily  to  differences  in 
the  precision  of  measurement.  Thus,  it  would  also  appear  that  the  basic 
mechanism  of  crack  extension  is  the  same  for  both  large  and  small  cracks 
thereby  supporting  the  view  that  the  small  crack  effect  is  due  to  differences 
in  the  mechanical  driving  force  for  long  and  small  cracks. 

The  possible  influence  of  periodic  test  interruptions  and  surface 
replication  on  the  the  growth  of  small  cracks  was  also  examined  by  varying  the 
sampling  interval  between  10,000  and  50,000  cycles.  However,  so  as  not  to 
confound  the  interpretation  of  these  results  by  the  above  mentioned  influence 
of  crack  growth  increment,  all  results  were  compared  using  data  analyzed  with 
a  common  cycle  interval  of  50,000  cycles.  As  shown  in  Figure  5-17,  the 
frequency  of  test  interruption  and  measurement  had  little  influence  on  the 
measured  growth  rates.  The  largest  deviation  in  rates  was  about  a  factor  of 
two  and  occurred  at  the  highest  aK  values  examined.  The  source  of  this  dif¬ 
ference  in  rates  is  not  currently  understood.  However,  if  the  deviation  is 
due  to  the  interruption  of  cycling  (rather  than  to  the  replication  process 
itself)  the  slower  rates  should  be  typical  of  those  which  occur  in  actual 
discs  since  these  are  also  subjected  to  periodic  unloadings  during  service. 

5 . 6  Small  Versus  Long  Crack  Growth  Rates  as  a  Function  of  Material , 

Environment,  Temperature  and  Loading  Variables 

fatigue  crack  growth  rates  for  small  and  long  cracks,  expressed  as  a 
function  of  :.k,  are  compared  in  Figures  5-18  thru  5-24  for  several  applied 
stress  levels,  environments,  test  temperatures,  and  material  conditions.  In 
nearly  ail  cases,  the  classical  small  crack  effect  was  observed  --  that  is. 


Figure  5-16.  Effect  of  crack  growth  increment  on 
computed  growth  rates  for  small 
cracks . 
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Figure  5-17.  Influence  of  frequency  of  measurement  and/or  test  interruption  on  the  growth 
rates  of  small  cracks. 
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Figure  5-18.  Comparison  of  large  and  small  crack  data  in 
CG  Astroloy  at  25°C  for  various  relative 
applied  stress  levels  and  environments. 
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Figure  5-19.  Comparison  of  large  and  small  crack  growth  rate 
data  in  CG  Astroloy  at  600°C  for  various  relative 
applied  stress  levels  and  environments. 
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Figure  5-20.  Comparison  of  large  and  small  crack  growth  rate 
data  in  Waspaloy  at  600°C  for  various  relative 
applied  stress  levels  and  environments. 
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Figure  5-21.  Effect  of  alloy  composition  on  the  crack  growth 
rates  in  both  small  and  large  cracks  at  600°C. 
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Figure  5-22.  Comparison  of  large  (■)  and  "small 
(□,  a  =  215-1000  um)  crack  growth 
rates  in  FG  Astroloy  at  200°C. 
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Figure  5-23.  Compcrison  of  large  and  small 
rates  in  FG  Astroloy  at  600°C 
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Figure  5-24.  Effect  of  temperature  on  the  growth  rates  of 
small  (symbols)  and  large  (lines)  cracks  in 
CG  4. troloy. 


small  cracks  exhibited  faster  rates  than  did  large  cracks,  at  a  given  aK 
value,  and  grew  at  aK  values  below  the  long  crack  threshold  stress  intensity, 

Kth* 

5.6.1  Effect  of  Applied  Stress  Level 

Figure  5-18  and  Figures  5-19  and  5-20  show  the  influence  of 
applied  stress  on  the  kinetics  of  small  crack  growth  at  temperatures  of  25°C 
and  60CPC,  respectively.  The  applied  stress  level  in  these  figures  is 
expressed  in  terms  of  the  relative  stress,  S,  defined  as  the  ratio  of  the 
applied  stress  to  material  yield  strength  (Smax/Sys).  The  CG  Astroloy  data  at 
25°C  in  Figure  5-18  shows  that  S  values  of  0.8  arid  1.1  give  similar  results, 
while  growth  rates  at  S  =  1.6  are  significantly  faster,  at  a  given  value 
of  aK.  On  the  other  hand,  as  shown  in  Figures  5-19  and  5-20,  growth  rates  in 
both  CG  Astroloy  and  Waspaloy  at  600°C  are  relatively  insensitive  to  S  for 
values  ranging  from  1.5  to  1.8.  Elevated  temperature  data  are  not  available 
at  lower  S  values  due  to  the  difficulty  in  initiating  and  growing  cracks  at 
600"C,  while  data  at  higher  S  values  were  not  obtained  since  S  =  1.8  was  the 
limit  for  global  elastic  response  of  the  bend  specimens.  If  data  were 
available  over  a  wider  range  of  S  values  at  600°C,  results  similar  to  those 
obtained  at  room  temperature  would  be  expected. 

The  above  results  suggests  that  the  influence  of  S  on  small  crack 
growth  kinetics  is  small  at  low  S  values,  larger  at  intermediate  S  values,  and 
saturates  at  higher  S  values.  The  latter  behavior  is  consistent  with  the 
alteration  of  the  local  stress  state  due  to  localized  yielding  at  the  surface 
of  the  bend  spec  men  as  described  in  Section  5.2.  This  localized  yielding 
causes  an  increase  in  S  to  be  counteracted  by  a  decrease  in  the  local  mean 
stress  (Figure  5-4)  and  concomitant  decrease  in  the  local  stress  range  which 
is  effective  in  opening  the  crack. 

It  is  significant  to  point  out  that  the  dependence  of  the  small 
crack  growth  kinetics  on  applied  stress  level  suggests  that  aK  provides  an 
inadequate  characterization  of  the  mechanical  driving  force  for  small 
cracks.  This  occurs  in  spite  of  the  fact  that  the  overall  specimen  response 
is  elastic  and  is  presumably  due  to  the  fact  that  the  small  scale  yielding 
requirement  of  linear-elastic  fracture  mechanics  also  requires  that  the  extent 


of  plasticity  be  small  compared  to  the  crack  size  --  a  condition  whicn  is 
violated  in  most  small  crack  experiments,  as  discussed  in  greater  detail  in 
Section  7. 

5.6.2  Effect  of  Environaent 

figures  5-18  and  5-19  alsc  provide  a  comparison  of  small  crack 
data  obtained  in  two  different  environments  --  vacuum  versus  air.  The  air 
data  were  obtained  using  the  three-point  bend  specimen  described  in  Section 
5.1,  while  the  vacuum  data  were  obtained  using  the  single  edge  notch  specimen 
described  in  Section  7.1.  No  measurable  difference  in  rates  were  observed. 
This  fact  is  believed  to  be  due  the  relatively  high  cyclic  frequency  of  10  Hz 
which  was  employed,  thereby  precluding  any  significant  interaction  with  the 
environment,  even  at  600°C.  This  interpretation  is  consistent  with  the  fact 
that  no  significant  difference  was  observed  between  the  growth  ^ate  of  long 
cracks  in  air  and  vacuum  tests,  below  1X10'4  mm/cycle  (Section  4).  At  higher 
growth  rates,  the  long  crack  data  obtained  in  air  are  about  a  factor  of  two 
faster  than  the  vacuum  data. 

5.6.3  Effect  of  Material  Variables 

Figure  5-21  compares  small  and  large  crack  data  from  CG  Astroloy 
and  Waspaloy  at  600"C.  The  long  crack  data  for  both  alloys  fall  within  a  re¬ 
latively  narrow  scatter  band.  However,  if  one  were  to  carefully  examine  these 
long  crack  results,  a  small  difference  between  the  two  alloys  actually  occurs 
which  is  consistent  with  the  influence  of  grain  size  which  was  discussed 
previously  in  Section  4.  In  contrast,  no  measurable  difference  was  observed 
between  the  small  crack  results  in  the  two  alloys  at  600°C.  This  observation 
is  consistent  with  the  room  temperature  results  of  Brown  et  al  (5.13]  on  these 
same  alloys  which  showed  small  crack  results  to  be  relatively  insensitive  to 
grain  size,  as  well  as  y'  morphology. 

Results  on  FG  Astroloy  at  200’C,  Figure  5-22,  show  good  agreement 
between  the  small  surface  crack  data  and  long  crack  data.  The  lack  of  the 
crack  size  effect  in  this  case  is  likely  due  to  the  fact  that  the  surface 
crack  length  (2a  =  430  to  2000  urn)  was  always  significantly  larger  than  the 
12  t.m  grain  size  of  this  material.  This  conclusion  is  in  agreement  with  the 
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room  temperature  data  of  Brown  et  al  (5.13|  which  showed  that  a  small  crack 
effect  in  FG  Astroloy  of  this  same  grain  size  was  limited  to  cracks  which  were 
smaller  than  about  300  um. 

As  shown  in  Figure  5-23,  results  on  FG  Astroloy  at  600°C  also  did 
not  exhibit  the  classical  small  crack  effect.  Only  tne  smallest  crack  sizes, 
corresponding  to  the  lowest  AK  values,  exhibited  rates  which  were  faster  than 
the  long  crack  data.  Curiously,  as  the  surface  crack  enlarged  the  rates  de¬ 
creased  below  those  of  the  long-thru  cracks.  Instead  of  merging  with  the  long 
crack  data.  This  difference  In  rates  Is  not  likely  to  be  due  to  the  differ¬ 
ence  in  environment  in  view  of  the  similarity  in  rates  observed  In  air  versus 
vacuum  in  CG  Astroloy  at  600°C,  Figure  5-19.  The  overall  trend  in  the  above 
results  is  remarkably  similar  to  the  room  temperature  results  of  Lankford  et 
al  (5.16|  on  PM  alloy  IN-100  having  a  grain  size  of  3  to  7  um.  However,  the 
current  results  appear  to  be  inconsistent  with  the  room  temperature  data  of 
Brown  et  al  [5.13]  which  showed  FG  Astroloy  to  exhibit  rates  which  were  com¬ 
parable  to  those  in  CG  Astroloy,  both  of  which  were  significantly  faster  than 
the  long  crack  rates.  Obviously,  additional  experiments  are  needed  to  resolve 
this  dlcotomy  in  the  small  crack  data  in  fine  grained  material.  Unfortun¬ 
ately,  these  data  are  extremely  difficult  to  acquire  due  to  the  small  crack 
sizes  involved  and  the  difficulty  in  initiating  cracks  in  fine  grained  mate¬ 
rials. 


5.6.4  Effect  of  Temperature 


Figure  5-24  summarizes  the  effect  of  temperature  on  the  growth  of 
small  and  large  cracks  in  CG  Astroloy.  At  room  temperature,  only  small  crack 
results  at  S  =  1.6  are  included  so  as  not  to  confound  the  comparison  with  the 
previously  discussed  effect  of  applied  stress  level.  As  one  might  expect,  the 
large  crack  growth  rates  increase  with  increasing  temperature  over  most  of 
the  aK  regime.  The  major  exception  to  this  trend  occurs  in  the  near-threshold 
reoime  where  the  600°C  data  cross-over  the  200°C  data  and  consequently  ex¬ 
hibits  a  higher  A«th  value.  This  cross-over  of  the  near-threshold  growth 
rates  with  Increasing  temperature  is  consistent  with  recent  observations  in 
other  Ni-base  superalloys  (5.17,  5.18].  At  low  aK  values,  the  small  crack 
growth  rates  do  not  exhibit  a  clear  trend  with  respect  to  temperature.  How¬ 
ever,  as  the  small  cracks  enlarge,  they  merge  with  the  long  crack  data  and 


The  complexity  of  the  temperature  dependence  for  both  long  and 
small  cracks  Is  believed  to  be  related  to  changes  In  the  character  of  the  de¬ 
formation  with  increasing  temperature.  Specifically,  as  Indicated  in  Section 
6.3,  the  change  in  fracture  surfaces  for  both  small  and  long  cracks  from  those 
having  highly  angular  facets  to  those  which  are  smoother  and  flatter  suggests 
that  slip  becomes  more  homogeneous  with  increasing  temperature.  These  results 
appear  to  be  consistent  with  observations  made  during  low  cycle  fatigue  of  Ni- 
base  alloys  which  shown  the  deformation  to  become  more  homogeneous  at  elevated 
temperature,  presumably  due  to  the  thermal  activation  of  cross-slip  and  climb, 
particularly  for  materials  with  coarse  grains  and  fine  precipitates  [5.19  - 
5.22].  Thus,  although  slip  is  still  predominantly  planar,  it  ran  apparently 
occur  on  intersecting  <111>  planes  to  produce  a  wavy  appearance.  These  pro¬ 
cesses  are  also  consistent  with  the  observed  time-dependent  deformation  of 
small  cracks  which  arrested  at  6QQ°C  after  growing  through  the  first  few 
grains  (Section  5.4). 

Although  the  above  explanation  of  the  complex  temperature 
dependence  of  crack  growth  in  these  alloys  is  plausible,  certain  details  of 
the  deformation  mechanism  attending  crack  growth  remain  unclear.  For  example, 
as  shown  in  Section  7,  the  deformation  fields  at  small  cracks  are  relatively 
homogeneous  at  all  temperatures,  thus  the  relation  between  the  crystallo¬ 
graphic  nature  of  the  crack  path  and  the  near-tip  deformation  field  is  not 
fundamentally  understood. 

In  spite  of  this  incomplete  mechanistic  understanding,  the 
current  results  clearly  demonstrate  that  the  small  crack  effect  is  not  just  a 
low  temperature  phenomenon,  thus  it  can  bp  of  potential  practical  importance 
to  life  prediction  in  aeroengine  components.  This  issue  is  examined  in 
Section  9. 

5.7  References 

5.1  The  Application  of  3-Dimensional  Finite  Element  Methods  to  Fracture 
Mechanics  and  Fatigue  Life  Prediction,  A.  C.  Pickard,  Engineering 
Materials  Advisory  Services,  Ltd.,  West  Midlands,  U.K. 

5.2  U.  S.  Lindholm,  K.  S.  Chan,  S.  R.  Bodner,  R.  M.  Weber,  K.  P.  Walker  and 
B.  N.  Cassenti,  "Constitutive  Modeling  for  Isotropic  Materials  (HOST)  - 
Second  Annual  Status  Report,"  NASA  CR-174980,  July  1985. 


5.3  JoDean  Morrow  and  G.  M.  Sinclair,  “Cycle-Dependent  Stress  Relaxation," 
ASTM  STP  237,  1958,  pp.  83-109. 

5.4  K.  S.  Chan  and  T.  A.  Cruse,  "Stress  Intensity  Factors  for  Anisotropic 

Compact-Tension  Specimens  with  Inclined  Cracks,"  Engineering  Fracture 
Mechanics,  Vol.  23,  No.  5,  1986,  pp.  863-874. 

5.5  ASTM  E  647:  "Standard  Test-Method  for  Constant-Load-Amplitude  Fatigue 
Crack  Growth  Rate  Above  10"°  m/cycle,"  Annual  Book  of  ASTM  Standards, 
Vol.  03.01,  1986,  pp.  714-736. 

5.6  K.  S.  Chan,  J.  E.  Hack,  and  G.  R.  Leverant,  "Fatigue  Crack  Growth  in 

MAR-M200  Single  Crystals,"  Met.  Trans..  Vol.  18A,  April  1987,  pp.  581- 
591. 

5.7  K.  S.  Chan  md  G.  R.  Leverant,  "Elevated-Temperature  Fatigue  Crack 

Growth  Behavior  of  MAR-M200  Single  Crystals,"  Met.  Trans.,  Vol.  18A, 
April  1987,  pp.  593-602. 

5.8  A.  C.  Pickard,  M.  A.  Hicks  and  R.  H.  Jeal,  "Fatigue  '87  Overview 

Lecture:  Applications  of  Fatigue  Analyses:  Aircraft  Engines,"  in 

Fatigue  '87,  Eds.  R.  0.  Ritchie  and  E.  A.  Starke,  Jr.,  Engineering 
Materials  Advisory  Service,  U.K.,  1987. 

5.9  Rolls-Royce,  unpublished  results. 

5.10  A.  T.  Stewart,  "The  Influence  of  Environment  and  Stress  Ratio  on  Fatigue 
Crack  Growth  at  Near  Threshold  Stress  Intensities  in  Low  Alloy  Steels," 
Eng.  Fract.  Mech.,  Vol.  13,  1980,  pp.  463-478. 

5.11  R.  0.  Ritchie,  S.  Suresh  and  C.  M.  Moss,  "Near-Threshold  Fatigue  Crack 
Growth  in  2-*s  Cr-Mo  Pressure  Vessel  Steel  in  Air  and  Hydrogen,"  J.  of 
Engineering  and  Materials  Technology,  Trans.  ASME,  Series  H.,  Vol.  102, 
1980,  pp.  293-299. 

5.12  W.  G.  Clark  and  S.  J.  Hudak,  Jr.,  "Variability  in  Fatigue  Crack  Growth 
Rate  Testing,"  J.  of  Testing  and  Evaluation.  JTEVA,  Vol.  3,  No.  6,  1975, 
pp.  454-476. 

5.13  C.  W.  Grown,  J.  E.  King  and  M.  A.  Hicks,  "Effects  of  Microstructure  on 
Long  and  Short  Crack  Growth  in  Nickel  Base  Superalloys,"  Metal  Science, 
Vol.  18,  July  1984,  pp.  374-380. 

5.14  W.  J.  Pardee  et  al,  "Quantitative  Nondestructive  Evaluation  (NDE)  for 
Retirement-For-Cause,"  Rockwell  International  Science  Center,  Report 
No.  5C5283. 16FR,  Final  Report  on  DARPA  Contract  No.  MOA  903-80-C-0641, 
1984. 

5.15  0.  L.  Oavidson  and  J.  Lankford,  "The  Breakdown  of  Crack-Tip 
Microstructure  Ouring  Fatigue  Crack  Extension  in  Aluminum  Alloys,"  in 
High  Strength  Powder  Metallurgy  Aluminum  Alloys,  Eds.  G.  J.  Hi  1  deman  and 
M.  J.  kccr fik ,  i MS/A IME ,  1986. 


5.16  J.  Lankford,  T.  S.  Cook  and  G.  P.  Sheldon,  "Fatigue  Microcrack  Growth  in 
a  Nickel-Base  Superalloy,"  Int.  J.  of  Fracture.  Vol.  17,  April  1981,  pp. 
143-155. 

5.17  J.  L.  Yuen,  P.  Rox  and  W.  0.  Nix,  "Effect  of  Oxidation  Kinetics  on  the 
Near  Threshold  Fatigue  Crack  Growth  Behavior  of  a  Nickel  Base 
Superalloy,"  Met.  Trans.,  Vol.  15A,  September  1984,  pp.  1769-1775. 

5.18  R.  H.  Van  Stone  and  D.  D.  Krueger,  "Closure  and  Near-Threshold  Crack 

Growth  in  Nickel-Base  Superalloys,"  Paper  Presented  at  the  19th  National 
Symposium  on  Fracture  Mechanics,  San  Antonio,  TX,  June  1986. 

5.19  M.  Clavel,  C.  Levaillant  and  A.  Pineau,  "Influence  of  Micromechanisms  of 
Cyclic  Deformation  and  Elevated  Temperature  on  Fatigue  Behavior,"  in 
Creep-Fatigue  Environment  Interactions,  Eds.  R.  M.  Pelloux  and  N.  S. 
Stoloff ,  TMS/AIME,  Warrendale.  PA,  1980,  pp.  24-45. 

5.20  R.  V.  Miner,  J.  Gayda  and  R.  0.  Maler,  "Fatigue  and  Creep-Fatigue 

Deformation  of  Several  Nickel-Base  Superalloys  at  650°C,"  Met.  Trans., 
Vol.  13A,  October  1982,  pp.  1755-1765. 

5.21  B.  A.  Lerch,  N.  Jayaraman  and  S.  D.  Antolovich,  Mater.  Sci.  Eng., 

Vol.  66,  1984,  pp.  151-166. 

5.22  J.  E.  King,  "Fatigue  Crack  Propagation  in  Nickel-Base  Superalloys  -- 

Effects  of  Microstructure,  Load  Ratio  and  Temperature,"  Mat.  Sci.  and 
Tech. ,  Vol.  3,  September  1987,  pp.  750-764. 


6.0  MECHANISMS  OF  SMALL  CRACK  INITIATION  AND  GROWTH 


This  section  presents  information  on  the  basic  mechanisms  of  fatigue 
crack  initiation  and  growth  in  the  Nl-base  alloys  examined  in  this  study. 
Results  are  summarized  in  three  parts:  First  the  role  of  pre-existing  pores 
and  persistent  slip  bands  on  the  initiation  and  early  growth  of  small  cracks 
is  considered.  Next,  the  orientation  relationship  between  these  cracks  and 
the  underlying  crystallography  within  individual  grains  is  determined.  These 
results  not  only  serve  to  elucidate  the  initiation  process,  but  also  provide 
information  essential  to  understanding  subsequent  growth  of  small  cracks  since 
the  displacement  modes  experienced  by  these  cracks  are  controlled  by  their 
orientation  within  the  crystallites.  Finally,  the  morphology  of  the  fracture 
surfaces  are  examined  as  a  function  of  the  primary  material,  loading  and  en¬ 
vironmental  variables  studied  In  this  program.  This  information  proved  to  be 
particularly  useful  in  understanding  the  similarities  and  differences  in  the 
growth  mechanisms  of  small  and  long  cracks,  as  well  as  the  complex  temperature 
dependence  of  their  growth  rates. 

6.1  The  Roles  of  Pores  and  Persistent  Slip  Bands  In  Crack  Initiation  and 

Early  Growth 

Many  of  our  basic  concepts  of  fatigue  crack  initiation  and  early  growth 
of  microcracks  have  come  from  mechanistic  studies  of  model  system,  such  as 
single  crystals  of  copper  or  iron-silicon  —  for  example,  see  Refs.  [6.1  thru 
6.4).  These  basic  studies  have  emphasized  the  role  of  persistent  slip  bands 
and  the  subsequent  formation  of  surface  intrusions  and  extrusions  in 
nucleating  cracks.  The  primary  objective  of  this  portion  of  the  current  work 
was  to  assess  whether  or  not  this  crack  initiation  mechanism  is  applicable  to 
Ni-based  superalloys.  Although  this  effort  concentrated  on  CG  Astro loy  at 
room  temperature,  general  features  of  the  observations  are  believed  to  be 
applicable  to  FG  Astroloy  and  Waspaloy. 

The  exoerimental  procedures  were  the  same  as  those  used  to  jtrdy  tir 
lora:  ueformat ion  at  microcracks  (Section  7.1),  with  one  important  modifica¬ 
tion.  Namely,  replication  was  not  used  during  the  current  experiments  so  as 
not  to  damage  the  extrusions,  since  these  features  are  known  to  be  fragile 


loading  to  about  80  percent  of  its  expected  fatigue  life,  then  observed 
directly  within  the  SEM  at  various  tilt  angles. 

In  all  cases,  fatigue  cracks  were  found  to  initiate  within  persistent 
slip  band,  occasionally  these  bands  were  formed  at  surface  pores.  Observation 
from  eight  cracks  may  be  classified  into  5  categories,  each  illustrated  by  a 
figure,  as  follows: 

1.  Pore  initiated  crack  showing  no  extrusions  or  debris.  A  low  magnifica¬ 
tion  view  of  two  fatigue  cracks  that  initiated  at  adjacent  pores  and 
linked  together  through  subsequent  growth  is  shown  in  Figure  6-l(a), 
while  Figure  6- 1 (b)  shows  a  detail  of  a  crack  traversing  a  pore.  The 
cubic  shape  of  v'-precipitates  at  the  bottom  of  the  pore  suggests  that 
the  crack  plane  is  of  the  <11 1>  type. 

2.  Pore  initiated  crack  showing  debris  or  tongue-like  extrusions  and  Mode 
III  displacements.  A  lower  magnification  view  is  shown  in  Figure  6-2(a) 
at  a  lower  tilt  angle,  while  Figure  6-2(b)  is  a  higher  magnification 
view  at  a  higher  tilt  angle.  The  figure  shows  debris,  or  extrusions, 
coming  from  the  corner  of  the  pore  and  from  a  small  crack  nearby. 
Evidently,  the  small  crack  showing  the  more  profuse  and  continuous 
extrusions  initiated,  but  did  not  grow.  In  contrast,  the  larger  crack 
evidences  debris,  less  continuous  in  nature  than  at  the  small  crack,  at 
isolated  locations  along  its  length. 

3.  Slip  band  initiated  cracks  showing  Mode  III  opening  but  no  debris  or 

extrusions.  A  low  magnification  view  of  the  whole  crack  is  shown  in 
Figure  6-3 (a ) .  The  direction  of  the  loading  axis  is  unknown  for  this 

crack,  observations  from  similar  cracks  showed  that  they  were  usually  at 
an  angle  of  40-55°  to  the  loading  axis  (see  Section  6.2).  The  white 
object  projecting  above  the  surface  on  the  figure  may  also  be  seen  on 
the  left  side  of  Figure  6-3(b)  which  is  a  view  looking  down  the  crack 
from  a  grazing  angle  to  the  surface.  Note  that  this  crack  is  open, 
similar  to  the  cracks  shown  in  Figures  6-1  and  6-2,  and  that  some 
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FIGURE  6-1.  Fatigue  crack  initiation  at  surface  pores  in  CG  Astroloy 
showing  (a)  an  overview  of  the  crack  with  some  debris 
along  the  lower  part  of  the  crack,  and  (b)  a  higher  mag¬ 
nification  view  of  one  of  the  initiation  sites. 
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FlbUKL  6-2 .  Two  views  of  fatigue  crack  initiation  at  a  surface  pore 
showing  debris  at  t he  corner  of  the  pore  and  along  the 
crack  further  away.  The  tilt  of  (a)  is  lower  than  (b), 
which  shows  better  the  somewhat  "tongue- 1 i ko"  nature  of 
the  debris.  Note  the  large  amount  of  debris  extruded 
from  the  very  small  crack  adjacent  to  the  main  crack. 
Ibis  crack  apparently  initiated  but  did  not  grow. 
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Slip  band  initiated  crack  showing  Mode  III,  and  having  associated 
debris.  An  overall  view  of  the  crack  is  shown  in  Figure  6-4(a).  The 
initiation  site  was  along  one  of  the  slip  bands  evident  in  the  one  grain 
showing  extensive  slip.  A  higher  magnification  view  of  this  region  is 
shown  in  Figure  6-4(b),  which  shows  the  crack  crossing  two  boundaries. 
At  one  boundary,  the  crack  continued  across  the  adjacent  grain  in  the 
same  direction,  but  finally  changed  direction  drastically,  continuing 
along  the  grain  boundary.  The  other  end  of  the  crack  also  changed 
direction  radically,  but  debris  accompanied  this  part  of  the  crack  and 
nearly  obscured  the  crack  path  over  a  short  region.  This  debris  does 
not  have  the  nature  of  organized  extrusions,  which  have  also  been  termed 
"tongues"  such  as  shown  accompanying  the  small  crack  in  Figure  6-2. 
Presumably  the  debris  generated  by  this  crack  resulted  from  the  change 
in  direction  of  the  crack,  which  caused  crack  face  rubbing.  The 
Mode  III  opening  of  the  crack  is  likely  to  have  ejected  this  debris  from 
within  the  crack. 

Figure  6-5  shows  the  relationship  between  the  slip  lines  and  the  y'- 
precipitates.  These  precipitates  are  cube-shaped,  with  sides  of  the 
cube  lying  along  {100};  thus  the  square  shapes  of  the  precipitates 
indicate  that  the  surface  of  observation  is  on  or  near  the  [100}.  Slip 
lines  cut  across  these  precipitates  on  the  diagonal,  indicating  slip  on 
{111}  which  are  inclined  to  the  surface  at  an  angle  of  approximately 
35  degrees.  Slip  could  be  along  either  of  the  two  <110>  Burgers  vectors 
lying  on  the  {111}  and  inclined  to  the  surface  at  45  degrees,  or  the 
third  <110>  Burgers  vector  in  the  plane  of  the  surface.  The  topography 
of  the  surface  shown  Figure  6-5,  which  was  also  confirmed  by  stereo¬ 
photographs,  indicates  that  at  least  part  of  the  slip  was  out  of  the 
plane  of  the  surface.  Therefore,  some  Mode  III  displacements  very 
likely  accompanied  the  crack  of  Figure  6-4. 

5.  Slip  initiated  cracking  accompanied  by  tongue-like  extrusions.  The  only 
crack  clearly  having  the  same  type  cf  extrusions  as  those  reported  in 
single  crystals  is  shown  in  Figure  6-6.  Photographs  (a)  thru  (c)  a^e 
shown  at  progressively  larger  amounts  of  tilt,  each  of  which  more 


URf  6-5.  Relationship  between  the  ^'-precipitate  structure  and  si  ip 
lines  showing  that  slip  is  of  <111>  type  and  is  likely  tc 
result  in  both  in-plane  and  out-of-plane  displacements. 
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Three  views  of  d  s' ip  band  initiated  tatigue  crack  show'ng 
extrusions  having  tongue-like  shapes.  These  micrographs 
were  made  at  a  decreasing  angle  of  incidence,  with  (c) 
being  at  a  grazing  angle  to  the  surface.  This  crack  appears 
to  have  initiated  within  a  single  grain  at  the  intersection 
of  two  slip  systems  and  shows  little  evidence  for  Mode  III 
d i sp 1 acements . 


Figures  6-4  and  6-5.  This  is  believed  to  be  partly  due  to  the  tilt  used 
in  the  photographs  and  partly  due  to  *he  fact  that  the  surface  plane  of 
the  specimen  does  not  appear  to  be  a  low  index  plane,  as  determined  from 
the  shape  of  the  precipitates.  The  grain  in  which  this  crack  initiated 
extends  from  the  line  of  blocky  y'-precipitates  in  the  left-center  of 
Figure  6-6(a)  to  the  line  of  blocky  y'-precipitates  very  near  the  top  of 
the  photograph.  Two  different  slip  systems  are  evident  within  this 
grain,  one  near  the  lower  left  boundary,  which  has  resulted  in  the  slip 
lines  pointing  vertically  in  the  photograph,  and  the  other  running 
approximately  parallel  to  the  crack  plane  at  the  upper  right  of  the 
figure.  The  crack  appears  to  have  formed  approximately  at  the  inter¬ 
section  of  these  two  slip  systems.  Although  one  might  expect  this  in¬ 
tersection  to  be  a  twin  boundary,  the  shapes  of  the  precipitate 
particles  do  not  support  this  expectation.  If  this  were  the  case,  it 
would  parallel  the  observation  of  Neumann  16.4] ,  who  found  cracks  in 
single  crystals  of  copper  growing  in  twin  boundaries  which  were  accom¬ 
panied  by  tongue-like  extrusions. 

In  conclusion,  it  has  been  shown  that  fatigue  cracks  in  CG  Astroloy 
start  due  to  slip  initiated,  either  in  favorably  oriented  grains,  or  at  sur¬ 
face  pores.  For  cracks  initiated  within  grains,  slip  on  [Ill]  crystallo¬ 
graphic  planes  probably  occurs  along  several  <110>  directions  simultaneously, 
causing  both  in-plane  and  out-of-plane  displacements.  Under  some  circum¬ 
stances,  this  can  cause  the  appearance  of  extrusions  of  material  having  a 
tongue-like  morphology,  but  such  extrusions  are  by  no  means  necessary  for 
crack  formation.  Debris  may  also  be  extruded  from  small  cracks  at  or  near 
initiation  sites,  but  this  debris  is  a  consequence  of  crack  formation  rather 
than  the  cause  of  it.  Debris  apparently  occurs  due  to  crack  face  rubbing, 
probably  accompanied  by  Mode  III  displacements.  Tongues  may  be  considered 
just  another  form  of  debris,  or  they  may  be  caused  by  some  more  organized 
mechanism,  such  as  has  been  detailed  by  Neumann  | 6 . 4 ] . 

6.?  The  Crystallography  of  Microcrack  Initiation  and  Growth 

The  primary  objective  of  this  portion  of  the  work  was  to  examine  the 
orientation  relationship  between  microcracks  and  the  underlying  crystal lo- 


channeling  technique  16.5]  was  used  to  determine  this  information.  The  test 
procedure  was  similar  to  that  given  in  Section  7.1.  However,  just  prior  to 
commencing  the  electron  channeling  work,  the  specimen  was  given  an  additional 
brief  electropolish  to  remove  any  oxide  which  might  have  formed  during  the 
experiment;  the  material  removed  by  this  process  was  insufficient  to  destroy 
the  cracked  grains.  Since  the  average  grain  size  of  this  material  was  30  ym 
and  the  spot  size  used  for  making  selected  area  electron  channeling  patterns 
(SACP)  was  15  ym,  information  was  obtained  from  within  individual  grains.  The 
fact  that  SACP  could  be  measured  meant  that  the  dislocation  density  in  the 
cracked  grains  was  below  about  10^/sq.  cm.  Previous  estimates  have  shown  that 
this  dislocation  density  corresponds  to  a  cumulative  strain  on  the  surface  of 
no  greater  than  about  10  percent  [6.6],  Naturally,  the  localized  strain 
within  the  slip  band  is  expected  to  be  much  higher. 

It  was  possible  to  obtain  SACP  of  sufficient  clarity  to  determine  the 
orientation  of  10  grains  containing  cracks.  A  summary  of  the  orientations  of 
the  grain  normals  is  shown  in  the  stereographic  triangle  of  Figure  6-7.  This 
figure  also  shows  a  face  centered  cubic  channeling  map  (made  from  a  copper 
single  crystal)  which  was  used  for  orienting  the  SACP  from  Astroloy.  Orienta¬ 
tions  of  surface  normals  of  grains  forming  cracks  were  grouped  near  <100>, 
<001>  and  <113>. 

By  knowing  the  angular  rotation  between  the  SEM  micrograph  of  the 
cracked  grain,  the  orientation  of  the  SACP,  and  the  geometric  relationship  of 
the  crack  to  the  loading  axis,  it  was  possible  to  derive  relationships  between 
the  crack  plane,  loading  axis,  grain  surface  normal  and  slip  direction. 
Definitions  of  these  angular  relations  are  shown  in  Figure  6-8.  The  associated 
angles  derived  for  each  of  the  oriented  grains  are  given  in  Table  6-1,  with 
examples  shown  in  Figures  6-9  and  6-10.  The  grain  shown  in  Figure  6-9  was 
chosen  as  an  example  because  the  crack  initiated  at  the  pore,  seen  as  a  black 
hole  in  the  backscattered  electron  image  of  Figure  6-9(a).  As  shown  in 
Figure  6-9(c),  the  orientation  of  that  grain  is  near  the  <110>.  The  remainder 
of  the  angular  relations  are  shown  in  Figure  6-9(d).  The  other  mode  of  crack 
initiation,  that  caused  by  slip  band  formation  within  a  grain,  is  illustrated 
in  Figure  6-10.  Note  that  in  both  illustrations,  one  of  the  slip  directions 
was  always  within  40  degrees  of  the  surface  normal,  indicating  that  Mode  III 


FIGURE  6-8.  Definition  of  geometric  relations  between 

L  -  loading  direction 
l  -  crack  direction 
n  -  surface  normal 
p  -  crack  plane  normal 
t  -  crack  plane  tangent 
b  -  slip  direction 

plane  of  <ji  is  defined  by  n,  t, 
plane  of  $  is  defined  by  t,  b, 
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TABLE  6-1 


ORIENTATION  RELATIONSHIPS  FOR  CRACKS  IN  CG  ASTROLOY 
(See  Figure  6-8  for  Definition  of  Orientation  Parameters) 


Crack 

Number 

Length 

No.  Grains 
Crossed 

end 

8 

L-L 

y 

n-t 

♦ 

y> 

4 

240 

1 

1 

55 

33 

31 

3 

2 

68 

33 

31 

5 

102 

2 

1 

64 

0 

7 

7 

167 

2.3 

1 

48 

6 

31 

3 

78 

1.2 

1 

64 

17 

31 

6 

59 

1 

1 

53 

0 

13 

8 

102 

1 

1 

47 

7 

7 

2 

2 

85 

0 

5 

10 

51 

2.3 

1 

43 

9 
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FIGURE  6-9.  Determination  of  crack  growth  plane  and  slip  direction  by 
selected  area  electron  channeling,  (a)  backscattered 
electron  image  of  crack,  showing  grain  pattern  and  location 
of  the  channeling  pattern  shown  in  (b).  Note  there  is  a 
rotation  between  the  pattern  and  image  which  has  not  been 
shown  in  this  illustration,  (c)  shows  the  location  (large 
dot)  of  the  surface  normal  to  the  grain  within  a  stereo¬ 
graphic  triangle,  and  (d)  is  a  sketch  of  the  angular  rela¬ 
tions  derived  for  this  crack.  For  further  definition  of 
the  vectors  shown,  consult  Figure  6-8. 
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FIGURE  6-10.  Determination  of  crack  growth  plane  and  slip  direction  by 
selected  area  electron  channeling,  (a)  backscattered 
electron  image  of  crack,  showing  grain  pattern  and  location 
of  the  channeling  pattern  shown  in  (b).  Note  there  is  a 
rotation  between  the  pattern  and  image  which  has  not  been 
shown  in  this  illustration,  (c)  shows  the  location  (large 
dot)  of  the  surface  normal  to  the  grain  within  a  stereo¬ 
graphic  triangle,  and  (d)  is  a  sketch  of  the  angular  re¬ 
lations  derived  for  this  crack.  For  further  definition  of 
the  vectors  shown,  consult  Figure  6-8. 


Clustering  of  the  surface  normals  of  grains  containing  fatigue  cracKs 
occurs  in  several  regions  of  the  standard  stereographic  triangle.  This  sug¬ 
gests  that  there  are  preferred  grain  orientations  for  crack  initiation  in 
Astroloy.  Fatigue  cracks  in  Astroloy  were  observed  to  initiate  at  slip  bands 
in  grains  which  do  not  contain  pores,  as  well  as  at  slip  bands  emanating  from 
pores;  therefore,  the  role  of  preferred  grain  orientation  in  initiation  of 
fatigue  cracks  in  this  alloy  is  somewhat  unclear.  However,  this  can  be 
examined  by  considering  the  Taylor  factors  of  individual  grains  within  which 
fatigue  cracks  initiated. 

The  appropriate  value  of  the  Taylor  factor,  M,  for  a  given  grain 
orientation  can  be  computed  on  the  basis  of  the  Bishop-Hill  approach  [6.6, 
6 . 7 1 .  In  particular,  the  Taylor  factor  calculations  have  been  performed  using 
[6.7,  6.8) : 

M  =  ^ —  [ -B  +  M^22  +  2  +  2G  de^  +  2H  de^]  (6-1) 

in  which  deyy  is  the  incremental  plastic  strain  In  the  principal  loading  (y) 
axis;  A,B,F,G,  and  H  represent  the  stress  states  capable  of  activating  the 
same  shear  stress  on  five  or  more  slip  systems.  Although  Astroloy  is  com¬ 
prised  of  y 1  (l_i2  structure)  In  a  y  (fee)  matrix,  the  predominate  slip 
systems  at  room  temperature  are  (111}  <110>,  i.e.,  they  are  similar  to  those 
in  cubic  metals  such  as  aluminum  and  copper.  The  number  of  these  stress 
states  is  fifty-six  for  cubic  metals  which  slip  on  {111}  <110>  slip  systems. 
The  values  of  A,B,F,G,  and  H  are  given  in  Ref.  [6.7].  The  plastic  strain  com¬ 
ponents  (de^,  de^t  etc...)  in  Eq.  (6-1)  are  obtained  by  transforming  the 
imposed  plastic  strains  (dexx,  dc^,  and  dezz)  into  the  cube  axes  (1,  2,  and  3 
axes)  of  the  crystals.  In  this  case,  the  loading  direction  is  the  y-axis  and 
the  z-axis  is  parallel  to  the  surface  normal  of  an  individual  grain.  Accord¬ 
ing  to  Bishop-Hill  [6.6,  6.7|,  the  appropriate  stress  state  for  Eq.  (6-1)  is 
the  one  which  maximizes  the  external  work.  The  value  of  M  for  a  given  grain 
orientation  is  therefore  obtained  by  maximizing  the  right-hand  side  of  Eq.  (6- 
1). 

Table  6-2  shows  the  results  of  the  Taylor  factor  calculation  for 
selected  grain  orientations  for  fee  crystals  deforming  by  (111)  <110>  slip 

unHor  f  Ka  nlano.  c>  r.a  in  rnnH  i  f-  i  nn  ( r\<r  =  -Hr  Hr  =  H  \  Tf  chnulH  ho  nntoH 


TABLE  6-2 


TAYLOR  FACTOR  CALCULATION  FOR  SELECTED  GRAIN  ORIENTATIONS  UNDER  PLAIN 
STRAIN  DEFORMATION  (de„=0).  The  average  value  of  M  for  all  possible 
orientations  Is  M  =  3.53. 


Orientation 

X 

Y 

z 

M 

-EAmI 

°1 

6°  off  (0211 

6°  off  (201) 

(Tl2 1 

2.0 

0.565 

°2 

4“  off  (133! 

2°  off  (71T) 

(Oil) 

2.258 

0.640 

°3 

4°  off  (T331 

2°  off  [7X11 

(Oil) 

2.258 

0.640 

°4 

(0011 

(010) 

[0011 

2.449 

0.693 

°5 

l  llO] 

(001) 

[1101 

2.449 

0.693 

°6 

[  111  1 

[1101 

[1121 

4.082 

1.155 

°7 

[T12] 

(111  I 

[TTll 

3.674 

1.040 

°8 

illoi 

(1101 

[0011 

4.899 

1.386 

°9 

(1121 

(1101 

[1TT1 

3.266 

0.924 

that  the  Taylor  factors  for  some  of  the  orientations  in  Table  6-2  were  pre¬ 
viously  reported  by  Chin  et  al  [6.8J.  The  present  results  are  identical  to 
those  reported  by  Chin  et  al.  Also  included  in  Table  6-2  are  some  of  the 
"weakest"  orientations  (0^,  O2,  and  O3)  which  result  in  low  M  values.  The 
lowest  possible  value  for  M  is  2,  and  it  is  observed  in  the  grain  orientation 
0^  The  M  values  for  orientations  02  (and  O3)  and  the  [0011  orientation  are 
2.26  and  2.45,  respectively.  These  "weakest"  orientations,  summarized  in  the 
standard  stereographic  projection  shown  in  Figure  6-11,  are  unique  because  the 
plane  strain  condition  can  be  satisfied  by  activating  the  highest  stressed 
slip  system  only.  The  predicted  "weakest"  orientations  are  compared  in 
Figure  6-12  with  the  orientations  of  grains  within  which  small  cracks  were 
initiated.  The  results  indicate  that  most  of  the  observed  orientations  are 
within  5  degrees  of  the  "weakest"  orientations;  the  worst  cases  are  within 
10  degrees. 

A  conclusion  which  can  be  drawn  from  the  results  shown  in  Figure  6-12  is 
that  slip  in  coarse-grained  Astroloy,  and  thus  fatigue  cracks,  tends  to  ini¬ 
tiate  within  grains  having  the  "weakest"  orientations.  This  seems  to  be  true 
for  slip  initiating  from  pores  or  within  the  grains.  This  behavior  can  be 
easily  understood  since  in  the  present  experiments  the  applied  maximum  stress 
was  above  the  proportional  limit,  but  below  the  macroscopic  cyclic  yield 
stress  (S  =  0.8),  thus  only  grains  with  low  Taylor  factors  can  exceed  the 
critical  resolved  shear  stress  required  to  activate  slip.  Once  nucleated 
within  slip  bands,  small  fatigue  cracks  propagate  in  a  microplastic  field 
since  they  are  entirely  embedded  within  slip  bands  which  extend  across 
individual  grains.  Fatigue  crack  propagation  can  cease  (either  permanently  or 
temporary)  or  continue  as  the  microcracks  approach  the  grain  boundaries, 
depending  on  whether  or  not  the  crack  tip  can  propagate  slip  onto  contiguous 
grains. 

In  order  to  gain  insight  into  the  effect  of  statistical  variation  of  the 
yield  strength  of  individual  grains  on  crack  initiation,  Taylor  factors  were 
computed  for  grains  whose  orientations  were  selected  randomly.  The  number  of 
grains  was  varied  from  one  to  twenty,  and  each  calculation  was  repeated  for 
fifty  sets  of  randomly  selected  grain  orientations.  Results  of  M  values  nor¬ 
malized  by  the  average  M  value,  M,  for  all  possible  orientations,  shown  as  a 
function  of  number  of  grains  in  Figure  6-13,  indicate  that  a  large  statistical 
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12.  Standard  stereographic  triangle  showing 
orientations  with  low  Taylor  factor 
values. 
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FIGURE  6-13.  Variation  of  normalized  Taylor  factors 
(M/M)  with  the  number  of  grains. 


variation  in  yield  strength  is  expected  in  the  plastic  zone  when  only  one 
grain  is  interrogated  by  the  crack  tip.  This  statistical  variation  of  yield 
strength  decreases  as  the  number  of  grains  affected  by  the  crack  tip  is 
increased.  Thus,  the  effect  of  grain  orientation  is  expected  to  play  a 
prominent  role  in  crack  initiation,  as  well  as  in  the  subsequent  erratic 
growth  of  small  fatigue  cracks. 

In  conclusion,  fatigue  cracks  in  Astroloy  at  room  temperature  initiate 
preferentially  in  grains  which  were  weakest  by  virtue  of  the  fact  that  their 
crystallographic  orientations  favor  {111}  <110>  slip.  These  grains  were  found 
to  have  surface  normals  near  the  ( 100 }  [110]  and  [113]  directions  and  exhibit 
the  lowest  Taylor  factors.  This  preferred  orientation  was  also  observed  for 
the  one  out  of  seven  cases  (in  the  current  experiment)  where  a  crack  initiated 
from  a  pre-existing  surface  pore.*  Thus,  crystallographic  orientation  is  an 
important  factor  in  governing  where  the  crack  initiates,  as  well  as  its  mode 
of  growth.  However,  it  may  also  be  possible  for  metallurgical  defects  to 
initiate  cracks  in  less  favorably  oriented  grains,  depending  on  the  nature, 
geometry,  size  and  distribution  of  such  defects.  Manufacturing  defects  would 
also  be  expected  to  play  a  similar  ro’e  in  crack  initiation  in  actual  aero¬ 
engine  components  as  discussed  in  Section  5.4. 

It  is  significant  to  note  that  the  crystallographic  orientation  deter¬ 
mined  here  for  small  cracks  in  Astroloy  is  the  same  as  those  observed  near  the 
fatigue  threshold  in  long  cracks  ;n  Astroloy  using  fracture  surface  etching 
[6.9  -  6.11],  as  well  as  those  in  other  Ni-base  superalloys  determined  using  a 
variety  of  techniques  [6.1?  -  6.14], 

6.3  Fractoqraphy  of  Small  Cracks 

The  fracture  surface  morphology  of  small  cracks  was  characterized  as  a 
function  of  test  temperature,  alloy  composition  and  grain  size  by  examining 
the  fractured  bend  specimens  from  Section  5.  These  results  assist  in  under¬ 
standing  the  influence  of  these  variables  on  the  small  crack  growth  rates. 
More  importantly,  however,  the  results  are  compared  with  known  large  crack 
behavior  to  elucidate  possible  differences  in  the  crack  growth  mechanisms  in 
small  and  large  cracks. 


6.3.1  Fractography  Versus  Test  Temperature 


Test  temperature  had  a  pronounced  effect  on  the  fracture  surface 
morphology  of  small  cracks  as  indicated  by  the  results  in  Figure  6-14  on 
CG  Astroloy.  The  initial  fracture  surfaces  at  25°C  were  highly  crystal¬ 
lographic  exhibiting  sharp,  angular  facets.  In  contrast,  crack  growth  at 
600°C,  while  still  initially  crystallographic,  is  significantly  flatter  and 
smoother  than  that  at  25°C.  Consequently,  the  relationship  between  facet  size 
and  grain  size,  which  is  clearly  evident  at  25°C,  is  indistinguishable  at  the 
elevated  temperature.  At  both  temperatures,  the  extent  of  crystallographic 
cracking  decreased  as  the  cracks  enlarged,  although  a  transition  crack  size 
can  be  identified  in  each  case.  This  transition  crack  size  was  larger  at  room 
temperature  than  at  600°C  (a  =  400  um  versus  200  pm). 

The  above  change  in  the  fracture  surface  morphology  with  temper¬ 
ature  also  occurred  in  Waspaloy.  However,  as  indicated  in  Figure  6-15  the  ex¬ 
tent  of  the  sharp,  angular  faceting  at  25°C  was  noticeable  less  in  Waspaloy 
than  in  CG  Astroloy.  In  accord  with  subsequent  discussion,  this  difference  is 
believed  to  be  due  to  the  fact  that  the  room  temperature  Waspaloy  specimen  was 
tested  at  a  significantly  higher  level  of  applied  stress  (S  =  1.7  versus  1.1). 

The  dramatic  influence  of  test  temperature  on  the  nature  of  the  crystal¬ 
lographic  crack  growth  was  particularly  evident  in  cases  where  the  temperature 
was  changed  during  the  experiment.  This  is  shown  in  Figure  6-16  which  con¬ 
tains  higher  magnification  fractographs  of  the  crack  initiation  sites  in 
CG  Astroloy  and  Waspaloy.  In  both  cases  cracks  were  initiated  at  25°C  until  a 
crack  of  about  one-to-two  grain  sizes  was  detected,  then  the  temperature  was 
increased  to  600°C  for  the  remainder  of  the  experiment.  The  transition  from 
the  sharp  angular  facets  at  25°C  to  the  flatter,  smoother  facets  at  600°C  is 
clearly  evident.  Interestingly,  this  transition  corresponds  to  the  point  in 
the  experiment  where  the  cracks  often  arrested  for  a  significant  period  of 
time,  as  described  in  Section  5.4. 

The  above  changes  appear  to  be  consistent  with  observed  changes  in  the 
homogeneity  of  slip  in  Ni-based  alloys  with  temperature.  Specifically,  in¬ 
creasing  temperature,  as  well  as  strain,  is  known  to  activate  processes  such 
as  cross  slip  and  climb  which  thereby  serve  to  homogenize  the  deformation 
[6.16  -  6.20| .  Such  thermally  activated  processes  would  be  expected  to 


F I GURE  6-14.  Effect  of  test  temperature  on  fracture  surfaces  of  small  cracks 
in  CG  Astroloy. 


Effect  of  test  temperature  on  fracture  surfaces  of  small  cracks 
in  Waspaloy  versus  CG  Astroloy. 


6.3.2  Fractoqraphy  Versus  Alloy  Composition 

Figures  6-15  and  6-16  also  provide  a  comparison  of  the  influence 
of  alloy  composition  on  the  fracture  surface  morphology  since  the  grain  sizes 
are  similar  in  both  CG  Astroloy  (30  urn)  and  Waspaloy  (42  ym).  Fracture 
surfaces  in  both  alloys  are  strikingly  similar,  exhibiting  the  same 
temperature  dependence  as  well  as  change  from  faceted  to  non-facetted  crack 
growth  with  increasing  crack  size.  The  only  noticeable  difference  in  fracture 
appearance  is  that  the  facets  in  Waspaloy  appear  more  ductile  than  those  in  CG 
Astroloy,  Figure  6-16.  This  difference  Is  understandable,  since  Waspaloy  is 
the  weaker  of  the  two  alloys,  the  enhanced  strength  of  Astroloy  being  gained 
through  additional  precipitation  hardening. 

The  overall  commonality  of  fracture  appearance  of  Waspaloy  and  Astroloy 
is  consistent  with  the  fact  that  the  growth  rate  behavior  of  small  cracks  in 
these  two  alloys  is  similar  at  both  room  temperature  [6.21]  and  elevated 
temperature  (Figure  5-22). 

6.3.3  Fractoqraphy  Versus  Grain  Size 

As  shown  in  Figure  6-17  for  the  case  of  Astroloy,  the  primary 
role  of  grain  size  on  the  fracture  surface  morphology  is  to  alter  the  size  of 
the  crystallographic  facets.  This  is  clearly  evident  at  25°C,  but  less  so  at 
600°C  since  the  smooth,  flat  surfaces  produced  at  this  temperature  largely 
obscure  the  grain  boundaries.  Nevertheless,  as  the  crack  enlarges  and  non¬ 
crystal  lographic  cracking  becomes  more  predominant  at  600°C,  the  fracture 
surface  of  the  fine  grain  material  exhibits  a  noticeable  finer  structure  than 
that  of  the  coarse  grained  material. 

6.3.4  Fractoqraphy  of  Small  Versus  Large  Cracks 

The  above  mentioned  trends  in  fracture  morphology  with 
temperature,  alloy  composition  and  grain  size  are  similar  to  those  which  have 
previously  been  observed  for  long,  thru-thickness  cracks  in  Ni-base  alloys  in 
the  near-threshold  regime  [6.9  -6.111.  An  example  of  the  latter,  taken  from 
Ref.  [6.11],  is  shown  in  Figure  6-18  for  both  FG  and  CG  Astroloy  at  25 J C  and 
600  C.  Comparing  Figure  6-17  with  Figure  6-16,  and  accounting  for  the 


Effect  of  grain  size  and  temperature  on  the  fractoaraphy  of  small 
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FIGURE  6-18.  Effect  of  grain  size  and  temperature  on  the 
cracks  in  Astroloy  (after  Ref.  6.14).  Note 
crack  fractography  in  Figure  6-17. 


dependence  of  the  fracture  surface  Is  similar  in  both  crack  size  regimes. 
This  similarity  Is  particularly  noticeable  in  the  coarse  grained  material. 


As  discussed  previously  in  relation  to  the  observed  temperature 
dependence,  the  fracture  surface  morphoTogy  should  be  related  to  the  homo¬ 
geneity  of  slip,  which  is  known  to  change  with  both  temperature  and  extent  of 
deformation.  Thus,  a  more  exact  comparison  of  the  fracture  surfaces  of  small 
and  long  cracks  should  be  done  at  equivalent  values  of  the  mechanical  driving 
force  --  that  is,  at  equivalent  magnitudes  of  the  local  crack-tip  strain. 

In  long,  thru-thickness  cracks  In  Ni-base  superalloys  it  is  well 
documented  that  crystallographic  faceting  is  favored  at  low  a K,  while  non¬ 
crystal  lographic  (classical  striated)  growth  is  favored  at  high  aK  [6.9, 
6.221.  In  fact,  it  has  been  proposed  that  the  transition  between  the  two 
crack  growth  processes  occurs  when  the  cyclic  plastic  zone  becomes  equal  to 
the  grain  size  [6.221 .  However,  is  should  be  pointed  out  that  these  studies 
have  computed  the  cyclic  plastic  zone  size  based  on  the  applied  AK  value.  This 
procedure  is  likely  to  have  over-estimated  the  actual  cyclic  plastic  zone  size 
due  to  the  significant  crack  closure  which  occurs  In  large  cracks,  as  shown  in 
Section  4.4.  As  shown  in  Section  7,  small  cracks  also  exhibit  closure,  but  to 
a  much  lesser  extent  than  do  long  cracks.  This  difference  in  crack  closure 
behavior  together  with  inelastic  and  mixed-mode  effects  in  small  cracks, 
causes  the  mechanical  driving  force  for  small  and  large  cracks  to  differ. 
Thus,  a  more  sensible  comparison  of  fracture  features  should  be  made  in  terms 
of  an  equivalent  aK,  (defined  in  Section  8),  which  takes  these  factors  into 
account  for  both  small  land  large  cracks. 

Figure  6-19  illustrates  the  general  relationship  between  the 
applied  aK  and  the  equivalent  aK,  designated  AKeq.  As  indicated,  AKeq  is  less 
than  aK  for  long  cracks  near  the  threshold,  while  the  opposite  is  true  for 
small  cracks.  Consequently,  the  respective  AKeq  values  for  small  and  long 
cracks  are  much  more  similar  than  indicated  by  their  applied  aK  values.  This 
similarity  in  AKeq,  which  if  properly  formulated  should  characterize  the  near¬ 
tip  deformation  state  in  both  crack  regimes,  is  believed  to  account  for  the 
similarities  in  crack  growth  mechanisms  in  small  and  long  cracks. 

In  addition,  differences  in  constraint  from  neighboring  grains 
in  small  versus  large  cracks  may  favor  crystallographic  crack  growth  in  small 
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FIGURE  6-19.  Schematic  illustration  of  consolidation  of  small 
and  large  crack  growth  rate  data  in  terms  of 

AK^q . 


the  first  few  grains  they  can  deform  and  crack  in  an  optimum  manner  dictated 
by  the  particular  orientation  of  each  grain  (Section  6.2).  In  contrast,  both 
large  surface  cracks  and  long  thru  cracks  are  constrained  to  deform  and  crack 
in  concert  with  a  large  number  of  grains  thereby  favoring  slip  on  multiple 
planes  which  gives  rise  to  a  flatter  fracture  surface  having  fewer  sharp 
facets. 

In  conclusion,  the  fractography  of  small  cracks  is  consistent 
with  the  orientation  relationship  established  In  Section  6.2.  Perhaps  more 
importantly,  both  these  observations  are  consistent  with  known  crack  propaga¬ 
tion  mechanisms  in  long  cracks.  Thus,  the  resolution  of  the  "anomalous"  be¬ 
havior  of  small  cracks  appears  to  require  a  proper  reformulation  of  the 
mec./anical  driving  force  --  a  task  undertaken  in  Sections  7  and  8. 
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7.  CRACK-TIP  DEFORMATION  AND  CRACK  CLOSURE  IN  SMALL  CRACKS 


The  work  summarized  in  this  section  was  designed  to  elucidate  the 
physical  processes  which  contribute  to  crack  size  effects  in  Ni-base  alloys 
under  fatigue  loading.  This  was  accomplished  using  high  resolution  observa¬ 
tions  and  quantitative  measurements  obtained  within  a  scanning  electron  micro¬ 
scope  (SEM)  equipped  with  a  special  fatigue  loading  stage.  The  local  dis¬ 
placements  at  the  tips  of  small  cracks  were  quantified  by  combining  the  SEM 
observations  with  the  techniques  of  stereoimaging.  These  local  displacement 
measurements,  combined  with  crack  closure  measurements,  provided  the  basis  for 
the  crack  driving  force  determined  in  Section  8.  While  most  of  the  data  pre¬ 
sented  in  this  section  are  on  the  Ni-base  superalloys  of  primary  interest  in 
this  study,  supplemental  results  obtained  from  reanalysis  of  previously  avail¬ 
able  data  on  an  aluminum  alloy  are  also  presented. 

7.1  Experimental  Techniques 

The  local  deformation  and  closure  of  small  fatigue  cracks  was  measured 
within  a  SEM  using  either  of  two  cyclic  loading  stages  capable  of  operating  at 
either  room  temperature  or  elevated  temperature.  Details  of  the  design  and 
operation  of  these  units  can  be  found  in  Refs.  7.1  and  7.2. 

The  specimen  design  used  for  the  SEM  experiments  is  shown  in  Figure  7- 
1.  All  specimens  were  prepared  using  low-stress  machining.  In  addition,  to 
eliminate  any  residual  stresses  remaining  after  machining,  the  gage  section  of 
the  specimens  were  further  prepared  by  electropolishing  in  a  solution  of 
25  percent  nitric  acid  in  methanol  at  -30°C.  Typically,  this  procedure  re¬ 
moved  a  40  um-thick  layer  of  material  from  the  specimen.  This  procedure  also 
revealed  the  microstructure  of  the  material  thereby  providing  a  frame  of 
reference  for  the  observation  of  small  cracks,  as  well  as  surface  features 
necessary  for  stereoimaging. 

7.1.1  Initiation  of  Fatigue  Microcracks 

All  fatigue  cracks  were  naturally  initiated  by  loading  the  gage 
section  of  the  specimen  in  Figure  7-1  in  three-point  bending.  As  in  the 
experiments  of  Section  5,  this  procedure  was  used  to  localize  the  initiated 


FIGURE  7-1.  Specimen  design  used  to  study  small  cracks 
within  the  SEM.  Dimensions  are  in  mm. 


cracks  thereby  simplifying  their  detection.  For  expedience,  crack  initiation 
was  conducted  in  a  room  temperature  laboratory  air  environment  (50  percent 
R.  H.)  using  a  laboratory  servo-hydraul ic  test  machine  cycled  under  a  10  Hz 
sinusoidal  waveform  and  a  load  ratio  (R)  of  0.1.  The  maximum  (outer  fiber) 
bending  stresses  were  85  to  98  percent  of  the  proportional  limit  of  the 
material,  as  determined  from  strain  gage  measurements. 

Crack  initiation  was  detected  using  plastic  replicas  obtained 
periodically  under  both  minimum  and  maximum  load.  Replicas  were  metal  coated 
and  examined  using  optical  microscopy  at  500X.  This  procedure  usually  re¬ 
vealed  either  no  cracks  or  several  cracks;  rarely  did  only  a  single  crack  ini¬ 
tiate.  However,  multiple  cracks  were  usually  widely  separated  and  did  not 
interfere  with  one  another  until  considerable  growth  had  occurred.  As 
described  in  Section  6,  these  cracks  were  found  to  initiate  either  at  slip 
within  grains  or  at  surface  pores.  Occasionally,  the  slip  which  initiated  the 
crack  could  be  traced  back  to  the  very  first  loading  cycle. 

7.1.2  SEM  Experiments 

Following  the  initiation  of  a  detectable  crack,  the  specimens 
were  switched  to  tension-tension  loading  (R  =  0.1)  in  vacuum  --  either  in  a 
laboratory  test  machine  equipped  with  a  vacuum  rig,  or  within  the  cyclic 
loading  stage  of  the  SEM.  Ac  room  temperature  this  transition  occurred  with¬ 
out  noticeable  change  in  the  crack  growth  rate.  However,  when  this  transition 
involved  increasing  the  test  temperature  to  SOOT,  cracks  arrested  for  many 
thousands  of  cycles,  then  eventually  resumed  growth  without  warning.  This 
behavior  was  identical  to  that  which  was  observed  in  the  experiments  of 
Section  5  and  greatly  hindered  data  acquisition  at  the  elevated  temperature. 

To  overcome  this  difficulty,  the  initial  portion  of  the  600'C 
experiments  had  to  be  conducted  entirely  within  the  SEM  where  cracks  could  be 
constantly  monitored  at  high  resolution.  During  this  period,  the  cyclic  load 
range  was  gradually  increased  until  the  cracks  remained  open  for  a  significant 
portion  of  the  loading  cycle  and  crack  growth  was  observed.  This  procedure 
typically  resulted  in  a  30  percent  increase  in  the  stress,  relative  to  the 
precracking  stress.  During  this  procedure,  the  cyclic  frequency  was  main¬ 
tained  at  3  Hz  in  order  to  preclude  measurable  creep  deformation.  Once  cracks 


were  growing  at  600“C,  the  specimen  was  periodically  transferred  from  the  SEM 
to  the  laboratory  test  machine  (at  1  MPa  vacuum)  to  lengthen  cracks  prior  to 
reinsertion  into  the  SEM  for  further  observation. 

7.1.3  Stereoimaqlnq  Analysis 

Periodically,  during  growth  of  the  small  cracks,  photographs  of 
the  entire  crack  and  each  of  the  crack  tip  regions  were  made  at  both  minimum 
and  maximum  load.  Examples  of  these  photographs  are  shown  in  Figures  7-2  and 
7-3.  Cracks  in  FG  Astroloy  initiated  primarily  at  pores.  Figure  7-2,  whereas 
cracks  in  CG  Astroloy  initiated  primarily  at  persistent  slip  bands.  Figure  7- 
3. 

Photographs  were  analyzed  using  the  stereoimaging  technique 
(7.2),  which  is  an  optical  method  of  measuring  relative  displacements  by  com¬ 
paring  two  photographs  of  the  same  region,  each  obtained  at  a  different  dis¬ 
placement  state.  Using  this  procedure,  the  crack-tip  displacement  field  was 
quantified  using  high  resolution  photogrammetric  equipment.  From  these  mea¬ 
surements,  it  was  a  simple  matter  to  compute  crack  opening  displacements. 
Furthermore,  by  taking  gradients  within  the  displacement  field  three  elements 
of  the  symmetric,  in-plane  strain  tensor  were  also  computed  using  finite- 
strain  formulation.  These  strains  were  used  to  determine  the  magnitude  of 
crack-tip  strain  and  the  distributions  of  strain  parallel  and  perpendicular  to 
the  loading  direction  (7.3,  7.4).  Plastic  zone  sizes  were  also  computed  by 
extending  these  distribution  functions  to  the  cyclic  elastic  limit,  defined  as 
twice  the  cyclic  yield  stress  divided  by  the  elastic  modulus. 

The  stereoimaging  technique  was  also  used  to  determine  the  load 
at  which  the  crack  was  fully  opened  to  its  tip.  This  procedure  entailed 
photographing  the  crack-tip  region  at  minimum  load,  Pmin,  and  at  progressively 
higher  proportions  of  the  maximum  load,  Pmax*  By  first  observing  the 
photographs  of  Pmin  and  Pmax  in  a  stereoviewer,  it  was  possible  to  determine 
the  location  of  the  crack  tip  to  within  about  0.1  urn.  At  loads  below  that 
required  for  crack  opening  to  the  tip,  only  a  portion  of  the  crack  appeared 
open,  and  with  increasing  load,  this  distance  decreased  until  the  crack  was 
completely  open  to  the  tip.  The  load  at  which  this  occurred  was  defined  as 


FIGURE  7-3.  14  um  long  crack  in  Coarse  Grained  Astroloy  at  20°C.  Crack 

was  slip  band  initiated. 


It  was  usually  possible  to  determine  the  value 


the  crack  opening  load,  PQp. 
of  Pop  to  10-15  percent  of  Pmax. 

By  viewing  photographs  both  parallel  and  perpendicular  to  the 
loading  axis,  a  value  of  applied  load  corresponding  to  the  onset  of  Mode  I 
opening  (P0p),  as  well  as  the  onset  of  Mode  II  sliding  (Ps)  could  be  deter¬ 
mined.  Occasionally,  for  large  cracks  and  low  aK  values,  crack-tip  sliding 
was  found  to  precede  crack-tip  opening  (that  is,  P$  <  Pop);  however,  for  small 
cracks  this  was  not  the  case.  In  selected  cases,  Mode  III  displacements  were 
also  determined  by  comparing  photographs  of  the  specimen  inclined  at  two  dif¬ 
ferent  angles  with  respect  to  the  electron  beam. 

7.2  Local  Crack-Tip  Measurements  on  Small  Cracks 

The  crack  tip  displacements  strains  and  crack  opening  loads  were  deter¬ 
mined  for  small  surface  cracks  of  various  sizes  in  CG  Astroloy  at  25°C  and  in 
CG  and  FG  Astroloy  at  600°C.  Oetailed  results  are  presented  for  selected 
cases  to  illustrate  the  nature  of  the  results  obtained  and  summaries  of 
important  parameters  are  provided  for  the  remainder  of  cases. 

7.2.1  Crack-Tip  Displacements 

Figure  7-4  provides  an  example  of  a  typical  displacement  field 
determined  from  stereoimaging.  These  results  were  obtained  for  the  small 
crack  in  FG  Astroloy  at  600°C  shown  in  Figure  7-2.  All  displacements  are  with 
respect  to  an  reference  point  ahead  of  the  crack  tip  at  a  remote  location 
which  experienced  little  or  no  deformation.  These  displacements  are  indica¬ 
tive  of  a  biaxial  crack-tip  deformation  field  and  are  surprisingly  homogeneous 
in  view  of  the  highly  planar  slip  character  of  this  alloy.  Similar  results 
were  observed  in  CG  Astroloy  and  both  25°C  and  600“C. 

As  indicated  in  Section  6,  small  cracks  were  rarely  normal  to  the 
applied  loading  direction  due  to  their  preferred  orientation  within  the 
crystallite.  Thus,  cracks  were  generally  inclined  at  some  angle  with  respect 
to  the  original  x'-y'  coordinate  system  of  the  stereoimaging  measurements,  as 
shown  in  Figure  7-4.  In  order  to  resolve  displacements  into  Modes  I  and  II, 
the  measured  displacements  were  transformed  from  the  original  x'-y'  system 
into  an  x-y  system  as  shown  in  Figure  7-5. 
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FIGURE  7-5.  Original  ( X ’  -  Y‘)  and  transformed  (X-Y) 
coordinates  for  the  measured  displacement 
fields  at  small  cracks. 


An  example  of  Mode  I  and  Mode  II  crack  opening  displacements 
(COO)  in  the  x-y  coordinate  system  Is  given  in  Figure  7-6.  Once  again,  these 
results  are  for  the  small  crack  in  FG  Astroloy  shown  in  Figure  7-2  and  cor¬ 
respond  to  the  original  displacement  field  of  Figure  7-4.  The  magnitude  of 
the  Mode  II  displacements  in  this  particular  case  is  relatively  small  since 
the  crack  was  nearly  normal  to  the  loading  direction.  In  contrast,  the 
Mode  II  displacements  are  significant  for  a  crack,  such  as  that  shown  pre¬ 
viously  in  Figure  7-3,  which  is  inclined  at  an  appreciable  angle  to  the 
loading  axis.  For  example,  cracks  inclined  at  an  angle  of  40°-50°  to  the 
loading  axis  exhibited  nearly  equal  Mode  I  and  Mode  II  displacements,  as  shown 
in  Figure  7-7  for  a  crack  in  CG  Astroloy  at  600°C.  This  mixture  of  Modes  I 
and  II  is  in  good  agreement  with  results  obtained  from  two  dimensional  elastic 
analysis  of  inclined  edge  cracks  using  the  boundary  force  methods  [7.5] 

Perhaps  the  most  Interesting  feature  of  the  above  COD  values  is 
the  fact  that  they  are  a  linear  function  of  square  root  of  the  distance  from 
the  crack  tip  and  thus  can  be  described  by: 

COD  =  CQ/r  (7-1) 

where  CQ  is  the  slope  of  the  lines  in  Figures  7-6  and  7-7,  and  r  is  the  dis¬ 
tance  from  the  crack  tip.  The  above  r-dependence  of  the  COD  is  precisely  what 
one  would  expect  for  an  elastic  crack  where,  for  example,  the  Mode  I  displace¬ 
ment  ranges  are  given  by: 

COO  =  [8(l+v)(l-v)|  (7-2) 
E/2tt 

where  aK  is  the  Mode  I  stress  intensity  factor  defined  by  linear  elasticity 
and  v  is  Poisson's  ratio.  By  combining  Eqs.  (7-1)  and  (7-2),  it  would  appear 
that  aK  can  be  computed  from  the  slope  of  lines  found  in  Figures  7-6  and  7-7, 
provided  the  cracks  conform  to  the  linear  elastic  fracture  mechanics  concept 
of  small  scale  yielding  (SSV).  However,  based  on  the  size  of  the  measured 
plastic  zone  in  relation  to  the  crack  size,  this  is  not  likely  to  be  the  case 
for  small  cracks  as  described  in  Section  7.2.2.  Thus,  the  /r  dependence  of 
COD  appears  to  be  fortuitous. 
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FIGURE  7-6.  Crack  opening  displacements  for  FG  Astroloy 
at  600°C.  Upper  and  lower  figures  corre- 
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FIGURE  7-7.  Crack  opening  displacements  in  CG  Astroloy  at  600°C. 

Large  Mode  II  displacements  are  due  to  the  fact  that 
the  crack  is  Inclined  with  respect  to  the  loading 
axis. 


In  order  to  better  understand  the  above  observations,  Eqs.  (7-1) 
and  (7-2)  were  applied  to  the  analysis  of  large  cracks  in  7075  and  7091 
aluminum  alloys  which  conformed  to  the  SSV  condition.  Displacements  measured 
using  stereoimaging  were  again  found  to  exhibit  the  /r-dependence.  However, 
the  aK  values  computed  from  the  measured  displacements,  termed  aK(disp.) 
differed  from  those  computed  from  applied  load  and  crack  size  using  the 
standard  elastic  stress  intensity  factor  equations.  As  shown  in  Figure  7-8, 
the  difference  between  aK(disp.)  and  the  applied  aK  was  found  to  increase  with 
increasing  applied  aK.  The  reason  for  this  discrepancy  is  thought  to  be  re¬ 
lated  to  the  influence  of  the  fatigue  crack  extension  process,  and  associated 
crack-tip  plasticity  and  crack  closure,  on  the  near-tip  displacement  field. 
Thus,  even  when  SSY  prevails,  the  local  displacements  at  growing  fatigue 
cracks  do  not  conform  to  the  analytical  results  from  elastic  analysis  of 
monotonical ly  loaded,  stationary  cracks.  Nevertheless,  when  SSY  prevails, 
these  local  events  appear  to  scale  with  aK,  as  evidenced  by  the  ability  of  aK 
to  provide  a  geometry- independent  measure  of  the  crack  growth  kinetics  of 
large  cracks,  at  least  under  constant  amplitude  loading. 

7.2. 2  Crack-Tip  Strain  and  Plastic  Zone  Size 

An  example  of  the  maximum  shear  strain  distributions  at  the  tips 
of  a  small  crack  in  FG  Astroloy  at  600°C  (Figure  7-2)  are  given  in  Figure  7- 
9.  The  strain  field  is  surprisingly  homogeneous  in  view  of  the  slip  character 
of  this  material.  Thus,  multiple  slip  must  be  occurring  on  numerous  (111) 
planes  within  the  crack-tip  plastic  zone.  Although  stereoimaging  is  capable 
of  defining  the  average  crack-tip  strain  field,  it  cannot  resolve  slip  details 
which  occur  on  a  much  finer  scale.  Thus,  it  is  not  surprising  that  the 
criterion  controlling  which  of  these  (111)  slip  planes  the  crack  will  follow 
is  not  obvious  from  examining  the  strain  field.  Similar  strain  fields  were 
observed  in  CG  Astroloy  at  both  25°C  and  600°C  and  gave  no  indication  of  the 
change  in  fracture  morphology  with  temperature  which  was  evident  in  the  frac- 
tographic  results  of  Section  6.3. 

These  strain  distributions  for  small  cracks  in  Ni-base  super¬ 
alloys  are  similar  to  those  found  for  small  cracks  in  aluminum  alloys  (7.6), 
in  that  there  is  a  large  component  of  strain  roughly  parallel  to  the  direction 
of  loading.  This  distribution  of  strain  is  different  than  that  found  for 
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FIGURE  7-8.  Comparison  of  AK  determined  from  measured  crack 
displacements  with  applied  aK  for  large  cracks 
in  several  aluminum  alloys. 


large  cracks  (7.7,  7.8)  where  a  higher  level  of  strain  is  found  ahead  of  the 
crack  tip.  These  differences  in  strain  distribution  between  large  and  small 
fatigue  cracks  appear  to  depend  on  whether  or  not  SSY  prevails  and  are  be¬ 
lieved  to  be  fundamental  to  understanding  the  small  crack  problem.  Further¬ 
more,  the  similarity  in  strain  distribution  between  small  cracks  in  super¬ 
alloys  and  in  aluminum  alloys  indicates  the  commonality  in  crack-tip  mechanics 
for  these  two  materials. 

Strain  distributions  for  large  fatigue  cracks  have  been  found  to 
be  of  the  following  functional  form: 

1)  Perpendicular  to  the  loading  axis  ($  =  0°)  (7.7] 

Ae'(r)  =  A'  +  M'  in  (r+B')  (7-3) 

or 

Ae'(r)  =  — - -  and  m  =  l/(l+n')  (7-4) 

(r+B)m 

2)  Parallel  to  the  loading  axis  ($  =  90°)  [7.8]: 

“'<r>  -  7PFT  (7-5) 

where  r  =  distance  from  the  crack  tip,  Ae'(r)  =  (strain  value  at  r)/(strain 
value  at  r  =  0),  n1  =  the  cyclic  hardening  exponent  given  in  Section  4,  and 
A',  M',  8',  A  and  B  are  material  constants.  Considering  the  variability  in 
the  measured  strain  data,  it  is  not  possible  to  determine  whether  Eq. (7-3)  or 
Eq. (7-4)  provides  the  better  representation  of  the  data. 

For  small  cracks  in  both  CG  and  FG  Astroloy  at  25°C  and  600  C, 
Eq . (7-4)  has  been  found  to  fit  the  data  reasonably  well,  but  it  is  likely  that 
Eq. (7-3)  would  also  fit  since  these  two  functions  are  very  similar  near  the 
crack  tip.  The  validity  of  Eq.  (7-5)  for  small  cracks  has  not  been  examined. 

Crack-tip  plastic  zone  sizes  (rp)  for  0  =  0°  and  90’  were 

computed  by  fitting  the  measured  strain  data  for  small  cracks  in  Astroloy  to 
the  above  equations.  These  results  are  summarized  in  Table  7-1  along  with  the 
crack  sizes,  applied  aK  values  and  crack-tip  strain  ranges,  At(0).  Although 
the  validity  of  Eq.  (7-5)  has  not  been  verified  for  the  specific  case  of  small 
cracks  in  Astroloy,  it  has  been  assumed  to  be  valid  for  the  purposes  of 
comparing  plastic  zone  sizes.  The  plastic  zone  sizes  for  *  =  0  and  90  are 


TABLE  7-1 


CRACK-TIP  STRAINS  AND  PLASTIC  ZONE  SIZES 
FOR  SMALL  CRACKS  IN  ASTROLOY 


Crack 

*K  _ 

Crack-Tip 

V 

U* 

MPa/m 

Strain,  Ae(0) 

»  =  0° 

4>  =  90° 

C6  Astroloy,  25°C 

45 

4.9 

0.042 

5.6 

4.0 

5.9 

0.032 

13.1 

3.0 

67 

5.9 

0.120 

15.3 

14.0 

0.050 

10.4 

5.0 

71 

5.4 

0.0046 

7.1 

- 

0.0054 

7.9 

- 

72 

6.2 

0.016 

4.9 

1.0 

0.049 

10.4 

5.0 

102 

6.5 

0.068 

21.1 

7.0 

CG  Astroloy,  600°C 

42 

4.4  0.038 

5.9 

3.0 

0.029 

6.6 

2.0 

54 

5.1  0.060 

7.6 

5.0 

0.127 

12.1 

11.0 

66 

5.4  0.107 

30.5 

9.0 

0.115 

17.0 

10.0 

FG  Astroloy,  600°C 

106 

7.7  0.046 

20 

3.0 

0.057 

25 

4.0 

133 

8.8  0.093 

24 

8.0 

0.043 

28 

3.0 

230 

11.9  0  036 

29 

3.0 

0.071 

43 

6.0 

plotted  as  a  function  of  aK  in  Figures  7-10  and  7-11,  respectively.  In 
neither  case  does  rp  follow  the  parabolic  relationship  with  aK  expected  from 
analytical  SSY  results.  The  large  variation  in  rp  versus  aK,  particularly 
for  $  =  90°  (Figure  7-11)  occurs  because  of  the  lack  of  symmetry  observed  in 
the  deformation  field  of  actual  cracks. 

In  order  to  assess  the  extent  to  which  the  SSY  concept  applies  to 
these  small  cracks,  the  ratio  of  crack  size  to  plastic  zone  size  (a/rp) 
for  0=0  was  examined.  These  results,  given  as  a  function  of  crack  size  in 
Figure  7-12,  show  that  a/rp  was  typically  2  to  4,  but  overall  ranged  from  1  to 
8.  These  values  are  considerably  smaller  that  of  a/rp  >  50  which  is  generally 
thought  to  be  reguired  to  insure  SSY  [7.9).  Consequently,  an  elastic-plastic 
crack-tip  parameter  is  needed  to  adequately  describe  the  driving  force  for  the 
growth  of  small  cracks.  Such  a  parameter  is  considered  in  Section  8. 

7.2.3  Crack  Closure  in  Small  Versus  Large  Cracks 

The  stereoimaging  technique  was  also  used  to  measure  the  load 
required  to  fully  open  the  tip  of  small  cracks,  Pop.  These  values  were  used 
to  determine  relative  effective  stress  intensity  factors,  AKeff/AK,  using  the 
following  relationship  derivable  from  basic  crack  closure  definitions  (see 
Section  4.4) : 

AKeff/AK  =  (1  -  Pop/Pmax)/(l  -  R)  (7-6) 

Results  for  small  cracks  are  shown  in  Figure  7-13,  along  with  the  Soniak  and 
Remy  (7.10)  data  from  large  cracks  which  were  previously  discussed  in  Section 
4.4.  As  indicated,  values  of  A«eff/AK  for  small  cracks  are  relatively  con¬ 
stant  over  a  wide  range  of  Kmax  (or  aK)  and  are  not  measurably  sensitive  to 
material  condition  or  test  temperature.  These  small  crack  results  are 
described  by  the  following  simple  relationship: 

AKgf f /aK  =  6  (7-7) 

where  e  =  0.35  for  both  FG  and  CG  Astroloy  at  25°C  and  600°C.  In  contrast,  as 
indicated  in  Figure  7-13,  the  long  crack  results  in  Astroloy  are  sensitive  to 
both  Kmax  and  test  temperature. 
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FIGURE  7-10.  Cyclic  plastic  zone  ahead  (<fr  =  0°)  of  the 
crack  as  a  function  of  applied  AK. 
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FIGURE  7-11.  Cyclic  plastic  zone  sizes  to  the  side  (  $  -  90°)  of 
small  cracks  as  a  function  of  applied  AK. 
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FIGURE  7-12.  Ratio  of  crack  size  to  plastic  zone  size  (a/rp)  for 
small  cracks  in  Astroloy  of  various  sizes. 
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FIGURE  V-  :  Comparison  of  crack  closure  behavior  of  small 

and  large  cracks  at  various  temperatures.  The 
dashed  lines  are  from  the  large  crack  data  of 
Soniak  ana  Remy  given  in  Section  4.4. 


The  above  results  In  the  Ni-base  alloys  are  qualitatively  similar 
to  those  observed  recently  in  a  7075  aluminum  alloy  [7.11]  as  shown  in 
Figure  7-14.  In  this  case,  both  the  small  and  large  crack  results  were 
obtained  using  stereoimaging.  These  high  resolution  measurements  clearly  show 
that  significantly  more  closure  occurs  in  large  cracks,  than  in  small  cracks, 
as  K|nax  is  decreased  and  a K  approaches  the  large  crack  threshold  for  fatigue 
crack  growth. 

Combining  Eqs  (7-6)  and  (7-7),  and  noting  the  identity  Pop/Pmax  = 
^op/^max*  <Pves: 

VWx  »  1  -  (1-R)  6  (7-8) 

Thus,  for  a  given  R,  Kop/Kmax  is  also  constant  for  small  cracks.  This 
relationship  is  in  agreement  with  recent  small  crack  data  of  Soniak  and 
Remy  [7.10]  on  necklace  Astroloy.  The  Soniak  and  Remy  data  for  R  =  0.1  are 
shown  in  Figure  7-15.  These  data  are  consistent  with  the  current  results  in 
that  crack  closure  in  small  cracks  was  found  to  be  insensitive  to  tempera¬ 
ture.  The  deviation  from  Eq.  7-8  which  occurs  above  a  Kmax  of  about  20  MPa/m 
is  believed  to  be  due  to  the  fact  that  the  small  cracks  have  grown  to  the 
length  where  they  now  exhibit  closure  more  typical  of  large  cracks.  Further¬ 
more,  the  fact  that  these  data  have  a  slope  R  is  consistent  with  the  linear 
relationship  between  J.'br'/K  1/Fmax  shown  for  large  cracks  in  Figures 
7-13  and  7-14. 

For  comparison.  Figure  7-13  also  contains  the  current  results  on 
Astroloy.  Although  both  sets  of  results  exhibit  the  same  functional  relation¬ 
ship,  the  current  results  give  a  higher  K0p/*max  than  do  those  of  Soniak  and 
Remy.  The  larger  amount  of  crack  closure  given  by  the  current  data  are  be¬ 
lieved  to  result  from  the  fact  that  these  measurements  were  acquired  local  to 
the  crack  tip,  while  those  of  Soniak  and  Remy  are  based  on  compliance  measure¬ 
ments  made  remote  from  the  crack  tip.  This  interpretation  is  consistent  with 
recent  measurements  of  Su  and  Sharp  [7.12],  obtained  using  a  laser-based 
interferometric  displacement  gage,  which  showed  that  Kop/Kmax  for  sma11 
fatigue  cracks  in  a  2024-T3  aluminum  alloy  increased  significantly  as  the 
mea:urement  location  was  moved  toward  the  crack  tip. 
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FIGURE  7-14.  Crack  closure  behavior  of  small  and  large  cracks  in 
a  7075  aluminum  alloy. 
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versus  K^gx  data  on  small  cracks 
he  current  study  with  those  of 
.10].  Note  that  the  latter  data 
k  closure  behavior  above  K^gx  = 


7.3  References 

7.1  0.  L.  Davidson  and  A.  Nagy,  "A  Low  Frequency  Cyclic  Loading  Stage  for 
the  SEM,"  J.  of  Physics,  Vol.  Ell,  1978,  pp.  207-210. 

7.2  A.  Nagy,  J.  B.  Campbell  and  D.  L.  Davidson,  Review  of  Scientific 

Instruments,  Vol.  55,  1984,  pp.  778-782. 

7.3  D.  R.  Williams,  D.  L.  Davidson  and  J.  Lankford,  "Fatigue  Crack  Tip 

Plastic  Strains  by  the  Stereoimaging  Technique,"  Exp.  Mech.,  Vol.  20, 
1980,  p.  134. 

7.4  D.  L.  Davidson,  D.  R.  Williams  and  J.  E.  Buckingham,  "Crack  Tip  Stresses 

as  Computed  from  Strains  Determined  by  Stereoimaging,"  Exp.  Mech., 

June  1983,  p.  242. 

7.5  P.  W.  Tan,  I.  S.  Raju  and  J.  C.  Newman,  Jr.,  "Boundary  Force  Method  for 
Analyzing  Two-Dimensional  Cracked  Bodies,"  NASA  Tech.  Memo.  87725,  NASA- 
Langley  Research  Center,  May  1986. 

7.6  0.  Lankford  and  D.  L.  Davidson,  "The  Role  of  Metallurgical  Factors  in 

Controlling  the  Growth  of  Small  Fatigue  Cracks,"  in  Small  Fatigue 
Cracks,  Eds.  R.  0.  Ritchie  and  J.  Lankford,  TMS/AIME,  1987,  pp.  51-71. 

7.7  0.  L.  Davidson,  Eng,  Fract.  Mech.,  Vol.  25,  1986,  pp.  123-132. 

7.8  D.  L.  Davidson,  Met.  Trans.,  Vol.  18A,  1987,  in  press. 

7.9  W.  F.  Brown,  Jr.  and  J.  E.  Srawley,  "Plane  Strain  Crack  Toughness 

Testing  of  High  Strength  Metallic  Materials,"  ASTM  STP  410,  1969. 

7.10  F.  Soniak  and  L.  Remy,  "Behavior  of  Long  and  Short  Fatigue  Cracks  in  a 

Power  Metallurgy  Superalloy  at  Room  and  at  High  Temperature,"  Fatigue 
' 87 ,  Eds.  R.  0.  Ritchie  and  E.  A.  Starke,  Jr.,  EMAS,  1987,  pp.  351-360. 

7.11.  D.  L.  Davidson,  "Small  and  Large  Cracks  ir.  Aluminum  Alloys,"  Acta  Met. , 
to  be  published. 

7.12  X.  Su  and  W.  N.  Sharpe,  Jr.,  "An  Experimental  Study  of  Closure  Behavior 
of  Naturally  Initiated  Small  Cracks,"  ibid,  pp.  343-350. 


8.  A  PHENOMENOLOGICAL  MODEL  FOR  THE  GROWTH  OF  SMALL  CRACKS 


This  section  uses  the  local  crack-tip  displacement  measurements  of  the 
previous  section  to  determine  the  mechanical  driving  force  for  the  growth  of 
small  cracks  based  on  the  relationship  between  crack-tip  opening  displacement 
and  aJ,  the  cyclic  J-integral.  To  facilitate  comparison  with  large  crack  data 
these  results  are  expressed  in  terms  of  an  equivalent  stress  intensity  factor, 
A«eq.  Since  these  results  are  based  on  local  crack-tip  measurements  they 
f  contain  the  influence  of  mixed-mode  loading,  plasticity  and  crack  closure. 

However,  since  thesr  same  factors  cannot  currently  be  computed  for  growing 
i  fatigue  cracks  subjected  to  cyclic  loading,  AKeq  cannot  simply  be  determined 

from  knowledge  of  crack  size,  component  geometry  and  remote  loading.  Thus,  a 
phenomenological  model  is  formulated  In  order  to  relate  AKgq  to  the  applied 
!  aK.  The  ability  of  AKgq  to  produce  crack  growth  kinetics  which  are 

independent  of  crack  size  is  also  assessed  by  comparing  crack  growth  rates  for 
both  small  and  large  cracks  as  a  function  of  AKeq.  The  engineering  utility  of 
this  approach  is  subsequently  examined  in  Section  9. 

8. I  The  Relationship  Between  Cyclic  CTOO  an  aJ 

It  is  of  interest  to  relate  the  cyclic  CTOD  values  measured  by 
stereoimaging  to  a  crack-tip  field  parameter  in  attempting  to  establish  a 
driving  force  for  the  growth  of  small  fatigue  cracks  under  elastic-plastic 
loading.  Available  elastic-plastic  solutions  for  stationary  cracks  under 
monotonic  loading  are  briefly  reviewed  since  they  provide  guidance  for  the 
adoption  of  an  operational  relationship  between  the  cyclic  CTOD  and  aJ. 

Hu  1 1  and  McCl intock  [ 8 . 1 1  have  provided  an  analytical  solution  for  a 
stationary  crack  under  antiplane  shear  in  an  elastic-perfectly  plastic 
material  under  small  scale  yielding.  The  relationship  between  Rice's  J- 
integral  [8.2]  and  CTOD  for  this  case  is  as  follows: 

J 1 1 1  = 


! 


"/4<t0  dt} 


(8-1) 


where:  is  the  J-integral  in  Mode  III,  dt  =  the  crack  tip  sliding 
displacement  under  Mode  III  loading,  and  tq  =  the  yield  stress  in  shear.  It 
is  interesting  to  note  that  the  constant  tt/4  in  Eq.  (8-1)  is  about  60  percent 
greater  than  the  value  of  1/2  obtained  by  Bllby,  Cottrell  and  Swinden  [8.31 
using  a  strip-yield  model  [8.4] . 

For  the  case  of  Mode  I  loading,  the  asymptotic  solutions  of  Hutchinson 
[8.5]  and  Rice  and  Rosengren  [8.6[  define  the  plastic  crack-tip  field  for  a 
power-hardening  material  experiencing  proportional  loading  according  to 
deformation  theory  of  plasticity.  One  difficulty  with  this  so-called  HRR 
field  is  that  the  crack  opening  displacement  tends  to  zero  as  the  crack-tip  is 
approached.  To  overcome  this  difficulty,  Tracy  [8.7],  following  Rise's 
suggestion,  proposed  that  the  CTOD,  dt,  be  operationally  defined  as  shown  in 
Figure  8-1.  This  definition,  combined  with  the  displacement  equation  for  the 
HRR  field,  gives: 


J  =  (l/dn)  0y$  dt  (8-2) 

where:  =  the  tensile  yield  stress;  and  dt  is  the  operationally  defined 
CTOO.  The  factor  dn  depends  primarily  on  stress  state  and  strain  hardening 
exponent,  and  mildly  on  the  materials  yield  strain,  as  determined  by 
Shih  [8.8]  using  finite  element  analysis.  For  the  case  of  plane  strain  and 
the  range  of  cyclic  strain  hardening  exponents  for  the  materials  and 
temperatures  of  interest  to  the  present  study  (n  =  0.145  to  0.029),  dn  is 
between  1.4  and  2.3.  For  the  case  of  a  perfectly  plastic  material  under  plane 
stress,  dn  =  1  and  Eq.(8-2)  gives  results  which  are  identical  to  those  of  the 
Ougdale-Barenblatt  model  (8.9,  8.10], 

The  primary  limitation  of  the  HRR  field  and  Eq.  (8-2),  with  respect  to 
fatigue  crack  growth,  is  the  fact  that  they  are  based  on  deformation  theory  of 
plasticity  which  assumes  that  the  plastic  strain  is  only  dependent  on  the 
current  state  of  stress  and  thus  is  independent  of  the  path  followed  on 
arriving  at  this  stress  state.  This  condition  is  generally  thought  to  be 
inapplicable  at  the  tip  of  a  growing  fatigue  cracks  where  unloading  occurs 
behind  the  advancing  crack  and  plastic  deformation  occurs  at  the  crack-tip. 
Nevertheless,  in  spite  of  these  analytical  limitations,  Dowling  [8.11,  8.12] 


FrGu«e  8_L 


has  demonstrated  the  utility  of  the  cyclic  J-integral,  aJ,  for  correlating  the 
growth  rates  of  large  cracks  for  a  variety  of  specimen  geometries  and  loading 
conditions  ranging  from  small  scale  yielding  to  fully  plastic. 

Lamda  [8.13|  has  proposed  that  the  Aj-integral  can  be  formally  defined 
by  replacing  the  :tress  tensor,  strain  tensor,  and  displacement  vector  in 
Rice's  J-integral  (8.2]  by  their  respective  ranges.  Furthermore,  he  has  shown 
that  the  resulting  AJ-integral  is  a  path-independent,  crack-tip  field 
parameter  for  elastic-plastic  materials,  provided  the  constitutive  behavior 
can  be  described  in  terms  of  the  cyclic  strain  energy  density  as  follows: 

°1j  =  dwij/d'ij  (8-3) 

where  and  -  •  j  are  the  stress  range  and  strain  range,  respectively,  and  W 
is  the  cyclic  strain  energy  density.  Eq.  (8-3)  requires  that  the  stress  range 
only  be  a  function  of  the  current  strain  range  and  not  the  deformation 
history.  As  discussed  in  Ref.  [8.14],  these  conditions  are  met  by  the  cyclic 
stress -strai n  curve  which  is  based  on  the  stress  and  strain  ranges  associated 
with  stable  stress-strain  hysteresis  loops.  Consequently,  materials 
exhibiting  stable  cyclic  plasticity  appear  to  be  treatable  as  a  nonlinear 
elastic  material  through  deformation  theory  of  plasticity. 

To  evaluate  the  applicability  of  the  above  concept,  Chan  [8.15,  8.161 
has  used  measured  stress-strain  fields  from  both  large  and  small  fatigue 
cracks  to  compute  values  of  Lamda' s  aJ-integral  for  different  line-contour 
integration  paths  around  the  crack  tip.  These  results  showed  that  A  J  was 
•-datively  independent  of  the  contour  path,  for  distances  greater  than  several 
C1CD  values  from  the  crack-tip.  Moreover,  the  average  contour-generated  aJ 
values  were  found  to  be  in  good  agreement  with  those  computed  directly  from 
C1CR  us^ng: 


•  J  •  A  3  ;  (8-4) 

■.here  a  is  the  crack-tip  stress  range  and  a  is  the  cyclic  CTOD.  Eq.(8-4)  ;s 
significant  in  that  it  provides  a  means  of  relating  the  measured  s  values 
to  ,\j  wnile  avoiding  much  of  the  uncertainty  associated  with  applying  the 


analytical  solutions  --  Eqs.  (8-1)  and  (8-2)  --  for  monotonically  loaded, 
stationary  cracks  to  growing  fatigue  cracks. 

8.2  Phenomenological  Hodel  for  the  Growth  of  Small  and  Large  Fatigue  Cracks 

Local  CTOD  values  for  small  cracks  in  CG  and  FG  Astroloy  at  several 
temperatures  were  determined  from  the  measured  displacement  fields  of 

Section  7.  These  values  were  then  used  to  compute  aJ  using  Eq.  (8-4).  To 

facilitate  comparison  of  the  resulting  small  crack  growth  kinetics  with  those 
of  large  cracks,  aJ  was  converted  to  an  equivalent  aK  using  the  following 
relationship  (8.13): 

a«eq  =  (EaJ)1/2  =  (E  Ao  6r)1/2  (8-5) 

where  the  crack  tip  stress  range  (ao)  was  determined  from  the  measured  crack- 
tip  strain  range  using  the  appropriate  cyclic  stress-strain  relation  from 
Section  3,  and  6r  is  the  vector-resolved  value  of  the  cyclic  CTOD  from  Mode  I 

and  II.  The  latte1'-  were  defined  as  the  respective  COD  valu'j  at  0.5  um  from 

the  crack-tip.  A  summary  of  the  crack-tip  deformation  parameters,  along  with 
the  resulting  aKeq  values,  is  provided  in  Table  8-1. 

Eq.  (8-5)  was  also  used  to  analyze  available  small  crack  data  in  a  ?075 
aluminum  alloy  at  room  temperature  [8.17].  Resulting  aKgq  values  for  both 
7075  aluminum  and  for  Astroloy  are  given  as  a  function  of  a<ec)  in  Fioures  8-7 
and  8-3,  respectively.  Results  in  both  materials  indicate  that  ,Keq 
related  to  the  applied  aK  as  follows: 

aKeq  =  Kp  +  saK  (8-6) 

The  first  and  second  terms  in  Eq.  8-6  are  viewed  as  the  plastic  and  e'astk 
contributions,  respectively,  to  the  crack  driving  force.  It  is  important  to 
note  that  since  aKeq  was  determined  from  measured  cyclic  crack- t-p 
displacements,  it  contains  the  influence  of  crack  closure,  as  wp'!  as 
plasticity  and  mixed  modeness.  Consequently,  the  second  (elastic)  term  in 
Eq.  (8-6)  is  merely  aKpff.  Thus,  for  the  case  of  Astroloy,  the  value  of  -  i1- 
ra a\  chmiM  mi  nr  i  dp  with  the  independent  crack  closure  measurements  m 


FIGURE  8-2.  Equivalent  versus  applied  IK  for  small  and  large 
cracks  in  a  7075  aluminum  alloy. 
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FIGURE  3-3.  Equivalent  versus  applied  AK  for  small  and  large 
cracks  in  Astroloy  at  various  temperatures. 


Section  7  which  showed:  =  0.35  =  e.  This  is  verified  by  the  fact 
that  the  straight  line  through  the  Astroloy  small  crack  data  in  Figure  8-3, 
constructed  using  9  =  0.35,  is  shown  to  provide  a  reasonable  fit  to  the  data. 

For  comparison.  Figures  8-2  and  8-3  also  contain  the  appropriate  A<eq 
values  for  large  cracks.  Since  large  crack  data  are  obtained  under  Mode  I 
loading  and  small  scale  yielding  conditions,  AKeq  for  large  cracks  is  simply 
given  by: 


AKeq  =  AKeff  =  AK  -  aKc1  (8-7) 

where  AKc1  is  the  portion  of  the  applied  AK  range  over  which  the  crack  remains 
closed.  For  the  case  of  the  Astroloy  data  in  Figure  8-3,  a  value  of  aK..,  = 

I 

4.0  MPa/m  has  been  used  based  on  the  large  crack  results  of  Soniak  and 
Remy  [8.18)  given  in  Section  4.2. 

Eqs.  (8-6)  and  (8-7)  define  a  phenomenological  model  of  the  mechanical 
driving  force  for  the  growth  of  small  and  large  cracks,  respectively.  This 
model  is  schematically  summarized  in  Figure  8-4,  together  with  the  corre¬ 
sponding  crack  closure  behavior,  to  emphasize  the  relationship  among  the 
various  parameters  of  the  model.  In  its  present  form,  the  model  contains 
three  fundamental  parameters:  K0  (or  aKci),  which  character izes  crack  closure 
in  large  cracks;  8,  which  characterizes  crack  closure  in  small  cracks;  and 
AKp,  which  characteri zes  the  plastic  contribution  to  the  crack  driving  force 
for  small  cracks.  As  indicated  in  Figure  8-4,  the  transition  between  so- 
called  "small"  and  "large"  crack  behavior  occurs  in  the  model  at  a  specific 
value  of  AK,  termed  AK^r,  which  ^s  defined  by  the  above  three  parameters  as 
fol lows: 


AKtr  =  (AKp  +  AKC] )/( 1  -  3)  (3-8) 

The  above  concept  would  appear  to  be  contrary  to  the  commonly  held  notion  that 
the  transition  from  "small"  to  "large"  crack  behavior  should  occur  at  a  given 
crack  size.  This  issue  is  further  discussed  in  Section  19. 
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FIGURE  8-4.  Phenomenological  model  for  the  growth  of  small  and 
large  fatigue  cracks:  a)  is  the  crack  closure  be- 
havior  and  b)  is  the  local  mechanical  driving  force 
versus  applied  AK. 


8.3  Small  and  Large  Crack  Growth  Rates  as  a  Function  of  aKeq  Versus  aKeff 

The  growth  rate  data  on  small  cracks  obtained  using  the  three-point  bend 
specimens  (Section  5)  and  the  single-edge-notched  tension  specimens 
(Section  7)  were  analyzed  using  the  above  phenomenological  model.  in  so 
doing,  it  is  implicitly  assumed  that  the  model  for  the  crack  driving  force  is 
independent  of  specimen  geometry.  This  is  believed  to  be  a  good  assumption 
since  the  model  is  based  on  direct  measurement  of  local  crack-tip  deformation 
and  crack  closure,  rather  than  on  quantities  inferred  from  remote 
measurements. 

Because  of  the  popularly  held  belief  that  crack  closure  alone  is  respon¬ 
sible  for  crack  size  effects,  crack  growth  rate  data  were  also  examined  in 
terms  of  AKeff.  Within  the  context  of  the  current  measurements,  this  is 
equivalent  to  setting  AKp  =  0  in  Eq.  (8-6)  of  the  above  phenomenological 
model . 

Figure  8-5  shows  the  room  temperature  data  on  CG  Astroloy  as  functions 
of  both  AKeq  and  aKeff.  The  large  crack  data  in  each  case  are  given 
parametrically  for  a  range  of  tKcj  values  intended  to  reflect  the  uncertainty 
in  the  determination  of  these  values.  Best  estimates  of  aKc]  based  on  the 
room  temperature  large  crack  closure  measurements  of  Soniak  and  Remy  [8.161 
for  necklace  Astroloy  are  3-5  MPa.m.  Based  on  these  results,  A<eq  appears  to 
provide  a  better  consolidation  of  large  and  small  crack  growth  kinetics  than 
does  A«eff. 

Using  the  same  format  as  that  above.  Figure  8-6  presents  results  on 
CG  Astroloy  at  600rC.  In  this  case,  the  best  estimate  of  aK  } ,  based  on  the 
65CTC  large  crack  data  of  Soniak  and  Remy  (8. 16),  is  2-4  MPa^m.  Once  again 
A«eq  appears  to  provide  a  better  consolidation  of  small  and  large  crack  data. 
However,  in  this  case  there  appears  to  be  an  downturn  in  tie  growth  nate 
kinetics  at  about  3  MPa/m.  Such  behavior  is  indicative  of  an  intrinsic 
threshold  for  fatigue  crack  growth,  “K  /th.  Similar  results  on  Waspaloy  at 
600  C  are  shown  in  Fiqure  8-7.  In  presenting  these  results,  it  is  implicitly 
assumed  that  the  above  empirical  modei  Darameters,  developed  for  small  cracks 
in  Astroloy,  are  also  applicable  to  Waspaloy.  This  assumption  seems 
reasonable  based  on  the  fact  that  the  growth  rate  behavior  of  small  cracks  is 
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FIGURE  8-5.  Comparison  of  crack  growth  kinetics  of  small 
(□)  and  large  (— )  cracks  in  CG  Astroloy  at 
25°C  in  terms  of  a)  equivalent  and  b)  effec¬ 
tive  stress  intensity  factors. 
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FIGURE  8-6.  Comparison  of  crack  growth  kinetics  of  small 
(□)  and  large  (  — )  cracks  in  CG  Astroloy  at 
600°C  in  terms  of  a)  equivalent  and  b)  effec¬ 
tive  stress  intensity  factors. 
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FIGURE  8-7.  Comparison  of  crack  growth  kinetics  in  small 
(□)  and  large  (— )  cracks  in  Waspaloy  at 
600°C  in  terms  of  a)  equivalent  and  b)  effec 
tive  stress  intensity  factors. 


MATERIAL  CONSTANTS  FOR  THE  PHENOMENOLOGICAL  MODEL  FOR  THE  GROkfTH  OF  SMALL  AND  LARGE  CRACKS 
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relatively  insensitive  to  alloy  composition  and  microstructure,  as  discussed 
in  Section  5.5.3. 

Based  on  all  of  the  above  results  it  appears  that  a  reasonable  represen¬ 
tation  of  the  crack  growth  kinetics  for  both  large  and  small  cracks  can  be 
obtained  from  the  following  equation: 


da/dN  *  C(aKeq  -  AKeq/th)n 


(8-9) 


where  C  and  n  are  empirical  constants,  and  Keq/th  the  intrinsic  threshold 
stress  intensity  factor.  Values  of  C,  n,  Kgq/,th  ?nd  AKcl,  b  and  AKp  |  from 
Eqs.  (8-6)  and  (8-7) j  for  the  materials  of  interest  here  are  given  in  Table  8- 

2. 
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9.  EFFECT  OF  SMALL  CRACK  BEHAVIOR  ON  COMPONENT  FATIGUE  LIFE 


The  fatigue  lives  of  actual  rig-tested  turbine  discs  have  been  computed 
using  the  phenomenological  model  for  small  crack  growth  presented  in 
Section  8.  Predicted  versus  measured  lives  are  compared  in  order  to  assess 
the  utility  of  the  proposed  model.  The  predicted  lives  based  on  the  small 
crack  model  are  also  compared  with  those  obtained  using  the  conventional 
linear  elastic  fracture  mechanics  and  large  crack  data.  The  influence  of  the 
large  crack  threshold  stress  intensity,  AKth,  on  the  latter  predictions  is 
also  examined. 

9.1  Analysis  of  Rig-Tested  Discs 

Engine  disc  test  data  were  taken  from  two  sets  of  rig  tests:  one  for 
Waspaloy,  and  a  second  for  Astroloy.  The  Waspaloy  discs  were  manufactured 
identically  to  the  Waspaloy  material  used  in  the  current  program.  The 
Astroloy  discs  were  manufactured  via  a  I160°C  hip-plus-forge  route  resulting 
in  a  necklace  microstructure  having  both  5  urn  and  40  urn  grain  sizes.  This 
microstructure  was  originally  chosen  to  provide  optimum  crack  growth  rate 
behavior  in  both  the  small  and  large  crack  regimes  (9.11.  In  contrast,  work 
in  the  current  program  has  focused  on  fine  and  coarse  grained  variants  of 
Astroloy  in  order  to  obtain  a  more  fundamental  understanding  of  the  mechanism 
of  small  crack  growth.  Nevertheless,  due  to  the  observed  insensitivity  of 
small  crack  growth  to  microstructural  variations,  the  coarse  grained  results 
from  Section  7  are  assumed  to  apply  to  the  analysis  of  the  necklace  Astroloy 
rig  tests.  However,  since  large  crack  growth  rates  have  been  shown,  both  in 
Section  4  and  elsewhere  19.1-9.31,  to  be  sensitive  to  microstructure,  the  data 
from  Section  4  were  used  in  making  the  large  crack  predictions.  Specifically, 
the  large  crack  data  of  Figure  4-14  and  Figure  4-21  were  used  for  the  large 
crack  predictions  in  necklace  Astroloy  and  in  Waspaloy,  respectively. 

Discs  in  the  rig  tests  were  typically  cycled  at  zero-maximum-zero  speed 
over  a  period  of  12  seconds,  thus  the  applied  load  ratio  was  near  zero  and  the 
nominal  cyclic  frequency  was  5  cpm. 

The  calculation  of  applied  ,\K  for  use  in  tqs.  (3-6)  and  (8-7)  were  based 
on  Pickard's  1 9 . 4 1  stress  intensity  factor  solutions  for  corner  cracks  ana 
surface  cracks,  and  incorporated  disc  stress  distributions  which  accounted  fm 


residual  stresses  imparted  during  manufacturing,  particularly  those  generated 
by  heat  treatment.  The  analysis  Techr.'nijes  used  have  been  previously 
validated  against  rig  tests  which  contained  preintroduced  (large)  cracks  of 
known  initial  size  [9.5)  Several  disc  features  were  analyzed  --  that  is 
bores,  di aphrams  and  boit  holes  --  depending  on  the  location  in  the  disc  where 
the  particular  failure  originated. 

Since  none  of  the  discs  analyzed  were  precracked,  the  i.v’tial  crack  size 
'•/as  assumed  to  be  either  a  semicircular  corner  or  surface  crack  having  the 
Radius  or  the  average  grain  size  or  the  material  --  that  is,  40  um  for  both 
Waspalcy  and  Astro ,  .•  'his  is  considered  a  physical  ly  realistic  lower  limit 
of  the  integration  since  crack  growth  across  the  initial  grain  was  rarely 
observed  in  the  sma 1 1  c-'ack  exper  ' '?nts.  In  one  instance,  failure  in  an 
Astro loy  disc  clear' v  originated  f - om  a  60  um  ceramic  inclusion,  located  near 
a  free  surface.  This  particu'ar  ect  was  idealized  as  a  90  um  surface  crack 

based  on  its  near- surf  ace  loca-.  ’or.  end  the  fact  that  little  or  no  initiation 
life  was  expected. 

With  regard  to  the  implementation  of  the  phenomenological  model  for  the 
growth  of  sma'i  cracks,  the  transition  between  small  and  large  crack  behavior 
was  made  at  a  transition  yK+;r  defined  by  Eq.  (8-8) 

9.2  Comparison  of  Predicted  Versus  Measured  Disc  Lives 

The  ratios  of  predated  to  measured  lives  for  all  of  the  discs  analyzed 
are  summarized  in  Figure  j-1.  The  life  ratios  labeled  "Small  Crack"  were 

obtained  using  small  crack  kiretics  oe low  aK^  and  large  crack  kinetics  above 

AKtr,  whereas  those  labeled  "large  f>ack"  were  obtained  using  the  large  crack 

kinetics  th^oughovi  the  computation.  As  one  might  expect,  the  small  crack 

predictions  g:,e  onsisterPy  smaller  life  ratios  than  do  the  large  crack 
predictions,  a  f tNr-jh  *  hr  r  * 1  cl; :  ,o  difference  between  the  two  sets  of  predic¬ 
tions  varies  great',,  depending  on  he  particular  case  analyzed. 

On  average,  the  life  ratio  for  the  small  crack  predictions  was  0.77. 
T  's  value  is  consistent  with  th»  'art.  that  crack  initiation  undoubtedly  con¬ 
tribute-;  to  the  measured  cyclic  lives,  but  wai  was  not  included  in  the 
predictions.  For  example,  tor  actual  discs  the  initiation  life  has  been 
estimated  to  he  between  5  percent  and  10  percent  of  the  total  fatigue  life 
id.fi  An  average  life  ratio  m  the  ase  of  the  large  crack  predictions  is 


FIGURE  9-1.  Influence  of  the  small  crack  effect  on  the  cyclic 
life  ratio  for  various  disc  temperatures  ar.u 
failure  locations. 


somewhat  more  difficult  to  interpret  due  to  the  large  overprediction  which 
occurred  for  the  case  of  the  Astroloy  diaphragm  cracks  at  200°C.  Excluding 
these  overpredictions  gives  a  life  ratio  for  the  large  crack  predictions  of 
1.0,  however,  including  the  overpredictions  gives  a  life  ratio  of  between  26 
and  infinity.  The  uncertainty  in  the  latter  value  arises  from  the  fact  that 
the  long  crack  predictions  in  the  case  of  the  Astroloy  diaphragm  crack  at 
200°C  predicted  a  life  which  was  several  orders  of  magnitude  larger  than  the 
measured  life  since  the  initial  aK  in  the  computation  was  less  than  the 
measured  aKth  value  for  necklace  Astroloy.  Overall,  the  small  crack 
predictions  appear  to  provide  a  more  realistic  view  of  the  behavior  of  the 
rig-tested  discs. 

The  actual  predicted  versus  measured  lives  are  compared  In  Figure  9-2. 
As  indicated,  most  of  the  predictions  are  within  a  factor  of  ±2  of  the 
measured  lives.  The  shortest  life  was  for  the  case  of  bolt  hole  cracking  in 
Astroloy  at  600°C  and  was  due  to  the  presence  of  a  near-surface  60  urn  ceramic 
inclusion.  The  longest  life  was  for  one  of  the  diaphram  cracks  in  an  Astroloy 
disc  at  200° C.  In  general,  these  results  show  that  the  deviation  between  the 
small  and  large  crack  predictions  increases  as  the  life  increases.  This  is 
not  surprising  since  the  lower  stresses,  and  associated  stress  intensity 
factors,  for  these  cases  cause  more  crack  growth  to  occur  below  aKtr  in  the 
small  crack  model.  The  major  deviations  are  therefore  associated  with  the 
marked  reduction  in  the  growth  rates  of  the  large  cracks  as  the  long  crack 
threshold  is  approached. 

In  order  to  further  examine  the  influence  of  the  large  crack  threshold 
on  the  predicted  lives,  selected  discs  were  analyzed  using  the  measured  long 
crack  growth  rates  above  about  10”^  mm/cycle  to  obtain  a  power  law  relation 
which  was  extrapolated  below  aK^.  In  other  words,  the  near-threshold 
behavior  was  ignored  and  a  simple  Paris  law  was  used  to  describe  the  growth 
rates.  These  results,  labeled  "Large/no  AKth",  are  given  in  Figure  9-3  where 
they  are  compared  with  the  previous  small  and  large  crack  predictions. 
Although  the  Large/no  aKth  predictions  still  overpredict  the  measured  lives, 
the  deviation  is  now  only  about  a  factor  o'  two.  Thus,  it  appears  that  a 
first  order  engineering  approximation  for  small  crack  effects  can  be  obtained 
by  simply  ignoring  the  long  crack  thresholu  and  using  a  simple  Paris  Law  to 
describe  the  fatigue  c^ack  growth  kineti;s.  The  utility  of  this  approach  as 
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FIGURE  9-2.  Comparison  of  measured  disc  lives  with  those  predicted 
using  small  versus  large  crack  growth  kinetics. 
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FIGURE  9-3.  Influence  of  various  crack  growth  rate  "1 
on  the  cyclic  life  ratio. 


wen  as  the  viability  of  the  phenomenological  small  crack  model  is  further 
discusses  further  in  Section  10. 
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10.  DISCUSSION 


Based  on  the  metallography,  selected  area  electron  channeling,  and 
fractography  of  Section  6,  it  appears  that  small  and  large  cracks  in  Ni-base 
superalloys  grow  by  a  similar  mechanism.  The  electron  channeling  patterns 
from  within  individual  grains  of  Astroloy  containing  small  cracks  showed  that 
the  preferred  crack  plane  was  [111}.  These  results  are  in  agreement  with 
those  of  Beevers  [10.1]  on  small  cracks  in  Astroloy  obtained  using  the  same 
technique.  Moreover,  this  preferred  cracking  orientation  is  consistent  with 
that  which  has  been  determined  for  large  cracks  using  a  variety  of  techniques 
[10.2-10.61.  These  results  are  also  consistent  with  the  slip  character  of  Ni- 
base  alloys,  as  well  as  the  crystallographic  facetting  observed  on  the  frac¬ 
ture  surfaces  of  both  small  and  large  cracks.  Furthermore,  the  fracture 
morphologies  for  small  and  large  cracks  also  exhibit  the  same  temperature  and 
applied  stress  dependence.  Thus,  it  appears  that  the  mechanism  of  crack  ex¬ 
tension  in  small  and  large  cracks  is  fundamentally  the  same  and  that  apparent 
differences  in  their  kinetics,  when  expressed  in  terms  of  a K,  are  due  to 
differences  in  mechanical  driving  force  and  not  to  actual  differences  in  their 
mechanism  of  extension. 

The  local  crack-tip  deformation  and  crack  closure  measurements  of 
Section  7  provide  a  direct  means  of  assessing  those  processes  which  influence 
the  crack  driving  force  for  small  cracks.  Based  on  the  ratio  of  crack  size  to 
plastic  zone  size  (typically,  a/rp  =  1  to  5),  it  is  clear  that  the  linear 
elastic  fracture  mechanics  (LEFM)  concept  of  small  scale  yielding  (SSY:  a/rp 
>50)  [10.7]  does  not  hold  for  small  cracks.  As  evidenced  by  the  data  obtained 
on  CG  Astroloy  at  room  temperature,  SSY  can  even  be  violated  when  the  applied 
stress  is  less  than  the  materials  yield  strength.  Violation  of  the  SSY  con¬ 
dition  is  consistent  with  enhanced  crack-tip  strains  and  altered  strain  dis¬ 
tribution  exhibited  by  small  cracks,  relative  to  those  which  exist  at  large 
cracks  at  the  same  aK.  In  view  of  this  situation,  {t  is  not  surprising  that 
the  crack  closure  behavior  of  small  and  la^ge  cracks  differs  so  significantly, 
as  shown  in  Section  7.4.  In  addition  to  these  effects,  small  cracks  in  coarse 
grained  Ni-base  superalloys  also  exhibit  significant  mixed-Mode  displacements 
due  to  their  preferred  cracking  planes  within  the  crystallite.  This  effect 
also  needs  to  be  included  in  the  crack  driving  force  for  small  cracks. 


As  demonstrated  in  Sections  7  and  8,  local  cyclic  CTOD  measurements 
appear  to  capture  all  of  the  above  effects.  Thus,  converting  these  local  CTOD 
measurements  to  aJ,  or  equivalently  to  AK0q,  provides  a  measure  of  the  crack 
driving  force  which  results  in  crack  growth  kinetics  which  are  independent  of 
crack  size,  within  experimental  scatter  (Section  8.3).  These  results  also 
provide  insight  into  why  previous  attempts  to  explain  small  crack  effects 
solely  on  the  basis  of  crack  closure  have  met  with  limited  success  [10.8- 
10.11].  In  these  instances,  it  was  generally  assumed  that  the  growth  of  small 
cracks  was  a  closure-free  process.  This  view  probably  developed  due  to  the 
fact  that  remote  compliance  measurements  were  of  insufficient  precision  to 
detect  crack  closure  in  small  cracks.  Consequently,  the  primary  difference 
between  small  and  large  cracks  was  attributed  to  the  closure  which  occurred  in 
large  cracks,  and  was  measurable  by  remote  techniques.  In  view  of  the  current 
results,  it  is  clear  that  the  erroneously  neglecting  crack  closure  in  small 
cracks  is  approximately  offset  by  erroneously  neglecting  the  plasticity  con¬ 
tribution  to  the  crack  driving  force.  When  crack  closure  is  more  precisely 
measured  --  as  in  the  present  study,  or  in  a  parallel  study  (10.121  --  it  is 
found  that  AKefrf  only  partially  accounts  for  the  observed  crack  size  effects. 

As  indicated  in  Section  8.2,  the  current  version  of  the  phenomenological 
model  for  the  growth  of  small  cracks  suggests  that  the  transition  from  "small" 
to  "large"  crack  behavior  occurs  at  a  specific  value  of  a K,  termed  AKtr.  This 
concept  is  inconsistent  with  the  popular  view  that  this  transition  should 
occur  at  a  given  crack  length,  atr.  However,  in  spite  of  numerous  attempts  to 
precisely  define  atr,  this  quantity  remains  elusive  --  for  example,  see 
Refs.  (10.13-10.151.  This  situation  suggests  that  the  transition  from  "small" 
to  "large  "crack  behavior  may  depend  on  additional  factors.  For  example,  in 
light  of  the  current  results,  it  is  conceivable  that  atr  could  also  depend  on 
the  extent  of  the  plasticity  contribution  to  the  crack  driving  force  (sKp),  as 
well  as  the  relative  closure  character istics  of  small  and  large  cracks 
(b  and  aKc1 ,  respectively)  --  all  of  which  can  differ  significantly  with 
material  and  applied  stress  level. 

With  regard  to  the  influence  of  applied  stress  level,  it  should  also  be 
pointed  out  that  one  limitation  of  the  current  version  of  the  proposed 
phenomenological  model  is  that  AKp  is  taken  to  be  constant  at  about  3  MPa-m. 
Intuitively,  it  seems  that  AKp  should  be  a  function  of  the  applied  stress 


range  (aS),  and  perhaps  the  load  ratio  (R).  Unfortunately,  we  currently  have 
insufficient  data  to  define  this  functional  relationship.  Assuming  AKp  = 
Kp(AS,  R)  then  it  is  also  likely  that  AKtr  =  AKtr(AS,  R).  In  order  to  verify 
this  additional  experiments,  conducted  at  higher  stresses  and  under  strain 
controlled  conditions,  would  be  required. 

In  spite  of  the  above  limitations,  the  proposed  phenomenological  model 
is  capable  of  providing  more  realistic  predictions  of  component  life  than  the 
classical  LEFM  approach  in  which  crack  growth  kinetics  are  defined  in  terms  of 
da/dN(A<)  and  AKth.  These  results  raise  important  questions  regarding  the 
usefulness  of  large  crack  threshold  data  for  fatigue  life  assessment  in 
situations  where  cracks  initiate  and  grow  from  a  free  surface. 

The  utility  of  AKth  can  also  be  questioned  based  on  the  fact  that  it  can 
be  extremely  i rreproduci ble  under  certain  conditions.  This  is  particularly 
true  for  Ni-base  superalloys  at  room  temperature  where  AKth  values  can  range 
from  about  7  to  20  MPa/m  with  significant  differences  also  occurring  in  crack 
growth  rates  to  aK  values  appreciably  above  AKth  [10.16-10.18].  Although  the 
data  generated  in  the  current  study  suggest  that  this  variation  in  AKth  i$ 
sensitive  to  grain  size,  other  studies  have  reported  fatigue  crack  growth 
rates  in  Astroloy  to  be  insensitive  to  grain  size  [10.16].  Interestingly,  the 
current  study  found  that  this  irreproducibi 1 ity  disappears  at  elevated  temper¬ 
ature  (Section  4),  thereby  suggesting  that  it  may  be  related  to  the  extremely 
rough  fracture  surfaces  and  associated  roughness- i nduced  crack  closure  which 
develop  at  room  temperature.  This  view  is  consistent  with  reported  dif¬ 
ficulties  in  the  use  of  the  compliance  technique  to  monitor  fatigue  crack 
growth  in  Al-Li  alloys  which  also  develop  extensive  roughness-induced  crack 
closure  [10.19-10.20).  Such  difficulties  may  contribute  to  the  recently  re¬ 
ported  geometry  dependence  of  near-threshold  fatigue  crack  growth  rates  in  Ni- 
base  superalloys  ( 10.21 1,  since  remote  crack  closure  measurement  are  known  to 
be  geometry  dependent  [ 10.221 . 

Roughness-induced  crack  closure  is  also  largely  responsible  for 
microstructural  effects  in  the  near-threshold  fatigue  crack  growth  rates  of 
large  cracks  [  10.23-10.25 1 .  In  the  case  of  small  cracks  where  crack  closure 
is  reduced  due  to  the  larger  crack  openings,  the  fatigue  crack  growth  rates 
become  relatively  insensitive  to  microstructure,  as  discussed  in 
Section  5.6.3.  Thus,  over-emphasis  of  the  near-threshold  large  crack  data  can 


result  in  improper  alloy  selection.  For  example,  this  appears  to  have  Deen 
the  case  in  Astroloy  where  the  necklace  microstructure,  having  duplex  grain 
sizes  of  5  um  and  40  urn,  is  commonly  selected  over  a  finer  grained  variant 
based  on  a  comparison  of  their  large  crack  growth  rates  in  the  near-threshold 
regime.  In  retrospect,  the  fine  grained  alloy  would  appear  to  have  been  the 
better  selection  based  on  the  insensitivity  of  small  crack  growth  rates  to 
microstructure,  combined  with  the  superior  strength  and  crack  initiation 
resistance  of  fine  grained  alloys. 


The  results  of  the  current  study  serve  to  establish  the  utility  of  AKeq 
as  a  driving  force  for  the  growth  of  both  small  and  large  fatigue  cracks. 
Ideally,  one  would  like  to  compute  the  local  value  of  AKeq  using  the  applied 
\K  as  prescribed  boundary  conditions.  In  theory,  it  is  possible  to  model  the 


mixed-mode,  elastic-plastic  crack  growth  process  under  cyclic  loading  using 
numerical  analysis.  In  fact,  such  techniques  have  been  used  to  model  aspects 
of  the  problem.  For  example,  the  finite  element  technique  has  been  applied  to 
the  analysis  of  anisotropic  effects  in  surface  cracks  (see  the  Appendix)  and 
crack  closure  in  large  Fatigue  cracks  [ 10.26  10.291,  while  the  boundary 
integral  equation  technique  has  been  applied  to  the  analysis  of  cracks  in  a 
previously  existing  plastic  field  i 10.30).  However,  due  to  the  complexity  and 
corresponding  cost  of  such  analyses  only  selected  cases  can  be  analyzed. 
Consequently,  these  results  are  essentially  numerical  "experiments"  and  thus 
are  difficult  to  ueneralize. 


Simpler  analyses,  such  as  those  based  on  the  Dugdale-Barenblatt 
idealization,  provide  a  viable  alternative  to  numerical  analysis.  These  have 
served  to  enhance  our  understanding  of  fatigue  crack  closure  in  large  versus 
small  cracks  (10.31,10.321,  and  interactions  of  small  crack  with  the  material 
microstructure  (10.33,10.34).  One  of  the  latter  cases  (  10.34)  is  particularly 
relevant  to  the  current  results  in  that  it  elucidates  the  contributions  to  oKp 
from  variations  in  the  yield  strength  of  individual  grains  due  to  their 
crystal lographic  orientation.  However,  one  of  the  practical  limitations  of 
Dugdale-Barenblatt  type  models  is  the  fact  that  the  solutions  become  unbounded 
when  the  applied  stress  equals  the  materials  yield  strength.  Unfortunately, 
this  is  the  region  of  greatest  practical  interest  --  both  for  analysis  of 
elevated  temperature  specimens  and  for  application  of  results  to  life 
prediction  in  turbine  discs. 
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11.  CONCLUSIONS 


Small  cracks  in  all  three  alloys  were  found  to  initiate  within  slip 
bands,  however  these  were  sometimes  found  to  emanate  from  pores, 
particularly  in  FG  Astroloy.  Selected  area  electron  channeling  patterns 
on  grains  containing  room  temperature  cracks  in  CG  Astroloy  showed  that 
the  formation  of  small  cracks  occur  preferentially  on  {111}  planes,  as 
might  be  expected  based  on  the  slip  character  of  this  material. 

The  classical  small  crack  problem  occurs  at  elevated  temperature  (200°C 
to  600°C),  as  well  as  room  temperature,  in  Ni-base  superalloys,  thus  the 
phenomenon  can  be  of  practical  significance. 

Based  on  the  similarity  in  the  fracture  surface  morphology  of  small  and 
long  cracks  as  a  function  of  temperature  and  applied  stress,  it  is 
concluded  that  the  mechanism  of  crack  propagation  is  similar  in  both 
cases,  thus  crack  size  effects  are  attributed  to  differences  'n 
mechanical  driving  force. 

The  crystallographic  nature  of  crack  growth  in  small  cracks  in  Ni-base 
superalloys  leads  to  mixed-mode  displacement  of  the  crack  faces.  The 
degree  of  mi xed-modeness  depends  on  the  specific  orientation  of  the 
crack  relative  to  the  loading  axis  and  on  the  magnitude  of  aK. 

The  kinetics  of  small  crack  growth  at  both  room  temperature  and  600°C 
are  not  measurably  different  in  Astroloy  and  Waspaloy  having  similar 
grain  sizes  (30  um  versus  42  urn,  respectively)  indicating  an  insen¬ 
sitivity  to  alloy  composition  and  microstructural  details.  In  contrast, 
limited  data  at  600"C  in  FG  Astroloy  did  not  exhibit  the  classical  small 
crack  effect  and  were  found  to  grow  significantly  slower  than  small 
cracks  in  CG  Astroloy. 

Strain  fields  at  the  tips  of  small  cracks  in  all  three  alloys  and  over  a 
range  of  temperatures  was  surprising  homogeneous,  relative  to  the  slip 
character  of  these  alloys,  indicating  that  multiple  slip  occurs  in  the 
highly  strained  crack-tip  region.  Consequently,  the  measured  strain 
fields  provide  no  clue  as  to  why  crack  growth  occurs  on  a  particular 
crystal lographic  plane. 

The  ratio  of  crack  size  to  measured  crack-tip  plastic  zone  size  (a/r  ) 
was  found  to  range  from  I  to  5,  consequently  the  linear  elastic  fracture 
mechanics  concept  of  small  scale  yielding  (a/r  >  50)  did  not  hold  for 
most  of  the  small  cracks  examined. 


Ill 


8.  Although  the  crack  opening  displacements  of  both  small  and  large  fatigue 
cracks  exhibit  the  same  /r-dependence  as  predicted  from  elastic  theory, 
meaningful  aK  values  cannot  be  determined  by  combining  these 
measurements  with  the  elastic  displacement  equations  for  a  static  crack. 
This  dichotomy  is  believed  to  be  caused  by  cyclic  plasticity  at  the  tip 
of  a  fatigue  crack  which  generates  a  local  stress  field  which  differs 
fundamentally  from  that  found  by  either  elastic  or  plastic  analysis  of 
monotonical ly-loaded  stationary  cracks. 

9.  Significant  crack  closure  occurs  in  both  small  and  large  cracks, 
aUhough  its  magnitude  and  dependence  on  fatigue  loading  variables 
differs  in  each  case:  [n  long  cracks  aKc1  =  KQ(1  -  R),  whereas  in  small 
cracks  aKc^  =  aAK,  where  «c  and  a  are  material  constants. 

10.  Formulating  the  driving  force  for  the  growth  of  fatigue  cracks  in  terms 
of  AKeq  which  accounts  for  crack  closure,  plasticity  and  mixed-modeness 
in  small  cracks,  and  crack  closure  in  long  cracks,  results  in  crack 
growth  kinetics  which  are  independent  of  crack  size. 

11.  Using  a  phenomenological  model  for  the  growth  of  small  cracks  based  on 
the  above  results  gave  better  predictions  of  disc  lives  compared  to 
conventional  linear  elastic  fracture  mechanics  and  long  crack  data.  The 
latter  approach  can  significantly  overpredict  disc  life  due  to  the 
misleadingly  high  A<th  values  exhibited  by  large  cracks,  compared  to 
those  in  small  cracks. 

12.  The  near  threshold  crack  growth  rates  exhibited  by  large  cracks  are  in 
general  not  applicable  to  life  prediction  of  actual  components  and 
furthermore  can  lead  to  improper  material  selection. 
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APPENDIX 


FINITE  ELEMENT  ANALYSES  IN  NON-ISOTROPIC  MATERIALS 

Linear  elastic  fracture  mechanics  has  been  found  to  be  applicable  to 
the  prediction  of  fatigue  crack  propagation  behavior  in  most  gas  turbine  disc 
materials  where  the  crack  in  question  is  long  compared  to  material  micro- 
structural  dimensions.  An  investigation  has  been  undertaken  to  identify 
whether  linear  elastic  fracture  mechanics  can  be  used  to  correlate  short  crack 
growth  behavior  when  account  is  taken  in  the  fracture  mechanics  solutions  of 
material  anisotropy  and  inhomogeneity.  This  has  been  achieved  by  constructing 
20  and  3D  finite  element  models  of  cracks  in  anisotropic  bodies  and  of  cracks 
in  local  regions  of  anisotropy  within  an  otherwise  isotropic  body.  The  latter 
case  represents  a  crack  growing  within  an  'individual  microstructural  unit. 

This  appendix  describes  the  method  developed  for  the  calculation  of 
crack  tip  stress  intensity  factors  for  anisotropic  materials.  The  first  part 
of  this  appendix  covers  work  on  2-dimensional  models  while  the  remainder 
covers  the  subsequent  work  which  developed  the  method  for  3-dimensional  models 
with  local  anisotropy. 

A.l  Two-Dimensional  Finite  Element  Models 

A  two-dimensional  model  was  constructed  using  8  node  isoparametric 
quadrilateral  plane  strain  elements.  Figure  A-l.  The  principal  model 
consisted  of  a  2  inch  (50.8  m)  perpendicular ,  edge  crack  contained  within  a 
4  inch  (101.6  mm)  square  region  of  anisotropic  material.  The  remainder  of  the 
52  inch  (1320.8  mm)  x  20  inch  (508  mm)  specimen  was  assumed  to  be  isotropic. 
Since  the  analyses  were  linear  elastic,  the  scale  of  the  model  was 
unimportant.  Additional  results  were  obtained  for  the  cases  where  the  whole 
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of  the  model  was  anisotropic  or  isotropic.  The  relative  dimensions  of  the 
model  were  chosen  to  allow  comparison  of  the  calculated  stress  intensities 
with  published  results  for  the  isotropic  case  { A. 1 ] .  The  crack  tip  elements 
were  isoparametric  quadri laterals  degenerated  to  3  sides  and  6  nodes  with  the 
mid-side  nodes  closest  to  the  crack  tip  moved  to  Ine  quarter  point  position  to 
model  the  elastic  singularity.  The  detail  close  to  the  crack  is  shown  in 
Figure  A-2.  The  model  was  elastically  loaded  by  a  fixed  displacement  of  the 
top  and  bottom  edges  thereby  preventing  rotation. 

The  anisotropic  region  was  modelled  with  the  elastic  properties  of 
SRR99,  a  single  crystal  nickel  base  turbine  blade  material.  The  isotropic 
material  was  represented  by  MAR-M-002,  an  equiaxed  nickel  base  turbine  blade 
alloy  of  a  similar  composition  to  SRR99.  These  materials  were  chosen  because 
of  the  availability  of  elastic  property  data.  It  is  expected  that  the  ratio 
of  anisotropic  to  isotropic  elastic  properties  for  a  disc  alloy  such  as 
Waspaloy  will  be  similar  to  those  considered  for  SRR99  and  MAR-M-002.  The 
elastic  properties  at  20’C  used  in  the  models  are  given  in  Table  A-l. 

The  finite  element  models  were  run  using  the  MARC  program,  level  K1.3. 
A.I.l  Elastic  Properties  of  Rotated  Anisotropic  Material 

SRR99  is  a  face  centered  cubic  nickel  base  material.  Considering 
Hooke's  Law  in  a  generalized  form. 

■  &  cu  ej  <#-1> 
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TABLE  A- 1 


ELASTIC  MATERIAL 


ANISOTROPIC ^MATERIAL 
(SR99  at  20~C) 

ELASTIC  STIFFNESS  GPa 


C 

c 

c 


11 

12 

66 


=250 
=156 
=  115 


ISOTROPIC  MATERIAL 
( MAR-M-002  AT  20  C) 

YOUNG ' S  MODULUS 
POISSON'S  RATIO 


PROPERTIES 


;=227.5  GPa 
=0.3 
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Provided  the  axes  of  elastic  and  global  material  symmetry  are  coinci¬ 
dent,  the  elastic  stiffness  matrix  can  be  represented  by  three  independent 


values, 

C11  =  C22  =  C33 
C12  =  C13  =  C23 
C44  =  C55  =  C66 


All  other  =  0 

Typical  values  of  for  SRR99  at  20°C  are  given  in  Table  A-l. 

In  the  case  where  the  global  and  elastic  material  axes  are  not  coinci¬ 
dent,  the  form  of  Hooke's  law  for  a  material  such  as  SRR99  which  exhibits  axes 
of  symmetry  of  order  4  reduces  to  a  stiffness  matrix  with  7  independent  con¬ 
stants. 
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The  full  6x6  3-dimensicnal  form  of  the  stress-strain  ieiationshio  has 
been  considered  to  allow  the  computer  code  written  for  rotation  to  be  used  for 


3-dimensional  models. 

To  obtain  the  global  elastic  properties  of  a  region  of  anisotropic 
material  with  known  local -properties,  the  following  transformation  was 
performed.  The  conventions  observed  are  given  in  Figure  A-3. 

Within  the  local  system. 

E  =  (  S  1 0 

Within  the  global  system 
e'  =  [S' | o 1 

and  the  elastic  compliance  matrices  (SI  and  [S' I  are  related  as  follows 

[S' |  =  N l S 1 NT 

Where  [N]  is  given  by 
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where  a.  ,  and  -y '  are  the  direction  cosines  of  the  angles  relating 
the  global  system  to  the  local  system  as  follows 


For  the  2-dimensional  model  considered,  the  matrix  of  direction  cosines 
i s  given  as  fol lows. 


C 

-S 

0 


0 
0 
1 J 


where  [N|  is  given  by 


S 

0 

0 

0 

-2SC 


C* 

0 

0 

0 

2SC 


0 

0 

1 

0 

0 

0 


C  =  Cose 
S  =  Sine 


0 

0 

0 

c 

s 

0 


0 

0 

0 

-s 

c 

0 


sc 

-sc 

0 

0 

0 

1 


The  variation  of  the  elastic  stiffness  For  SRR99  versus  rotation  of  the 
property  axes  about  the  common  z  axis  is  shown  in  Figure  A-4. 


A. 1.2  Calculation  of  Elastic  Strain  Energy  Release  Rate  1 G 1  Using  MARC 

The  elastic  singularity  at  the  crack  tip  was  modelled  using  singular 
isoparametric  elements  with  quarter  point  mid-side  nodes  (A.2).  MARC  has  the 
facility  to  determine  the  strain  energy  difference  by  means  of  a  differential 
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156  GPu 
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"DEGPt  r 

Stiffness  Rotation 


f 


stiffness  calculation  for  a  pre-defined  crack  advance  [A. 3).  The  finite 
element  model  with  the  perpendicular  edge  crack  was  run  with  rotated  local 
anisotropy,  rotated  whole  model  anisotropy  and  complete  isotropy.  The  values 
of  ' G '  for  the  col  inear  differential  crack  advance  and  the  remote  stress 
derived  from  the  MARC  program  output  are  given  in  Table  A-2. 


A. 1.3  Plane  Elastic  Anisotropy 

Plane  problems  in  terms  of  Hooke's  Law  are  reduced  to  a  form  involving 
a  maximum  of  six  independent  elastic  constants. 
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L  XY  J 

where  S^j  refers  to  either  the  plane  stress  compliance  a^-  or  the  plane  strain 
compliance  b^.  The  relationship  between  a ^ j  and  b ^ j  is  given  in  Table  A-3. 

The  method  employed  for  the  consideration  of  anisotropic  elasticity  was 
outlined  by  Sih  et  al  [A. 4].  The  analytic  functions  required  for  the 
expressions  describing  stresses  and  displacements  require  parameters  which  are 
the  conjugate  pairs  of  roots  of  the  following  characteristic  equation. 

"  ^16u^  +  (^12  +  ^66^  _  ^26u  +  ^22  =  ®  (A~3) 

The  general  expressions  for  stress  and  displacements  in  a  small  region 
around  the  crack  tip  for  the  symmetric  and  skew  symmetric  loadings  are  as 
follows,  where  r  is  a  radial  distance  from  the  crack  tip  rotated  at  an 
angle  e  anticlockwise  from  the  colinear  crack  direction. 
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ELASTIC  STRAIN  ENERGY  RELEASE  RATE  1 G ‘  AND  REMOTE  STRESS 
FOR  A  STRAIN  OF  0.3846E-3 
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TABLE  A- 3 


RELATIONSHIP  BETWEEN  PLANE  STRAIN  (bij)  AND 
PLANE  STRESS  ( a i j )  ELASTIC  COMPLIANCE 


Plane  strain:  b. 


Plane  stress:  a. 


1  -5  “ 


<1-1 


a33 

fi12 

=  “21  = 

a33 

2 

a22a33  *  a23 

—  h  - 

a16a33  "  a13a36 

a33 

“16 

“61 

a33 

a66a33  “  a362 

K 

=  b  _  = 

a26a33  '  a23a36 

“26 

"  “62 

A 
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°x  = 


/ ( 2nr ) 


0  =  Re 
y  J ( 2nr ) 


A 

uru2 

( - - - 

J ( cose+y2sin9 ) 

-L 

/ ( cosS+y^sine ) 

1 

(  1,1 
/(cos9+y2sin9 ) 

y2 

uru2 

/ (cos9+y1sin8 ) 

yly2 

1 

( - 

v/,(cos6+y1sin9 ) 

1 

ui-u2 

/ ( cos6+y2sin8 ) 

Displacements  parallel  and  perpendicular  to  crack  face. 


K1/(2r)  1 

u  =  — _ — . —  Re  - (VJ1p2Acose+u2sin6] 


Kx/(2r)  1 

v  =  -  Re  - ( y.q  /[cos9+y.,sine  ] 

/(tt)  yl'^2  12  2 


u2P1>/[cose+u1sine  ] 


(A-4) 


y2ql/^cose+ylsine  ^ 


°x  = 


ulw2 


/(2irr)  ^ ( cos9+y2sin9 ) 


/ (cose+y^sine ) 


°Y  =  T~~— '  Re 
Y  / (  2tir ) 


ll~u2  ^  (cose+y2sin0 )  /(cosS+y^ine ) 


T  =  - —  Re 

^  /(2nr) 


-  (_ - i - £ - ) 

pl~u2  *Mcose+u1sine)  /  ( cos9+y2sin9 ) 
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Displacements  parallel  and  perpendicular  to  crack  face. 


K2/(2r)  1 

u  =  -  Re  -  ( cos 6+U-, sine  ]  -  p.. /[cose+y,  sine  ] ) 

/(it)  ’■  c  1  i 


yl~y2 


(A-5) 


K2/(2r)  1 

v  =  -  Re  -  (q^/tcose+u^sine]  -  q, /[cose+y. sine ] 

/(it)  *  1 


yl"y2 


Pj  "  S11  yj  +  S12  "  S16yj 
qj  =  S12  uj  +  (S22/uj  ]  ~  S26 

For  mixed  mode  conditions,  the  strains  energy  release  rate  for  colinear 
crack  extension  is  given  by 


S22  KI  yl+y2  S11  KII 

6  =  — — -  Re  i  — — -  +  — — —  Im 

2  UX  U2  2 


|W 


(A-6) 


+  -  Kx  Kjj  S22  Re.[i/(P1U2)  )  +  S11Im(y1u2) 


For  the  orthotropic  case  when  the  crack  lies  in  a  plane  of  symmetry,  both 
Re(i/uj^2)  and  Im(^u2)  equal  zero  and  the  modes  of  cracking  are  independent 
and  can  be  written  as  follows. 


G  =  K. 


2  S11S22 


1 1  r 


(C  /C  J*  2S12  +  b66 
b22/Sll '  +  ~ - 


2  i  2S,  -  +  Sc<- 

+  k  c  _i±  /c  e  \  I  j.  12  66 

II  /,  S22S11  *  - - - 
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The  values  of  and  the  relevant  functions  of  and  u2  for  the  elastic 

conditions  considered  are  given  in  Table  A-4. 

Because  of  difficulties  encountered  in  solving  simu1 taneously  for  Kj 
and  K j  j  by  the  method  of  differential  crack  advance,  although  this  has  been 
done  for  isotropic  materials  [A. 5),  attention  was  given  to  the  method  of 
displacement  extrapolation, 

A. 1.4  Displacement  Extrapolation  for  Calculation  of  Kj  and  Kjj 

For  a  small  region  surrounding  the  crack  tip,  the  expressions  for 
displacements  parallel  and  perpendicular  to  the  crack  face  have  been  given, 
Eqns.  (A-4,  A-5).  Putting  e  =  180  degrees  such  that  V'  lies  along  the  crack 
face  modifies  those  expressions  as  follows. 

Symmetric  Loading 


u  =  Kj/t  2r/m )  Re 


,>‘lP2 


v  =  Kj/^r/iT)  Re 


(tJlq2 


U2P1 1 


^l1 
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Skew  Symmetric  Loading 


u 


Kn/(  2x/v )  Re 


pl> 


v  =  Kjj/ (  2r/ir )  Re 


Pj  =  S11  +  S12 


W1  "  y2 


- 


S16  uj 


qj  “  S12  uj  +  S22  ‘  S26 


( A- 8 ) 
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SOLUTIONS  TO  CHARACTERISTIC  EQUATION 
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Summing  the  corresponding  displacement  above  for  the  symmetric  and  skew 
symmetric  cases  gives  the  displacements  for  the  mixed  mode  condition. 
Substitution  of  the  isotropic  properties  into  those  expressions  leads  to  the 
foi lowing. 

2 

u  =  K^j/(2r/ir)  (  2  S  ^  =  1-  v  for  plane  strain 

E 

2 

v  =  KI/(2r/ir)  (2  ;  ^2  =  1-  v  for  plane  strain 

E 

Calculation  of  the  relative  movement  of  opposite  pairs  of  crack  face  nodes 
along  and  normal  to  the  crack  plane  and  their  substitution  into  the 
displacement  Equations  (A-7.A-8)  will  lead  to  values  of  Kj  and  Kn  at 
distances  ' r '  from  the  crack  tip.  Extrapolation  of  these  values  to  r  =  0 
gives  the  crack  tip  Kj  and  Kjj. 

To  allow  a  comparison  between  the  calculated  K  values  for  isotropic 
bodies  and  those  published  [A.l],  K  has  been  normalized  with  respect  to  Ko  = 
0/ (na) 

o  =  remote  stress 
a  =  crack  length 

The  values  calculated  for  K  and  K/Ko  are  given  in  Table  A-5.  The 
result  for  the  isotropic  case  is  within  the  bounds  expressed  by  Rooke  and 
Cartwright  [A.l],  i.e.,  ,  1.180+0.047  max. 

A. 1.5  Crack  Extension  and  Branching 

The  extension  of  a  crack  within  a  single  crystal  material  cannot  be 
assumed  to  lie  perpendicular  to  the  remote  stress  direction  or  col  inear  with 
the  original  defect.  Putting  aside  the  crystallographic  nature  of  mixed  mode 
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single  crystal  cracking,  consideration  has  been  given  to  the  hypothesis  that 
crack  propagation  will  occur  In  a  radial  direction  In  which  the  elastic  strain 
energy  release  rate  is  a  maximum.  The  value  of  ' G '  for  a  radial  crack  advance 
has  been  calculated  by  two  methods:  (a)  differential  crack  advance  within  the 
finite  element  program,  and  (b)  use  of  nodal  displacements  from  the  finite 
element  program  output  to  supply  stress  intensity  factor:  end  hence  'C  for 
crack  branching. 

The  differential  crack  advance  technique  was  used  for  the  models  with 
local  at. isotropy  rotated  by  30°  and  0°.  The  nodal  advance  used  In  this 
technique  Is  shown  in  Figure  A-5. 

The  method  using  nodal  displacements  involved  angular  stress  intensity 
factors  k.  The  convention  used  is  given  in  Figure  A-6.  It  was  suggested 
[A.6|  that  at  crack  branching  initiation,  the  angular  stress  intensity  factors 
are  given  by 


k  =  -  cos ( 9/2 ) 
1  2 


Kj.  ( 1  +  cos 6  )  -  3K.JJ  sine 


(A-9) 


kI;r  =  —  cos(e/2)  Kj  sine  +  (3  cose  -  l) 


Using  the  values  for  Kj  and  Kjj  obtained  from  the  displacement 
extrapolation  method,  kj  and  kjj  can  be  calculated  from  Eqn  (A-9). 
Replacement  of  K  by  k  in  Eqn.  (A-6)  gives  ' G '  for  values  of  crack  branching 
angle,  (A. 7).  This  method  has  the  advantage  that  only  a  single  run  of  the 
finite  element  model  is  require^ 

In  order  to  justify  the  use  of  the  displacement  extrapolation/angular 
stress  intensity  method,  the  results  of  a  comparison  between  that  method  and 
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FIGURE  A-5.  Extension  by  Node  Advance 
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FIGURE  A-6.  Definition  of  Angular  Stress 
Intensity  Factors 
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the  differential  crack  advance  method,  for  local  anisotropy  rotated  by  30°  and 
0°,  is  shown  in  Table  A-5.  The  results  are  plotted  in  Figure  A-7.  The  sets 
of  curves  for  1 G '  versus  branching  angle  for  the  models  considered  are  given 
in  Figures  A-8  and  A-9. 

A.1.6  Angled  Edge  Crack  Models 

As  an  introduction  to  the  modelling  of  angled  edge  cracks,  results 
previously  calculated  by  the  method  of  displacement  extrapolation  for  an 
isotropic  material  were  used  to  determine  the  strain  energy  release  rate  for  a 
range  of  crack  branching  angles.  For  cracks  angled  at  30°  and  45°  to  the 
direction  perpendicular  to  that  of  loading,  the  maximum  values  of  'G'  were 
obtained  for  crack  branching  angles  of  -30“  and  -40“ .  Using  the  values  of  Kj 
and  Kjj  obtained  by  the  method  of  Hellen  and  Blackburn  [A. 5]  intensity  factor 
method  gave  a  crack  branching  angle  of  -44°  for  the  45“  edge  crack  model. 

A. 1.7  discussion  on  2-01 wens ional  Analyses 

As  shown  in  Figure  A-7,  the  angular  stress  intensity  factor  method  gave 
results  very  similar  to  those  obtained  by  the  method  of  differential  crack 
advance.  The  former  was  favored  because  of  the  advantage  of  requiring  only  a 
single  program  run  for  each  model. 

Rotation  of  local  anisotropy  and  complete  model  anisotropy  gave  maximum 
strain  energy  release  rates  at,  or  within  4  degrees  of  the  col  inear  crack 
direction  using  LEFM.  On  the  basis  that  for  homogeneous  isotropic  materials 
crack  propagation  is  suggested  to  occur  in  the  direction  of  maximum  ‘ G ' ,  the 
present  work  indicates  that  crystallographic  mixed  mode  fracture  in 
anisotropic  regions  cannot  be  predicted  by  energy  considerations.  However, 
the  method  developed  is  capable  of  determining  the  crack  tip  stress  intensity 


factors  Kj,  Kjj  and,  for  30  cases  Kjjj  and  the  strain  energy  release  rate  for 
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FIGURE  A-7.  Strain  energy  release  rate  from  a)  differential  crack  advance  [deal  b)  angular 


local  aniso.  rotated  0 


whole  aniso.  rotated 


FIGURE  A-9.  Strain  Energy  Release  Rate  by  Angular  S.I.F.  Method 


a  crack  extension  on  a  particular  crystallographic  plane.  Although  K  is 
preferred  for  the  correlation  of  fatigue  crack  propagation,  questions  arise  as 
to  Its  suitability  and  consideration  should  be  given  to  local  crack  tip 
sensitivity  during  mixed  mode  fracture,  the  involvement  of  local  shear  strain 
and  whether  the  possible  significant  crack  tip  stress  system  affects  the 
assumed  material  homogeneity.  Use  of  Kjj  or  Kjjj  as  a  correlating  parameter 
for  shear  type  fracture  may,  however,  be  possible. 

A. 2  Three-Dimensional  Stress  Intensity  Factors  for  Cracks  in  Locally 

Anisotropic  Material 

This  section  describes  the  method  used  for  the  calculation  of  crack 
front  Stress  Intensity  Factors  in  3-dimensional  locally  anisotropic 
materials.  Two  finite  element  models  have  been  developed,  an  angled  edge 
through-crack  and  a  semicircular  surface  crack  in  uniaxial ly  loaded 
rectangular  prismatic  blocks.  The  models  were  considered  Isotropic  for 
comparison  with  published  results  and  then  modified  so  that  the  crack  was 
contained  within  a  region  of  local  anisotropy  to  simulate  the  structure 
sensitivity  experienced  by  small  cracks.  The  method  of  calculation  described 
can  be  applied  to  the  edge  and  semicircular  surface  cracks  for  any  orientation 
of  crack  and  material  anisotropy. 

A. 2.1  Three-Dluienslonal  Finite  Element  Models 

Two  finite  element  models  were  constructed  using  20  node  isoparametric 
brick  elements.  The  crack  fronts  were  represented  by  20  node  isoparametric 
brick  elements  with  one  face  degenerated  to  create  lb  node  wedges.  The 
elastic  singularity  at  the  crack  front  was  represented  by  movement  of  the 
midside  nodes  closest  to  the  crack  front  to  the  quarter  point  position, 
Figure  A-10.  The  quarter  point  node  coordinates  were  computer  generated  using 
FORTRAN  Code  written  specifically  for  this  work.  Both  models  were  developed 
using  PATRAN,  a  finite  element  pre/post  processing  package. 
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The  linear  elastic  models  were  uniaxial ly  loaded  by  a  fixed  uniform 
vertical  displacement  of  the  top  face  nodes,  zero  vertical  displacement  of  the 
bottom  face  nodes  and  a  limited  number  of  horizontal  nodal  zero  displacements. 
A. 2. 2  30  Degree  Angled  Edge  Through-Crack  Model 

Figure  A-ll  shows  the  48  finite  element  mesh  of  the  model  and  an 
enlarged  view  of  the  region  of  local  anisotropy.  The  crack  was  angled  at 
30  degrees  to  the  horizontal  X'  axis.  The  relative  dimensions  were  chosen  to 
allow  comparison  between  the  calculated  and  published  stress  Intensity  results 
for  the  completely  isotropic  case. 

A. 2. 3  Perpendicular  Semicircular  Surface  Crack 

Figure  A-12  shows  the  404  finite  element  mesh  of  the  model  and  an 
enlarged  view  of  the  crack  In  relation  to  the  bottom  half  of  the  region  of 
local  anisotropy.  The  angle  of  the  crack  and  the  relative  dimensions  of  the 
model  were  chosen  to  allow  comparison  between  the  calculated  and  published 
stress  intensity  results  for  the  completely  isotropic  case. 

A. 2. 4  Elastic  Properties  Of  Models 

As  with  the  2-dimensional  models,  the  anisotropic  ’’■egions  of  the  3- 
dimensional  models  were  represented  by  SRR99,  a  single  crystal  nickel  base 
turbine  blade  material.  The  isotropic  material  was  represented  by  MAR-M-002, 
an  equiaxed  nickel  base  turbine  blade  alloy  of  a  similar  composition  to 
SRR99.  The  elastic  properties  at  20°C  of  the  two  materials  used  in  the  models 
are  given  in  Table  A-l. 

The  finite  element  models  were  run  using  the  MARC  program,  level  Kl.3. 


FIGURE  A-12.  Surface  Crack  Mesh 


A. 2. 5  Crack  Orientation  Within  Anisotropic  Region 

The  models  with  local  anisotropy  were  designed  to  consider  a  crack 
along  a  locally  oriented  (111)  crystallographic  plane.  The  direction  of  crack 
advance  perpendicular  to  the  global  Y 1 Z *  plane  was  taken  to  be  the  locally 
oriented  (iTu).  Figures  A-13a  and  A-13L  show  a  general  and  plan  view 
respectively  of  the  semicircular  crack  in  relation  to  the  low  index  axes  of 
local  anisotropy.  An  identical  relationship  was  assumed  for  the  angled  edge 
through-crack  model  where  the  crack  front  was  parallel  to  the  global  Y '  Z ' 
plane  and  perpendicular  to  the  local  (110)  direction. 

Figure  A-14  shows  how  the  orientation  of  a  semicircular  surface  crack 
was  defined  in  terms  of  the  global  coordinate  axes  of  the  model.  The 
orientation  defined  in  Figure  A-14  can  be  applied  to  the  edge  through-crack 
with  the  angle  5  =  90°  across  the  crack  front. 

A. 2. 6  Elastic  Properties  Of  Rotated  Anisotropic  Material 

MARC,  the  finite  element  package  used  for  the  analysis,  has  the 
facility  to  orient  the  elastic  stiffness  matrix.  However,  because  the 
calculations  of  the  stress  intensity  farters  rpninred  a  knowledge  of  the  local 
elastic  compliance  along  the  crack  front,  the  rotation  of  anisotropy  was 
performed  by  computer  programs  written  for  this  work. 

To  obtain  the  global  elastic  properties  of  a  region  of  anisotropic 
material  with  known  local  properties,  the  transformation  described  in  Section 
A. 1.1  was  performed.  The  conventions  observed  for  the  rotating  coordinate 
systems  are  given  in  Figure  A- 15. 


A-32 


plan  of  (111) 
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A. 2. 7  Three-Dimensional  Elastic  Anisotropy 


Provided  attention  is  restricted  to  a  sufficiently  small  region  around 
the  crack  front,  problems  involving  elastic  anisotropy  can  oe  reduced  to  a  sum 
of  the  plane  and  antiplane  solutions.  The  method  used  for  3-dimensional 
elastic  anisotropy  was  derived  from  the  solution  to  plane  problems  described 
by  Sih  et  al  [A. 4].  Referring  to  the  notation  given  in  Figure  A- 16,  the 
analysis  of  the  plane  problem  assumed  the  r  axis  to  be  parallel  to  the  leading 
edge  of  the  crack  and  the  p  and  q  axes  to  be  parallel  and  normal  to  the  crack 
surface  with  the  pq  plane  being  a  plane  of  elastic  symmetry.  For  the 
semicircular  crack,  the  axis  p  was  an  extension  of  the  radius  of  the 
semicircle  and  q  was  the  (111)  direction. 

Plane  Elastic  Anisotropy 

The  stresses  and  displacements  close  to  the  crack  front  for  plane 
symmetric  and  skew- symmetric  loadings  are  given  in  Section  A. 1.3.  Figure  A- 17 
shows  the  notation  used  for  plane  extension. 

Antiplane  Elastic  Anisotropy 

The  analytic  functions  required  for  the  expressions  describing  stresses 
and  displacements  require  parameters  which  are  conjugate  pairs  of  roots  of  the 
following  characteristic  equation. 


a55u 


2 


44 


=  0 


The  expressions  for  stress  and  displacement  in  a  small  region  around 
the  crack  tip  for  the  shear  loading  are  as  follows,  where  r  is  a  radial 
distance  from  the  crack  front  rotated  at  an  angle  e  anticlockwise  from  the 
coplanar  crack  direction. 
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FIGURE  A-17.  Notation  for  Plane  Extension 
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Antiplane  Shear  Loadli 


/(2irr)  </(cos6  +  u3sin0) 


/  (  2 -nr ) 
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Displacement  parallel  to  the  crack  front  Is  given  by 

w  =  Km  /(2r/ir)  /(a44a55  -  a452)  IM[/(Cose  +  u3sine ) ]  (A-10) 

A. 2.8  Displacement  Extrapolation  For  Calculation  of  KI,  KII  and  Kill 

For  a  small  region  surrounding  the  crack  tip,  the  expressions  for 
displacements  parallel  and  perpendicular  to  the  crack  face  and  parallel  to  the 
crack  front  have  been  given,  Eqns.  (A-4.A-5.A-10) .  Putting  9  =  180  degrees 
such  that  1 r 1  lies  along  the  crack  face  modifies  those  expressions  as  follows 

Plane  Symmetric  Loading 


u  =  ( 2r/m  )  Re 


U1  "  u2 


(»l*2  "  y2pl) 


v  =  K^/ (  2r/ir )  Re 


y  1  '  y2 


(v1q2  -  p2qi) 


(A— 11) 


Plane  Skew  Symmetric  Loadli 


u  =  2r/tr )  Re 


v  =  K^/(2r/it)  Re 


y  1  '  w2 


(P2  -  PX) 


-  v2 


(A- 12 ) 


pj  =  Slluj  +  S12  "  S16wj 
qj  =  S12yj  +  (S22/uj)  '  S26 


Antiplane  Shear  Loading 

w  =  Km  /2r/w)  ^(a44a65  -  a452) 

Surominy  the  corresponding  displacements  above  for  the  symmetric  and 
skew  symmetric  cases  gives  the  displacements  for  the  mixed  mode  plane 
condition.  Equations  (A- 11)  to  (A- 13)  can  be  reduced  to  the  isotropic  case  by 
substitution  of  the  following  values. 
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(14) 
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(16) 

Calculation  of  the  relative  movement  of  opposite  pairs  of  crack  face 
nodes  parallel  to  the  normal,  the  binormal  and  the  tangent  to  the  crack  front 
(pqr  in  Figure  A- 16)  and  their  substitution  into  the  displacement  Equations 
(A- 11)  to  (A- 13)  will  give  values  of  Kj,  Kjj  and  Kjjj  at  a  distance  ' r '  from 
the  crack  tip.  Extrapolation  of  these  values  to  r  =  0  will  give  the  crack  tip 
stress  intensity  factors  K. 

A. 2. 9  Elastic  Strain  Energy  Method 

The  elastic  stain  energy  difference  was  calculated  in  MARC  by  a 
differential  stiffness  method  [ A. 3 |  and  the  Strain  Energy  Release  Rate  1 G 1  was 
calculated  by  dividing  tne  strain  energy  difference  by  the  area  of  crack 
advance.  The  mesh  distortion  associated  with  the  calculation  of  G  is  shown  in 
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Figure  A-18.  All  the  values  of  G  obtained  by  the  Strain  Energy  Method  assumed 
advance  of  the  crack  in  Its  existing  plane.  The  results  were  used  to  verify 
the  Kj,  Kir  and  Km  values  calculated  by  the  Displacement  Extrapolation 
(Section  A. 2. 8)  as  follows: 

Verification  Of  Calculated  K  Values 

The  method  of  differential  stiffness  for  the  calculation  of  G  leads  to 
problems  In  the  simultaneous  solution  for  Kj,  K j j  and  K|j|.  However,  the 
values  of  G  associated  with  each  of  the  three  modes  of  opening  can  be 
calculated  as  follows  I  A. 4] 


KI 

GI  =  -  -  S22  Im 


VII 


KI^1  +  V  *  KII 
PlV2 


II  “  ~T  S11  Im[  KII{yl  +  y2)  +  KIuly2  1 


G  -  K][I1 

Gm - 


a44a55 


2  /(a44a55  "  a45  ) 


( A- 17 ) 


(A-18) 


(A— 19 ) 


For  coplanar  crack  advance  and  general  anisotropic  conditions  the  value  of  G 
is  given  as 

G  =  Gj  +  G  j  j  +  G  j  j  j  (A-20) 

Substitution  of  the  values  of  Kj,  Kjj  and  Kjjj  from  displacement 
extrapolation  into  the  above  equations  allows  comparison  between  a  directly 
calculated  and  a  derived  value  for  G  and  hence  a  means  of  verifying  the  values 
of  K. 

In  addition,  for  cracking  on  the  (111)  plane,  a  plane  of  elastic 
symmetry  as  specified  in  the  derivation  of  the  characteristic  Equation  (A-4) 
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FIGURE  A-I8(a).  Differential  stiffness  energy  release 
calculation  -  mesh  distortion. 


FIGURE  A- 18(b) .  Energy  release  rate  calculation,  3D 
crack  front. 
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only  exists  at  two  points  around  the  semicircular  crack  front.  The  comparison 
between  the  calculated  and  derived  values  of  6  allows  assessment  of  the  extent 
to  which  the  plane  assumption  must  be  rigidly  applied. 

A. 2. 10  Results  of  3-01»ens1onal  Analyses 

The  computer  coding  developed  for  the  generation  of  elastic 
properties,  the  calculation  of  the  roots  of  the  characteristic  equations  and 
the  calculation  of  stress  intensity  factors  from  the  displacement 

extrapolation  method  was  used  for  both  the  isotropic  and  locally  anisotropic 

models. 

30  Degree  Angled  Edge  Through-Crack 

The  48  element  FE  model  used  for  this  analysis  was  set  up  as  a 

preliminary  trial  to  investigate  the  use  of  MARC  with  3-dimensional  locally 
anisotropic  crack  models.  As  a  consequence  of  the  coarseness  of  the  mesh  and 
the  resulting  scatter  of  the  rear-tip  nodal  displacements,  the  stress 
intensity  results  for  the  isotropic  case  were  lower  than  the  published  results 
for  angled  edge  cracks  [ A . 1 ,  A . 8 1  by  at  least  7  percent.  Having  developed  the 
techniques  for  handling  such  problems,  efforts  to  refine  and  rerun  the  job 
were  diverted  in  preference  of  the  semicircular  surface  crack  model. 

The  stress  intensity  factors  for  the  angled  edge  crack  were 

normalized  with  respect  to 

Kq  =  o/(ira  ) 

where  o  =  remote  uniaxial  stress 

a*  =  crack  length 
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Figure  A-19  shows  the  normalized  Mode  I  stress  intensity  factors,  for 
the  isotropic  and  locally  anisotropic  models  and  the  isotropic  plain  strain 
results  given  by  the  Rooke  and  Cartwright  [A.l]  and  Preston  and  Pickard  (A.8|. 

Figure  A-20  shows  the  normalized  Mode  II  and  Mode  III  stress 
intensity  factors  for  the  isotropic  and  locally  anisotropic  models.  The  plane 
strain  Mode  II  results  according  to  Rooke  and  Cartwright  [A.l]  and  Preston  and 
Pickard  [ A. 8 ]  are  also  shown. 

Figure  A-21  shows  G^g/Gg,  the  ratio  of  strain  energy  release  rates 
derived  by  displacement  extrapolation  and  calculated  by  the  energy  method. 
The  lack  of  symmetry  of  G^g/Gg  along  the  crack  front  of  the  anisotropic  model 
indicated  that  the  displacement  extrapolation  method  and  the  strain  energy 
method  were  unable  to  give  equivalent  results  with  the  coarse  mesh  provided. 

Semicircular  Surface  Crack 

The  stress  intensity  factors  were  normalized  with  respect  to 
Kq  =  2  o/( a/w) 

where  o  =  remote  uniaxial  stress 
a  =  crack  radius 

Although  the  method  of  analysis  was  applicable  to  angled  surface 
cracks,  it  was  decided  to  validate  the  model  for  the  isotropic  case  by 
comparison  with  published  results  (A.l,  A. 9,  A. 10)  and  so  a  perpendicular 
surface  crack  model  was  developed. 

Figure  A-22  shows  the  normalized  Mode  I  stress  intensity  factors 
for  the  isotropic  model  around  the  crack  front  for  6  =  0  to  90°C  compared  with 
the  full  curve  by  Rooke  and  Cartwright  (A.l)  and  surface  and  full  crack  depth 
results  by  Hellen  and  Blackburn  [A. 9)  and  Pickard  [A. 10).  According  to  Hellen 
and  Blackburn  [A. 9],  the  results  at  a  free  surface  should  be  treated  as  under 
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(14,22)  isotropic 


de  displacement  extrapolation 


FIGURE  A-21.  30  Degree  Angled  Edge  Through-Crack 
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FIGURE  A-22.  Semicircular  Surface  Crack 


plane  stress  and  their  result  and  that  given  by  the  current  analysis  are  shown 
in  Figure  A-22. 

Figure  A-23  shows  the  isotropic  and  locally  anisotropic  normalized 
Mode  I  stress-intensity  factors  given  by  the  current  analysis  for  6  =  0  to 
180°. 

The  normalized  K j j  and  Kjjj  values  for  the  locally  anisotropic  case 
are  given  in  Figure  A-24.  For  the  isotropic  case  the  K  j  j  and  Kjjj  values  were 
zero.  The  Mode  II  anisotropic  plane  stress  solution  at  the  free  surface  is 
also  given  in  Figure  A-24. 

Figure  A-25  shows  Gjg/Gg,  the  ratio  of  strain  energy  release  rates 
derived  by  displacement  extrapolation  and  calculated  by  the  energy  method. 
Although  it  is  recognized  that  the  energy  method  provides  a  weighted  average 
value  for  the  strain  energy  release  rate  close  to  the  free  surface,  ..hereas 
that  derived  from  the  displacement  extrapolation  K  values  is  specific  to  the 
surface,  the  two  values  of  G  were  much  closer  when  plane  stress  conditions 
were  assumed  for  the  Modes  I  and  II  solutions,  both  for  the  isotropic  and  the 
locally  anisotropic  cases. 

A. 2. II  Conclusions 

The  displacement  extrapolation  technique  was  used  to  determine  3- 
dimensional  stress  intensity  factors  in  locally  anisotropic  materials.  This 
technique  has  the  advantage  of  being  able  to  separate  the  values  of  K  for  each 
of  the  three  modes  of  opening  and  when  used  with  quarter  point  elements  at  the 
crack  front,  is  expected  to  yield  accurate  results.  However,  because  of  the 
material  anisotropy  the  3-dimensional  solutions  should  trictly  only  be  valid 
at  points  where  a  plane  of  elastic  symmetry  exists  simultaneously 
perpendicular  to  the  crack  plane  and  crack  front.  In  order  to  check  the 
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FIGURE  A -23.  Semicircular  Surface  Crack 
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FIGURE  A- 25.  Semicircular  Surface  Crack 


rigidity  of  this  symmetry  restriction,  the  elastic  strain  energy  release  rate 
G,  derived  from  the  K  values,  was  compared  with  the  G  calculated  by  the 
differential  stiffness  (energy)  method.  The  values  of  Gde/Ge  given  for  the 
perpendicular  semicircular  surface  crack  model  show  that  the  symmetry 
requirement  does  not  need  to  be  enforced  rigidly.  Confidence  in  the  locally 
anisotropic  K  values  being  consistently  higher  than  the  corresponding 
isotropic  values  leads  to  interesting  consequences  for  the  presentation  of 
'small'  crack  growth  data  in  relation  to  those  for  'long'  cracks. 

Although  the  effect  of  plane  stress  conditions  at  the  free  surface 
would  be  very  localized,  the  disparity  between  the  displacement  extrapolation 
and  the  differential  stiffness  (energy)  methods  was  reduced  by  such  an 
assumption.  Together  with  metallurgical  features  which  hinder  crack  advance 
along  the  free  surface,  the  plane  stress  condition  could  be  a  contributor, 
factor  in  the  successful  modelling  of  crack  shapes.  Accurate  component  iife 
prediction  depends  very  much  on  knowledge  of  the  crack  shape  and  accurate 
stress  intensity  solutions.  The  results  given  here  indicate  that  the  stress 

intensity  solutions  used  to  present  small  crack  data  may  require  revision. 

A. 3  Suronary 

This  report  considers  the  calculation  of  stress  intensity  factors  in 
a  material  containing  a  small  crack  in  a  region  considered  to  be  locally 
anisotropic.  The  finite  element  analysis  work  was  developed  in  two  stages,  a) 
2D  plane  solutions  and  b)  3D  solutions.  Both  stages  involved  the  calculation 
of  isotropic  solutions  to  check  the  models  and  methods  against  published 
results  before  the  introduction  of  local  anisotropy.  All  the  analyses  assumed 
linear  elasticity. 
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20  Finite  Element  Analysis 


A  perpendicular  edge  crack  model  was  constructed  using  8  node 
isoparametric  quadrilateral  plane  strain  elements.  The  cre_k  tip  was 
constructed  using  degenerate  crack  tip  elements  to  represent  the  elastic 
singularity.  Stress  intensity  factors  were  calculate  '  by  both  the  strain 
energy  method  and  the  displacement  extrapolation  method. 

The  modelling  of  crystal lographic  crack  oranching  was  investigated 
for  a  range  of  orientations  of  local  anisotropy.  The  strain  energy  release 
rate  peaked  at  or  very  close  to  the  direction  of  col  inear  crack  advance  and  it 
was  concluded  that  this  parameter  could  not  be  relied  upon  to  predict  the 
direction  of  crystallographic  crack  growth. 

3D  Finite  Element  Analysis 

Two  models,  an  angled  edge  through  crack  and  a  perpendicular 
semicircular  surface  crack,  were  constructed  using  20  node  isoparametric  brick 
elements.  The  crack  tip  region  was  constructed  using  degenerate  quarter  point 
elements  to  represent  the  elastic  singularity.  Strain  energy  calculations  for 
the  locally  anisotropic  model  were  used  to  verify  the  displacement 
extrapolation  results.  The  anisotropic  region  around  the  crack  was  orientated 
such  that  the  crack  ran  along  a  close  packed  plane. 

The  results  of  the  analyses,  for  a  limited  number  of  anisotropic 
orientation  and  geometric  configurations,  have  revealed  that  increased  stress 
intensities  (on  average,  5  percent  greater  than  the  isotropic  case.  Figure  A- 
23)  may  be  observed  for  particular  configurations  compared  to  those  derived 
using  isotropic,  homogeneous  models. 
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